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Final Report - Rapid Structural Analysis Methods from Open VSP Geometry 


Objective - Develop and validate advanced/rapid structural analysis methods from existing 
Open Vehicle Sketch Pad (VSP) analysis tools as a continuation of work performed under NIA 
IDIQ Contract NNL08AA08B Task Order NNL11AC44T. 

Work Performed - Rapid structural analysis for use in conceptual design of advanced concepts 
along with tool validations 

Statement of Work - 
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3.1 Developed rapid structural analysis methods from OpenVSP geometry. 

3.1. Develop a multi-section wing structural design and analysis method implemented in an open- 

source software code that extends the current OpenVSP and VSP Structural Analysis Module 
(SAM) SAM capabilities. This task shall include description of the theory applied, validation of 
the implementation, source code, and a user’s guide. 

3.2. Research a fuselage structural design and analysis method for the “Fuse2” component 
implemented in an open-source software code that extends the current OpenVSP and VSP SAM 
capabilities. This task shall include description of the theory applied. 

3.3. Develop a set of load cases that typically are active design constraints for the wing and fuselage 
of an airliner and demonstrate an automated design of the wing and fuselage structures with 
these cases. Document the regulation that the case is derived from, the implementation strategy, 
the application on a test case, and a means for identifying which constraint cases are active for 
each element of the Finite Element Model (FEM). 


Accomplishments 

3.1 MSWing - Task performed, theory developed and applied. Validation was targeted at 
Boeing/McDonnell Douglas Advanced Composites Technology Test Wing. Results 
documented in attached Briefing. 

VSP related grid quality issues and absence of multi-section wing unique notation impeded 
progress. MSWing notation issues resolved by writing translator that sequences skin, rib 
and spar notation to achieve element and node uniqueness required. The grid quality issue 
ws not resolved we were unable to implement the rigid body body connections required to 
continue on to a CalculiX solution. Work on the task continues and we will share the results 
upon completion. User's guide updated and attached. 

3.2 Fuselage research and structural design - task was removed from the contract 

3.3 Recommended loads cases cases developed and documented as follows: 

Design Load Conditions are limited to wings only with a focus on methodology development. 
For simplicity commercial transport assessments have been done at maximum takeoff weight 
(MTOW) and an in-flight symmetrical Design Ultimate Load (DUL) = 3.75 g or 1.5 times Design 
Limit Load (DLL) = 2.5 g. For future CD stress and/or mass prediction applications we 
recommend use of what we consider a more representative subset of design air and ground 
loads. For commercial transport applications our recommendations are based on Federal 
Aviation Regulation Part 25, Subpart C - Structure. For general aviation and larger Unmanned 
Air Systems (UAS) we recommend the Part 23 equivalent. For military aircraft we reccommend 
loads based on service standards and vehicle type. For mass estimate calibration applications 
we use actual vehicle handbook values. For example, when we evaluated the A6-E metal wing 
we used design loads from the USN A6-E Detail Weight and Balance Report. 
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Specific recommendations for air transport CD application include the following subset of Part 
25 requirements: 

a. Design factor of safety = 1.5: FAR 25.303 

b. Design limit load assessments based on yield stress - FAR 25.203(a) 

c. Design ultimate load assessments based on ultimate stress - FAR 25.203(b) 

d. Symmetric maneuvers with zero pitching acceleration - FAR 25.331(b) 

e. Flight maneuvering envelope and design airspeeds defined at four points as shown in 
Figure 3 for flaps up only: FAR 25.333, 25.335(a) and 25.335(c) 

f. Limit maneuvering load factors = +2.5 g and -1.0 g - FAR 25.337 

g. Flight design weights at 2.5g at MTOW and 2.25 g at zero fuel weight: FAR 25.343 

h. Landing loads at a level landing condition with lateral drift: FAR 25.475(i) 

i. Braked roll loads at a level roll out at MTOW with a vertical reaction factor of 1.0 and a 
coefficient of friction of 0.8: FAR 25.493(a) 

We also recommend application of additional Part 25 load conditions at a TBD fuel tank over- 
pressure with and crash load-induced fuel ram effect. For example, the Boeing ACT wing was 
designed for 15 psi fuel over-pressure and a 9g crash load 4 . Wing tank pressure and fuel ram 
effects are currently not implemented as VSP SAM loads. They can, however, be defined as 
design loads through CalculiX. 

Deliverables : 

- Final results briefing: attached 

- Source code: by separate attachment 

/s/ 

Armand J. Chaput 
Principal Investigator 
15 June 2015 
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Task 3.1 - Develop Rapid 
Structural Analysis Methods 
from Open VSP Geometry 


15 June 2015 
Final Results Briefing 

Armand J. Chaput 
University of Texas at Austin 

Department of Aerospace Engineering 
and Engineering Mechanics, University 
of Texas at Austin 
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VSP Structural Analysis Module (SAM) Objectives 


(1) Expand VSP capabilities - make higher order, physics 
based tools and methods available for conceptual design 
(CD) application 

(2) Integrate VSP structures module with an open source 
finite element method (FEM) structural analysis 
program - provide user friendly interactive environment 

- Currently focused on Calculix (available under terms of GNU General 
Public License as published by the Free Software Foundation) 

(3) Develop open source methods for FEM-based CD 
mass property (MP) estimation 

- Prior effort developed basic single trapezoidal section wing 
and tail structural analysis and mass estimation methods 

• i.e., limited to linear thickness and constant taper ratio 

- Current effort expands capabilities to multi-section wings 
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Background - Airframe Mass Property (MP) Estimation 


Airframe mass is driven by multiple requirements, many of 
which are not captured by traditional parametric methods 

- A problem for mass prediction involving new vehicles and/or 
new materials, processes and/or design requirements 

Geometry, loads, and other design requirements are primary 
drivers for load carrying airframe structural mass 

- FEM analysis can predict required primary structural mass 
with good accuracy compared to parametric methods 

Non-primary structural mass driven by other requirements 

- Some of which can be captured by CD FEM-based methods 

- Some of which cannot, e.g., clips, brackets, fasteners, etc. 

- Therefore, parametric adjustment factors are still required 

Even with factors FEM-based methods improve CD primary 
structure mass estimate quality vs. parametric methods 

- It doesn’t cost more money nor take any more time 
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VSP SAM Computational Architecture 


Vehicle Sketch Pad 

Parametric 
Internal Geometry 


Parametric 
External Geometry 


External and Internal Mesh 
Generation 


UT Input Executable (Java) 

Wing Trim 

Thickness and 
Material Properties 

Boundary Conditions 
and Load Cases 


CalculiX Input File 


UT Convergence Executable 
(Java) 


Thickness Iteration 


Solution Files 


Stress Convergence 

Mass Calculation 




CalculiX© 
FEM Input 

FEM Solution 


FEM Post Process 
and Graphics 

Output Files 
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Current VSP SAM CD Load Cases 



Gust envelope 



Positive High Alpha 



Positive Low Alpha 



Negative Low Alpha 




Nodes sized to 
most demanding 
load case 



Negative High Alpha 
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Proposed Load Cases for Conceptual Transport Aircraft 


Based on selected requirements from FAR Part 25, Subpart C - Structure 

• Design factor of safety = 1 .5: FAR 25.303 

• Design limit load assessments based on yield stress - FAR 25.203(a) 

• Design ultimate load assessments based on ultimate stress - 
FAR25. 203(b) 

• Symmetric maneuvers with zero pitching acceleration - FAR 25.331 (b) 

• Flight maneuvering envelope and design airspeeds defined at four 
points each at minimum controllable airspeed and maximum level flight 
speed at maximum positive and positive load conditions : FAR 25.333, 
25.335(a) and 25.335(c) 

• Limit maneuvering load factors = +2.5 g and -1 .0 g - FAR 25.337 

• Flight design weights at 2.5g at MTOW and 2.25 g at zero fuel weight: 

FAR 25.343 

• Landing loads at a level landing condition with lateral drift: FAR 
25.475(i) 

• Braked roll loads at a level roll out at MTOW with a vertical reaction 
factor of 1 .0 and a coefficient of friction of 0.8: FAR 25.493(a) 

Recommend consideration of additional fuel tank over-pressure with a 
crash load-induced fuel ram effect at 15 psi and 9g’s axial 
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Application to Boeing MDC Advanced Composite (ACT) Wing 



Load Case(s) 
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Advanced Composite Technology (ACT) Wing 


Wing carry- 
through not 
included in 
ACT test box 


ACT thickness 
not linear 


Mass analysis based on data from NASA 
CR-2001-210650-AST Composite Wing 
Program-Executive Summary 


Sources: multiple NASA and Boeing ACT Documents 



Mass Categories 

Ibm 

Fraction 

Upper and lower skins 

1656 

0.41 

Spar caps and stringers 

1238 

0.31 

Spar webs 

350 

0.09 

Stress based (subtotal) 

3244 

0.81 




Aero ribs and intercostals 

520 

0.13 

MLG rib blkhd 

75 


MLG pad up 

31 

0.03 

Bolts and nuts 

80 

0.02 

SOB Pad up 

25 


Access panel pad-up 

45 

0.02 

Total 

4020 

1.00 


Process - GrEP VARTM 


Spar web cr CDN - 27.0 ksi 
Ribs <j CDN = 37.9 ksi 
Skin o CDN = 36.5 ksi 
Min gage (given) 

= 0.22 in 





Actuator 
load (Ibf) 

2y/b 

1 

40500 

0.947 

2 

99750 

1.000 

3 

-3000 

0.634 

4 

21000 

0.704 

5 

15000 

0.384 

6 

-45000 

0.382 

7 

45000 

0.202 

8 

6000 

0.296 
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Single Section Constraint 


A-6E Wing Chord Distribution 


1.20 

1.00 

0.80 



0.40 

0.20 

0.00 



BL (in) 


ACT Wing Box t/c 



A-6E Wing Thickness Profile 



VSP structural module grid limited to single 
trapezoidal planform and stream-wise chord 
- Even CD level structures need more 


flexible multi-panel capability 


n ACT 

CDNS(psi) 


Case 1 

Case 2 

Diff 

Skin 

36.46 

47.31 

30% 

Ribs 

37.85 

31.94 

-16% 

Spars 

27.02 

34.04 

26% 


- In order of priority we need (1) non-linear 
t/c, ( 2 ) non-linear chord, ( 3 ) non-stream- 
wise CUt Capabilities (2013 Chaput email) 
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Other Generic Boeing Transport Wings (hand iterated) 


Aircraft or 
article 

W DG [lbs] 

DUL 

N z 

Half 

span S ref 
[ft 2 ] 

Sweep 

[deg] 

Half 

span [ft] 

Half 

span 

AR 

TR 

(t/c) r 

(t/c), 

Dihed 

[deg] 

Wing 

Fuel 

Engine 
pylon or 
store (per 
side) 

Engine or 
store 
mass (ea) 

747-100 

713,000 

3.75 

2790 

37.5 

98.5 

6.96 

0.265 

0.1794 

0.078 

7 

Y 

2 

8608 

737-200 

100,800 

3.75 

502.5 

25 

45.4 

8.21 

0.220 

0.126 

0.112 

6 

Y 

1 


727-300 

160,000 

3.75 

793.5 

32 

55.2 

7.67 

0.265 

0.154 

0.09 

3 

Y 

0 




B747-100 VSP 
Simplified to 15 ribs 


i 


8.04e+006 
7.65eO06 
?.27e+006 
6.89e*006 
6 . 51e+006 
6 . 12e+006 
5 . 74e+006 
5 . 36e+006 
4.97e+006 
4 . 59e*006 
4 . 21e+006 
3 . 83e+006 
3.44e*006 
3 . 06e+006 
2 . 68e+006 
2 . 30e+006 
1.91e+006 
1.53e+006 
1 . 15e+006 
7 . 65e*005 
3.83e+005 


o cdn — 46.5 ksi 



B747, Min gage = 
0.1” exc. ribs = 0.36” 


NASA 110392 
Wing Weights 

PDCYL 

(Ibm) 

Load- 

carrying 

structure 

(Ibm) 

Primary 

structure 

(Ibm) 

Total 

structure 

(Ibm) 

B-727 

8688 

8791 

12388 

17860 I 

B-737 

5717 

5414 

7671 

10687 

B-747 

52950 

50395 

68761 

88202 

B-720 

13962 

11747 

18914 

23528 | 

DC-8 

22080 

19130 

27924 

35330 

MD-11 

33617 

35157 

47614 

62985 | 

MD-83 

6953 

8720 

11553 

15839 

L-1011 

25034 

28355 

36101 

46233 i 



B727-100 VSP 
Full 26 ribs 


4.57e*006 
4.33e*006 
4 . 09e+006 
3 . 85e+006 
3 . 60e+006 
3 . 36e+006 
3.12e*006 
2.88e+006 
2.64e*006 
2.40e+006 
2 . 16e+006 
1.92e*006 
1 . 68e+006 
i . 44e+006 
1.20e*006 
9 . 61e+005 
7 . 2ie+005 
4 . 81e*005 
2.40e*005 


Ocd/v = 27 - 4 ksi 



B727, Min 
gage = 0.1” 



B737-100 VSP 
Simplified to 14 ribs 


2.32e+006 
2 . 20e*006 
2.09e+006 
1.97e+006 
1.86e+006 
1.74e+006 
i . 62e+006 
1.51e+006 
1 . 39e*006 
1.28e+006 
1 . 16e+006 
1 . 04e+006 
9 . 28e+005 
8 . 12e+005 
6.96e+005 
5.80e+005 
4.64e+005 
3 . 48e+005 
2 . 32e+005 
1 . 16e+005 


Gqdn — 29.8 ksi 



B737, Min gage = 
0.1” exc ribs = 0.15” 


Wing Boxes Notional and Extrapolated to Centerline 


Ardema, M. D., Chambers, M. C., Patron, A. P., Hahn, A. S., Hirokazu M., Moore, M. D., “Analytical Fuselage and 
Wing Weight Estimation of Transport Aircraft”, NASA Ames Research Center, NASA Technical Memorandum 
110392, California, 1996 
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VSP SAM Mass Estimate Correlations 


Aircraft or 
article 

MTOW 

(Ibm) 

Sref 

(sqft) 

Wing 

loading 

(psf) 

Span 

(ft) 

AR 

LE Sweep 
(deg) 

TRl 

(nd) 

(t/c) r 

(t/c), 

DUL n z 

Wfdgw/Sref 

(psf) 

Wing 

Mounted 

Engine 

1/2 Wing 
box mass 
(Ibm) 

Ribs (Ibm) 

Skins 

(Ibm) 

Spars 

(Ibm) 

Weighted 

CDNS(psi) 

Rib 

CDNS 

(psi) 

Skin 

CDNS 

(psi) 

Spar 

CDNS 

(psi) 

B747-21P 

713000 

5469.0 ' 

130.4 

197.1 

6.96 

38 

0.26 

0.18 

0.08 

488.9 

JT9D-3A 

25483 

1221 

17238 

7023 

46.50 

46.50 

46.50 

46.50 

B737-200 

100800 

1005.0 

100.3 

90.8 

8.21 

25 

0.22 

0.13 

0.11 

376.1 

JT8D-7A 

2726 

280 

1717 

730 

29.80 

29.80 

29.80 

29.80 

B727-100 

160000 

1587.0 

100.8 

110.3 

7.67 

32 

0.26 

0.15 

0.09 

378.1 


4264 

527 

3007 

730 

27.40 

27.40 

27.40 

27.40 

A-6E 

60705 

529 

114.8 

53 

5.31 

29.5 

0.857 

0.90 

0.059 

673.2 

n/a 

1052 

161 

654 

237 

34.40 

12.36 

46.18 

16.88 

ACT Test 

181000 ’ 

1385 

130.7 

129.4 

12.1 ’ 

29* 

0.63** 

0.42** 

0.17** 

490.1 


3997 

582 

2633 

782 

42.48 

31.94 

47.31 

34.04 

X-56 

374 

56 

6.7 

28.0 

14.0 

22 

1 



31.7 


41.12 

3.41 

24.50 

13.20 

29.74 

9.00 

40.80 

14.60 


*based on MDX **based on test box 
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Concluding remarks - VSP Structural Design and Analysis 


1 . FEM design and analysis can be accommodated during 
conceptual design without adding onerous requirements for 
higher levels of design detail 

- Serious structural design issues can be assessed, identified and 
resolved without slowing down the pace of design 

2. Design and analysis methods can be applied intelligently by 
designers who are not structural specialists 

- Tool specific details have been pushed into the background 

3. FEM model results correlate with actual mass property data 
through use of mass-correlated o CDN 

4. VSP SAM MSWing capability impeded by V2 grid quality 
including at connections and absence of V3 structures module 

5. Otherwise VSP SAM methods are ready to replace parametric 
mass estimates for primary load carrying wing structure 
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Study Intent & Background 

► Design a vehicle concept that is VTOL 
and Cruise Efficient 

° Zero ground infrastructure 
° 24 hr mission 

► Performed Conceptual Design Study: 

° Mission seemed possible, but would require 
new technology and advanced vehicle 
configuration 

► Funded to build and flight test a 50% 
scale prototype 

° Purpose of funding is to: 

• Flight demonstration of robust inbound and outbound 
transitions 

• Flight demonstrate cruise efficiency (not full 24 hrs) 



201 1 - Early Concepts for Long Endurance VTOL 




201 2 - 35% Scale Prototyping and Flight Testing 



Dos Samara 


Trifecta (AKA Baseline) 


Splitwing 



201 3 - Greased Lightning (GL-1 0) 


• A new design direction was taken to take advantage of new distributed electric 
propulsion benefits 



Greased Lightning (GL 



General Arrangement 




General Arrange 










nboard Arrangement 


Variable Pitch 
Forward Flight Props 


Batteries housed 
in each nacelle 


Fixed Pitch 
Vertical Flight Props 
(shown folded) 



Retracted Landing Skids 


Qty 2 Diesel Engines 
8 hp each 

Turning Alternators 



Satellite Communication 
System 


Aeronautics Systems Analysis 
Branch, SACD, LaRC 


Aerodynamic Controls 




oE Model 


ranc 





Aeronautics Systems Analysis 
Branch, SACD, LaRC 





Control Strategy 


► Vertical Flight 
° Pitch 

• Fore/Aft Propeller Thrust Modulation 

° Roll 

• Left/Right Propeller Thrust Modulation 

° Yaw 

• CW/CCW Propeller Thrust Modulation 

• Flaperons 

► Forward Flight 
° Pitch 

• Fast -> Elevator 

• Slow -> Tail Rotation 

° Roll 

• Flaperons 

° Yaw 

• Rudder 




Aeronautics Systems Analysis 
Branch, SACD, LaRC 


GLARF 















Final Geometry 


Pressure (Pa) 

_ 1000.0 

600.00 

200.00 



Aeronautics Systems Analysis 
Branch, SACD, LaRC 


4 Flying Models 

► Pilot Proficiency Trainer 

Purpose: Provide a tilt wing aircraft that can hover, transition, and fly in wing borne flight IOT train the pilot with an 
aircraft that transitions. 

° Weight: less than 8 lbs 

° Endurance: greater than 3 minutes of hover time 
° Wingspan: TBD 

° Control System: None, but will use rate gryos IOT make easier to fly 

► GLARF (Greased Lightning Almost-Ready-to-Fly) 

Purpose: 1) Test different experimental control law designs IOT reduce risk for Greased Lightning flight testing. 2) 
Pilot Proficiency 

° Weight: ~1 2 lbs IOT test on the zip line. 

° Endurance: greater than 3 minutes of hover time 
° Wingspan: ~6 ft 

° Control System: Partial ASROV (2 lbs) 

► Hovering Plywood Bird (aka Iron Bird) 

° Purpose: 1 ) Validate performance of wiring harness. 2) Validate ASROV hardware in hovering only flight. 

° Weight: 35 to 55 lbs 

° Endurance: greater than 3 minutes of hover time 
° Wingspan: ~1 0.5 ft 
° Control System: Full ASROV (3.5 lbs) 


Greased Lightning 

Purpose: Provide a flight demonstration of a VTOL long endurance concept IOT gather interest and support for future 
work. 

° Weight: 55 lbs 

Endurance: Goal -> 80 minutes 
J/Vingspan: ~1 0.5 ft 

tern: Full ASROV (3.5 lbs) 


Aeronautics Systems Analysis 
Branch, SACD, LaRC 


Control System Development 


All 4 of the models 

Each of these models were to have a custom designed ASROV Control system 
installed in each one. This system was designed by the Dynamic Systems and 
Controls Branch. This was designed to be a high intensity and very in depth 
controller unit. The board had custom designed components just for this 
project. The time that we were working on this custom unit, I spent time 
learning COTS units that could achieve the same end goal of flying this 
vehicle concept through a transition. 

I started to focus my attention on multiple different paths for a COTS unit to 
accomplish this task. In all my years in UAV’s a good place to start for 
answers is RC forums online, and talking to people in the industry. I found 
that there are multiple VTOL UAV controllers, but they only do hover and 
forward flight, not transition from one to the other. 

I then tried starting a few different VTOL control boards that could possibly 
hp a rp<;parrh tnnl fnr thp A’sROV engineers. 



ASROV Control Unit 



APM 2.5 / PIXHAWK Controller Board I 


I started doing hover only flights with the APM 2.5 on the 
GLARF (251b) model as shown in the picture. This model 
hovered great with it APM, but no capability for forward 
flight transition. 

Then proceeded to investigate using the APM for 
transition work on a quick 2 day build foam airplane. I 
learned very quickly the issues in the code for using it as 
a transition control board. 

The open source software gives you 
flexibility to program anything, but 
it takes too long to modify the code. 

The pixhawk board was now being 
installed as a data acquisition board only, 
not to control the model in any fashion. 

The pixhawk has a much higher frequency 
data acquisition rate, that’s why we used it 
only for data collecting. 


yJPPUPlLQT §: 
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KK Flight Controller Board 


During the research stage of finding different control boards to use I came across 
the KK VTOL Mixer Control Board. I found this board on a RC Forum a man that 
was designing/building VTOL UAV’s that transition from hover to forward flight. 
He designed a model that is similar to GL-1 0 called the XC-1 42, and the full scale 
model of this was built and flown in the 1 960’s. I knew that this transitional 
control board would work for our VTOL design, and proceeded ahead testing it. 
The board started being tested on a 21b foamie all the way up to the 551b GL-1 0 
model, and have had success all the way through!!! 



KK Control Board 




Subscale XC142 with KK Board 
during transition period 






GLARF / GL-1 0 Flight Test Cards 


CARD MANEUVER 


Hover Directional Stability 
SFF Pitch Trim & Authority 
Full Transition / Tail Rotation Schedule 
Fwd Flight Stalls 

Roll/Yaw Oscillation Investigation 
Higher Weights 

Pitch and Directional Stability w/ Reducing 
Tail sizes 


te: Flight test cards are written for specific vehicles, but follow these descriptions 
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Turbulence Model 
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Roll Angle Vs. Flight Mode 
(6 Full Transition Flights) 
Sept 16 2014 
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GI-10 Full Transition!! 



GL-1 0 Testing 

► Prove our efficiency calculations and the 
transition corridor of the model at the 50% 
scale 

► Using COTS controllers we used a more rapid 
control research development 

► Prove the VTOL capabilities and long 
endurance can be combined 



GL-10 Completes FULL transition!!! 



Transition Corridor Testing 

► We expanded the transition corridor by 
reducing the tail volume coefficent to validate 
we still have 



Removal of the horizontal tail “add on’s” 
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Red line is where the alpha was at 60deg nose 
up due to less tail volume, and the trim 
condition was off in the controller. I followed 
proper protocol and saved the model. 




Check out the Youtube VIDM 

https: // www.youtube. com /watch?v=kXq!26sF5uc 
986k VIEWS!!!!! 
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Chapter 1 

Introduction 


Flight has always captured man’s imagination. This is evidenced by the great variety of aerial 
vehicles that exist today. Everything from fixed-wing to rotorcraft; satellites to spaceships; 
mono-wing to quadrotor. However, despite the wide variety of hying vehicles, not one of them 
has attained eternal flight. Accomplishing this feat is one of the great challenges still facing 
the aviation community. 


Motivation 

Achieving eternal flight opens the doors to atmo- 
spheric satellites. Existing satellites have a great 
number of capabilities that enrich our lives; however, 
their distance from the surface of the earth precludes 
certain types of transmission capabilities. Once eter- 
nal flight is achieved, that vehicle can serve the same 
role as ordinary satellites, but its close proximity 
will allow for real time two way communications, 
like wireless broadband internet. And with active 
controls, atmospheric satellites would not be con- 
strained to geosynchronous orbits, like our existing 
satellite technology. 

Objective 

Many projects are under way to achieve this goal; 
however, most of these research efforts follow the 
same design methodology, and have exhausted the 
limits of this particular design. This concept in- 
troduces a completely new aerial vehicle structure, 
which uses the best features of fixed-wing and rotor- 
craft designs. Combining the best features of differ- 
ent classes of aircraft, expands the capabilities be- 
yond what either one can achieve on its own. 


Outline 

This research introduces a revolutionary concept ve- 
hicle that shows great promise in achieving eternal 
flight capabilities. This chapter introduces the con- 
cept vehicle and outlines its features and benefits. 

Section 1.1 presents the Executive Summary 
which outlines the past and present state of satel- 
lite technology. Current efforts are under way to ex- 
pand these capabilities, but an alternative approach 
is required to achieve the objective. 

Section 1.2 describes the Goal of Eternal Flight 
Different types of satellites provide various functions, 
but a technology breakthrough is essential to satisfy 
the existing commercial interest. 

Section 1.3 outlines the Mission Objective which 
decomposes the problem. It illustrates the majority 
of existing concepts are based on Human Powered 
Aircraft designs, and describes several competing de- 
signs currently underway. 

Section 1.4 presents the Technical Challenge of 
eternal flight. It shows that improvements in energy 
capture and energy storage alone are not enough to 
attain the goal. 

Section 1.5 introduces a new Concept Approach 
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called the Centrifugally Stiffened Rotor (CSR). It 
provides an overview of the concept vehicle and de- 
scribes the individual components. 

Section 1.6 performs the CSR Feasibility Analy- 
sis. It outlines the benefits of centrifugal stiffening, 
and illustrates how the CSR achieves high Lift to 
Drag ratios with low parasitic drag. 

Section 1.7 conjures up Unique New Missions ap- 
plicable to this concept vehicle. Low disk loading, 
large swept areas, and a stationary hub, make this 
attractive for several different types of projects. 

1.1 Executive Summary 

Satellites have fundamentally changed our lives, pro- 
viding dramatic improvements to our quality of life 
that could have been accomplished in no other way. 
New research efforts are attempting to satiate de- 
mand for increased capabilities, but no existing plat- 
form has fully achieved eternal flight. This research 
presents a novel alternative approach, which exceeds 
existing capabilities, outpaces current research ef- 
forts, and may be the answer to attain eternal flight. 

1.1.1 Existing and Future Demand 

Over the past 50 years, satellites provided worldwide 
telecommunications to connect the world, geoposi- 
tional navigation to permit precise location determi- 
nation, and imaging and sensor systems that have 
let us develop a depth of understanding of our planet 
as well as the entire universe. While incredibly ex- 
pensive, and by far the most expensive aerospace 
products on a dollar per pound basis, their unique 
ability to provide a “God’s eye view of the world” 
also provides productivity that is difficult to match. 

Yet there is need for even greater capabilities, 
and at lower costs. Society has developed an in- 
creased appetite for omnipresent communications 
and surveillance capabilities, but with the added de- 
sire for these services to be achievable in real-time as 
on-demand services. These new capabilities are not 
well met by existing satellites, but instead require a 
new solution that can achieve several orders of mag- 
nitude lower signal latency, with greater operational 
flexibility, while at dramatically lower cost. 


1.1.2 Present Research Efforts 

The question this study asks is whether a new form 
of satellites is possible, and proposes a unique ad- 
vanced concept approach which has never been pre- 
viously explored. Low Earth Orbit (LEO) satellites 
are about 50x closer than Geosynchronous (GEO) 
satellites; Atmospheric Satellites (AT-Sats) are an- 
other 50x closer than LEO satellites. 

AT-Sats are aerial platforms that operate at al- 
titudes of around 60,000 feet where air density and 
winds aloft offer the operating conditions requiring 
the lowest power without weather concerns. Several 
concept approaches have been developed over the 
past decade through efforts such as DARPA Vulture 
and ISIS. Each of these DARPA programs exceeded 
$150 million, and had the intent of achieving a large- 
scale flight demonstrator. Vulture researched High 
Altitude Long Endurance (HALE) fixed wing con- 
cepts that utilized ultra-lightweight structure tech- 
niques pioneered in Human Powered Aircraft (HPA) . 
ISIS researched HALE airships that use Lighter- 
Than-Air (LTA) structures and buoyant gases. 

Both approaches applied advanced solar cells to 
provide a renewable energy source with the objec- 
tive of staying aloft for years at a time. Both of 
these programs were cancelled due to lack of feasibil- 
ity in closing the day-night energy balance required 
for staying aloft without excessive altitude loss. Be- 
cause of this, new concept approaches are desired 
that can achieve lower cruise power requirements to 
remain airborne throughout the year, even at high 
latitudes, while accommodating seasonal wind and 
solar flux variations. 

1.1.3 Alternative Approach 

A novel vehicle concept approach was developed 
with the objectives of achieving a 50% reduction 
in the structural mass while maintaining a robust 
structural configuration, as well as enabling a 33% 
increase in the Lift to Drag ratio. The approach 
uses the tensegrity principle to attempt to put as 
much of the structural loads into tensile load paths, 
where advanced carbon composites provide incred- 
ible strength to weight, while minimizing compres- 
sion load paths. The method to achieve this is to 
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rotate several wings around a stationary hub at very 
large diameters via a tether at quite low airspeeds, 
but with a concentration of the large battery mass 
required for night energy storage at the outer ex- 
tremities. The concept is called the Centrifugally 
Stiffened Rotor (CSR) because the resulting cen- 
trifugal forces applied to the tether and wing limit 
the bending through in-plane stiffening. 

Three rotor-wings rotating around the hub pro- 
vide maximum concentration of centrifugal force, 
with a minimum number of wings to form a stable 
planar summation of forces and moments. Because 
of this structural stiffening, wing aspect ratios (the 
ratio of span to chord length) can be dramatically 
increased while using thin airfoil sections to achieve 
extremely low levels of induced and parasite drag. 

Traditional rotors achieve poor aerodynamic effi- 
ciency because of poor, triangular span loading due 
to a linearly decreasing velocity distribution across 
the rotor that results in high induced drag. How- 
ever with this concept, the large radius of rotation of 
the rotor- wings with the tether comprising the inner 
80% enables the rotor- wings to achieve elliptic span 
loading. The tether takes advantage of advanced 
aerospace grade compressed Spectra providing in- 
credible strength to weight at small cross-sectional 
tether diameters that result in less than 10% of the 
drag being associated with the tether. 

Aerodynamic efficiency is further increased by 
taking advantage of a highly synergistic Outboard 
Horizontal Stabilizer (OHS) configuration layout to 
achieve thrust instead of drag from the tail pitch 
trim surface. Another favorable aerodynamic char- 
acteristic of the CSR concept is that since the outer 
rotor- wings are producing the lift and the inner hub 
is stationary, large payload volumes can be incor- 
porated into the center hub without incurring drag 
from the “cruise speed” required for lift generation. 

1.1.4 CSR Evaluation 

New tools were required to analyze this concept since 
this rotating multibody approach is far outside of ex- 
isting aeronautic systems analysis tool capabilities. 
Initial performance and aerodynamic analysis results 
of the CSR concept indicate compelling advantages 
for this approach in comparison to existing Vulture 


HPA flying wing solutions. LTA airships were not 
compared because prior studies have shown those 
approaches are unable to provide station keeping ca- 
pability for a large portion of the year due to winds 
aloft, and the high parasite drag associated with the 
large volumes required at high altitude for buoyancy. 

While the CSR concept performance results are 
compelling and indicate feasibility to achieve day- 
night energy balance, the key concept risk identi- 
fied early on is the ability to maintain control of the 
multibody dynamics that are inherent in the design. 
Therefore, the majority of project scope focused on 
performing detailed control modeling and simulation 
of the system with non-ideal, real world disturbances 
which introduce dynamic instabilities. The control 
efforts are discussed in detail, as this is a very unique 
multibody controls problem with methods developed 
that can be applied to other problems of interest, 
such as asteroid capture with tethered systems. 

In order to validate the developed control algo- 
rithms, a subscale prototype testing effort has been 
initiated. The CSR subscale experiment will estab- 
lish control feasibility, but was outside of the scope 
of this Phase I study and therefore is only outlined 
at the end of this report. The CSR controls effort 
has led to a PhD topic by the lead control researcher, 
which will be completed in the following months as 
part of his final dissertation. 

1.2 Goal of Eternal Flight 

A variety of existing satellites fulfill different technol- 
ogy roles. Corporations have shown a sincere desire 
to expand the realm of satellite capabilities, to meet 
their ever expanding needs. 

1.2.1 Technology Breakthroughs 

Google initiated a “Solve for X” conference in 2012 
(https://www.solveforx.com), to engage in radical 
moonshot visions where Huge Problems intersect 
Breakthrough Technologies. Such ideas are already 
being developed in the little known Google-X Lab, 
and have started to achieve huge impacts in technol- 
ogy and mission areas that are critical to aerospace. 

The Google “Driverless Car” is one of these 
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breakthrough attempts to turn-on- it ’s-head the pre- 
conceived notions of what is not only possible, but 
what we will soon wonder how we ever lived without. 
Google has already achieved what critics said would 
be a ‘hard stop’ to such idealistic visions by making 
robotic cars legal in 4 states; with almost every auto 
company now mirroring their autonomy goals. 

The “Solve for X” conference has many themes 
congruent to NASA’s mission; such as Cheaper Satel- 
lite Pictures , Low Power Broadband Wireless Every- 
where , Inexpensive Renewable Energy , and Reversing 
Climate Change. What if an aerospace mission con- 
cept had the potential to achieve all these moonshot 
goals? What if NASA could show the potential to 
dramatically impact many of the top national chal- 
lenges with a single mission concept path and tech- 
nology suite? How could this be possible...? 

1.2.2 Expanding Needs 

The first satellite was placed into orbit in 1957, and 
within 15 years satellites became highly capable plat- 
forms providing essential services through standard- 
ized satellite bus designs. Satellites have fundamen- 
tally changed our lives, providing dramatic improve- 
ments to our quality of life that could have been 
accomplished in no other way. Over the past 50 
years, satellites provided worldwide telecommunica- 
tions connecting the world, geopositional navigation 
permitting precise localization, and imaging and sen- 
sor systems that let us develop a depth of under- 
standing of our planet as well as the entire universe. 

While incredibly expensive, and by far the most 
expensive aerospace products on a dollar per pound 
basis, their unique ability to provide a “God’s eye 
view of the world” provides productivity that is dif- 
ficult to match. Yet there is need for even greater ca- 
pabilities, and at lower costs. Society has developed 
an increased appetite for omnipresent communica- 
tions and surveillance capabilities, with the added 
desire for these services to be achievable in real-time 
as on-demand services. These new capabilities are 
not well met by existing satellites, but instead re- 
quire a new solution that can achieve several orders 
of magnitude lower signal latency, with greater op- 
erational flexibility, while at dramatically lower cost. 


1.2.3 Geosynchronous Satellites 

Geosynchronous (GEO) satellites are able to remain 
over a fixed location above the Earth for continuous 
surveillance or communication across a large cover- 
age area because of their extreme distance. Because 
of their distance away from the Earth’s surface, a 
long transmission delay (or signal latency) occurs. 
This delay makes the use of GEO satellites for two 
way communications impractical; however, can pro- 
vide a large, constant coverage area with one way 
communications with great effectiveness. Due to the 
large transmission distances, the power required is 
considerably increased due to power loss being in- 
versely proportional to the transmission distance. 
With such a broad coverage area, frequency reuse 
is extremely limited. Coverage at high latitude lo- 
cations is challenging due to very low looking angles 
from the satellites, which increases signal obstruc- 
tion and losses from building, trees and mountains. 

1.2.4 Low Earth Orbit Satellites 

Low Earth Orbit (LEO) satellites are significantly 
closer proximity to Earth, with a decrease in cover- 
age area, but lower signal latency. However, LEO 
satellites are required to constantly orbit, and are 
unable to maintain a fixed reference location above 
the Earth. This necessitates large fleets of LEO 
satellites to maintain continuous coverage over any 
specific location, with very high associated capital 
costs for both the satellites and the cost to launch 
and insure the satellites. While the signal latency is 
reduced, it is still extremely challenging to use for 
two-way digital data transmission with an inability 
to achieve broadband speeds. 

1.2.5 Atmospheric Satellites 

Atmospheric satellites (AT-Sats) are a new breed 
of aircraft that operate at altitudes around 60,000 
feet, where air density and winds aloft offer operat- 
ing conditions that require the lowest power without 
weather concerns. AT-Sats operate over relatively 
fixed locations for continuous coverage with a sin- 
gle vehicle. Because of their proximity to Earth, the 
signal latency is extremely low, along with the trans- 
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mission power required (enabling the entire payload 
to shrink in size). The coverage area is greatly re- 
duced, essentially down to the size of a single urban 
area with the most direct lines of communication or 
the size of several states for the full coverage area. 
Since the aircraft is able to takeoff and land, the 
cost of deployment (and recovery) is very low, with 
the added ability to upgrade and refurbish satellite 


payloads as equipment ages. Likely the total system 
costs of AT-Sats are several orders of magnitude less 
than other satellite solutions, however current High 
Altitude Long Endurance (HALE) aircraft are only 
able to maintain flight for several weeks at a time, 
and lack the ability to provide the highly reliable and 
persistent capabilities that satellites offer. 



Orbital period 


15 hours 


5 hours 


Iridium 

Mm 10 Mm 


Hubble 


Radius of 


20000 


10000 miles 


'10000 mile: 


W-\ss 1 -- 

25000 km/h 


15000 mph 


-20000 km/h 


pOOO mph — - 


Geostationary 

20 hours \E a rth Orbit 

COMPASS 

MEO stellites 

GLONASS 


Orbital speed 


Figure 1.1: Existing GEO/LEO satellite orbit: height above Earth, period and speed. 
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Satellite 

GEO 

LEO 

AT-Sat 

Example 

XM 

Ecostar 

Motorola 

Iridium 

AV Global 
Observer 

Distance above 
Earth (mi) 

22,230 

500 

10 

Coverage 
Diameter (mi) 

7926 

2300 

500 

Fleet 

Required 

1 

50-60 

1 

Satellite 

Cost 

$300M 

$50M x50 

A 

to 

o 

£ 

Launch & 
Insurance Cost 

$100M 

$25M x50 

< $1M 

Ideal Signal 
Delay 

>250 

5-40 

0.1 


Table 1.1: GEO, LEO, and AT-Sat comparison. 


1.2.6 Commercial Interest 

Satellite solutions vary in their strengths and weak- 
nesses, with each type providing value. The or- 
bits of existing satellites is depicted in Figure 1.1. 
GEO satellites offer a continuous capability with 
a single satellite, LEO offers 50x closer proxim- 
ity for improved single latency, power requirements, 
and imaging resolution. AT-Sats offer another 50x 
closer proximity than LEO satellites for even bet- 
ter latency, power, and resolution; while also of- 
fering reusability and extremely low capital costs. 
The summary in Table 1.1 provides a comparison 
between the different varieties of satellites. 

Referring back to two of the highest priorities of 
Google’s “Solve for X” challenges, Cheaper Satellite 
Pictures and Low Power Broadband Wireless Every- 
where ; AT-Sats offer an incredible match to these 
desired goals. Just within the period of this Phase 
I study, both Google and Facebook have entered 
into AT-Sat development, with each purchasing 
aerospace companies focused on aircraft for this mis- 
sion. Google purchased Titan Aerospace in April, 
2014 (http://techcrunch.com, 2014/04/14, Google 
Acquires Titan Aerospace, the Drone Company Pur- 
sued by Facebook ), with Facebook buying Ascenta 
(a spin-off of the QinetiQ Zephyr project) in March, 


2014 (http:/ /www. theguardian.com, 2014/03/28, 
Facebook Buys UK Maker of Solar Powered Drones 
to Expand Internet ). 

Google also initiated test flights from their 
Google-X Lab Project Loon this year, which utilizes 
a large network of balloons that drift with the winds 
while using limited directional control by changing 
the buoyancy (and resulting altitude which have dif- 
ferent wind speeds and directivity). The Project 
Loon approach is severely limited as it requires un- 
controlled flight over countries, and is therefore only 
being tested in the southern hemisphere. The bal- 
loons also only have a flight duration of 50 to 100 
days before requiring descent and payload pickup, 
with the trajectories starting at the most southern 
latitudes and spiraling their way eventually to the 
equator with the prevailing winds. Clearly these in- 
vestments showcase the interest that exists in devel- 
oping alternative platform solutions to satellites. 

1.3 Mission Objective 

Several research groups, both commercial and gov- 
ernment sponsored, are attempting to design feasible 
AT-Sat solutions which achieve eternal flight. The 
problem reduces down to an energy balance relation- 
ship, where HPA techniques dominant the majority 
of current design paths. 

1.3.1 Problem Decomposition 

The basic problem of AT-Sats is decomposed into 
Power Required = Vehicle Drag x Flight Velocity . 
Thus, to minimize the power required to remain 
aloft, both the vehicle drag and flight velocity need 
to be minimized. 

Minimizing drag means flying at the highest 
cruise altitudes where the Lift to Drag (L/D) ra- 
tio is maximized. Altitudes where the air density 
is low, but not so low that there is insufficient den- 
sity to maintain flight. High structural efficiency is 
achieved with lightweight motors and system com- 
ponents, as well as with high efficiency components 
that experience the least energy losses. 

The key determining factor for vehicle feasibility 
is based on achieving Energy Collected in Daylight > 
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Energy Expended During Night Operation + Energy 
to Maintain Flight. Maintaining this energy balance 
is most difficult during the winter months at high 
latitudes, where there is low solar flux and low in- 
cidence angles on the solar cells. Wind velocities 
are minimized by operating at altitudes between 60 
and 70 thousand feet. However, even operating at 
these altitudes, when considering greater than 99% 
of days, the winds aloft can reach 120 miles per hour. 

1.3.2 Human Powered Aircraft Approach 

Conventional Human Powered Aircraft (HPA) con- 
cepts are able to achieve L/D ratios of about 35 at 
optimal aspect ratios of around 25. The aspect ratio 
is a dimensionless parameter that tracks the span to 
chord ratio of the wing and is defined as span squared 
divided by the wing area. L/D is an effective metric 
for tracking the overall aerodynamic efficiency that 
includes the parasite drag (skin friction drag associ- 
ated with wetted area) and induced drag (both 2D 
and 3D pressure and vorticity drag sources resulting 
from creating lift). 

The resulting payload mass fractions of the HPA 
inspired aircraft is less than 5%, with about 45% of 
the mass required for the structure and about an- 
other 45% allocated for the propulsion system, all 
while using the most advanced solar cell and bat- 
tery technologies. The growth factor (a dimension- 
less weight sensitivity parameter that compares the 
resulting increase in vehicle weight that results from 
increasing the payload weight by one pound) of such 
concepts are quite high, which is an indicator that 
the feasible design space is quite restricted. 

1.3.3 Best In Show 

The current world record holder for High Altitude 
Long Endurance (HALE) flight was set four years 
ago by the QinetiQ Zephyr 7 aircraft. Figure 1.2 
illustrates this vehicle, along with the other con- 
cept approaches currently pursued. The Zephyr is 
a HPA inspired structure that maintained high al- 
titude flight for a two week period in midsummer 
in New Mexico. The aircraft was 110 pounds, with 


a wing span of 74 feet, and had a wing loading of 
less than 1 pound per square foot. The payload was 
less than 5 pounds, and while the Zephyr achieved 
70,741 feet during flight, it was not able to main- 
tain altitude during night operation. Essentially, the 
aircraft used its altitude as potential energy during 
night operation, when energy was being used from 
the batteries which had been stored from the solar 
cells during day light operation. 

AT-Sat operational products cannot utilize such 
altitude energy storing. Lower altitudes limit com- 
munication line of sight and weather effects impede 
the ability to station keep. The Zephyr is highly re- 
flective of all the other aircraft concept approaches 
currently pursued and provides many indicators re- 
lating to future concept feasibility. 

The Zephyr was only able to maintain flight for 
two weeks because they used the most advanced 
lithium-sulfur batteries. The batteries were capable 
of about 530 Watt hours per kilogram, which is more 
than twice the specific energy of any other lithium 
rechargeable batteries at that time. Lithium-sulfur 
chemistry only permits up to approximately 20 
recharge cycles before seriously degrading their spe- 
cific energy. Also, they are only capable of providing 
very low power discharges rates below a 1 C rating, 
where C is the battery rating that relates to how 
quickly a battery can discharged in comparison to its 
charge time, and a 1 C rating means the discharge 
rate can only be as high as the charging rate. 

By far, the two largest mass fractions for the 
Zephyr were the batteries and wing spar. Since HPA 
structures have done everything possible to achieve 
ultra-lightweight structures for conventional flying 
wing configurations, the major focus for research 
programs such as DARPA Vulture has been to opti- 
mize the energy system. These efforts currently fo- 
cus on high specific energy regenerative fuel cells (en- 
ergy storage) and high efficiency/ lightweight solar 
cells (energy capture). The Zephyr project has just 
recently been purchased by Airbus (http: //defense- 
update. com, 2014/04/25, Airbus Invites Partners 
to Test Mission Payloads with Zephyr 8 Rseudo- 
Satellite’) and a new Zephyr 8 is under development 
as a follow-up next generation airframe. 
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Figure 1.2: Current HALE AT-Sat approaches under development. 


1.3.4 DARPA Concepts 

Several concept approaches have been developed 
over the past decade through efforts such as DARPA 
Vulture and ISIS. Each of these DARPA programs 
exceeded $150 million, and intended to achieve a 
large-scale flight demonstrator. 

Vulture researched HALE fixed wing concepts 
that utilized ultra-lightweight structure techniques 
pioneered in HPA, where concepts follow design and 
construction techniques pioneered by Paul Macready 
with the Gossamer Albatross 35 years ago. The air- 
craft have very low wing loading which provides sub- 
stantial surface areas for solar cells, and minimizes 
the cruise velocity and required power. The struc- 
tures are extremely delicate with Mylar surfaces for 
much of the wing surface, but they still have huge 
root bending moments absorbed through these thin 
wings. This approach results is a relatively small 
payload and has great sensitivity to winds. 

ISIS researched HALE airships that use Lighter- 
Than-Air (LTA) structures. LTA concepts utilize 


enormous hydrogen/ helium filled vehicles to achieve 
a difference in buoyancy between the internal vol- 
ume, and the already very low density air at high 
altitudes. Such large volumes require equally large 
surface areas. If winds aloft could be ignored, such 
approaches would be ideal. However during the ma- 
jority of the yearly seasons, and at the majority of 
latitudes, LTA approaches yield such high drag that 
the power requirements at even modest winds exceed 
a feasible solution. 

Both approaches applied advanced solar cells to 
provide a renewable energy source with the objec- 
tive of staying aloft for years at a time. Both of 
these programs were canceled due to lack of feasibil- 
ity in closing the day-night energy balance required 
for staying aloft without excessive altitude loss. Be- 
cause of this, new concept approaches are desired 
that can achieve lower ‘cruise’ power requirements to 
remain airborne throughout the year, even at high 
latitudes, while accommodating seasonal wind and 
solar flux variations. 
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1.3.5 Commercial Concepts 

Clearly a great deal of momentum and interest is 
building in AT-Sat vehicle development, with a focus 
by all (except the Google Loon Project) on conven- 
tional HPA conventional wing configurations. Sev- 
eral variants of wing-tail configurations are shown in 
Figure ??, but each concept maintains a highly sim- 
ilar HPA structural approach while span-loading the 
vehicle weight across the span as much as possible 
to decrease the wing bending moments. 

Lockheed utilizes a bottom strut across the fuse- 
lage element to achieve greater effective wing spar 
depth; however, without trussing between the fuse- 
lage elements this approach does not provide sig- 
nificant improvement. Furthermore, wire crisscross- 


ing trusses between each fuselage element at flight 
speeds would dramatically increase drag. 

Aurora utilized a combinatory Z-wing arrange- 
ment. Each wing element would dock together in 
flight with a hinge element that does not transmit 
bending moments. Also, each wing element would 
fly at inclination angles better suited to capture the 
solar flux for energy storage. 

Boeing utilized an Outboard Horizontal Surface 
(OHS) to minimize the induced drag and augment 
the pitch trimming capability. Only Boeing was se- 
lected for the Phase II DARPA Vulture program. 
However, after detailed studies, this concept showed 
questionable feasibility and was descoped into a fo- 
cused effort to advance the energy storage and cap- 
ture technologies. 
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Figure 1.3: Pareto sensitivity frontier of energy required compared to energy available. 
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1.4 Technical Challenge 

A number of improved technologies are required to 
achieve feasible AT-Sat operational capability. Key 
advances are required to increase reliability, decrease 
payload weight, and expand current energy collec- 
tion and capture efficiency. 

1.4.1 System Reliability 

One such topic area includes achieving satellite-like 
reliability. Enabling continuous flight for a year or 
more requires confidence across all aircraft compo- 
nents. This reliability must be accomplished in all 
systems including: energy storage, control surface 
actuators, flight control avionics, and sensors that 
provide an autonomous flight capability. 

1.4.2 Light Weight Antennas 

Another technology area includes ultra-lightweight 
antenna systems. While targeting eternal flight, 
these improved systems can minimize the payload 
weight. Phased- array multi-beam systems of 10 or 
even 100 beams will likely be required for cell phone 
providers to support 100,000 concurrent lines for mil- 
lions of customers within the coverage area. 

1.4.3 Energy Capture and Storage 

Performing technology sensitivity analysis of the 
conventional wing configurations helps to under- 
stand why DARPA descoped the Vulture program. 
Examining all these sensitivities is beyond the scope 
of this Phase I study. However, NASA/TP-2007- 
214861 (High Altitude Long Endurance UAV Anal- 
ysis of Alternatives and Technology Requirements 
Development) is an excellent source that provides 
comprehensive understanding. 

A key result from this study is a Pareto tech- 
nology frontier chart that shows the sensitivities of 
the two primary energy parameters: collection and 
storage. Figure 1.3 presents the Pareto sensitivity 
frontier of the energy required compared to energy 
available for completing a day/night cycle during the 
worst winter day conditions. They are specified as a 
function of the specific energy of the energy storage 


system and the solar cell efficiency. The infeasibility 
of this mission concept at current technology levels 
is clearly seen with only about 34% energy available 
versus energy required for accomplishing a day /night 
cycle on the worst winter day. 

Currently, lightweight thin-film solar cells are at 
32% efficiency (Micro Link Devices) and the best 
lithium rechargeable batteries that are capable of at 
least 365 cycles (one year of operation for an AT- 
Sat) are nearly 280 Whr/kg (Panasonic/NCR cells 
with the 18650 form factor). 

If energy storage and solar cell technology are the 
only improvements provided for future AT-Sat vehi- 
cles, the required feasibility to achieve a 100% energy 
balance requires battery specific energies greater 
than 1000 Whr/kg and solar cells that exceed 75% 
efficiency. Therefore, it is highly desirable to not 
only focus technology development efforts on im- 
proved energy, but also on a new concept approach 
that could dramatically lessen the energy required 
for sustained high altitude flight. 

1.5 Concept Approach 

The proposed concept has never been attempted be- 
fore. This vehicle is a new approach that dramati- 
cally lowers structural weight, greatly reduces drag, 
and provides a far more robust structure in weather 
than existing HPA structures. The approach uses 
three extending and retracting tethered wings which 
rotate around a stationary central hub, and the ro- 
tation balances all the forces and moments within 
the system. Figure 1.4 depicts the vehicle, and a 
detailed description is provided below for each of 
the three major components: the rot or- wing satel- 
lite components, the central stationary hub, and the 
tethers that connect them together. 

1.5.1 Rotor-Wing Satellites 

All the controls, and the majority of the vehicle mass 
is located within the rotor-wing satellites. This sec- 
tion provides the general description of these struc- 
tures, why three arms are used, and how they impart 
controls on the entire CSR system. 
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Figure 1.4: Centrifugally Stiffened Rotor (CSR) concept, with zoomed in view of the rotor- wing. 


General Description 

Three rotor-wing satellites are coupled together 
through a high tensile, small diameter cable tether. 
This number provides the least number of rotor- wing 
satellites necessary to achieve a stable and control- 
lable tip path plane, while achieving the maximum 
tether centrifugal force. The rotor- wing satellites are 
located at the ends of each tether, and provide all 
the lift, propulsion and control for the CSR system. 

Approximately 80% of the entire vehicle system 
mass is located in these components, which is domi- 
nated by the battery mass required to provide energy 
during night time operation without altitude loss. 
The battery mass is located at the outer fuselage to 
maximize the centrifugal force and provide stiffening 
of the rotor- wing structure and alleviate the bending 
moment that is exerted on the rotor- wing spar. 

Each wing-body is driven by a propeller with op- 
timum twist and pitch for the rotation rate. Since 
the rotor-wing thrust force is a tip driven system, 
there is no torque transmitted back to the central 
hub. This propulsion technique is opposite that of a 
conventional helicopter rotor, where a torque is ap- 
plied to the central shaft. This requires the need for 
a counter-torque on the traditional helicopter body 
to prevent the helicopter body from rotating in the 


opposite direction. There have been many exam- 
ples of tip-driven helicopters in the past that have 
avoided the need for a tail rotor. 

Why Three Arms? 

More than three satellites could be used, but it 
would result in a poorer performing system. A sys- 
tem of three components maximizes the rotor-wing 
Reynolds number along the chord, which is an im- 
portant factor for achieving the lowest drag with 
HALE aircraft. Additional rotor- wings would also 
spread the battery mass to more components. This 
results in a lower centrifugal force, which is the key 
method for alleviating the bending moment seen by 
the rotor- wing spar. 

Centrifugal stiffening also allows for relatively 
thin wing airfoils to achieve increased laminar flow. 
These thinner airfoils are only possible due to the 
inherent wing stiffness associated with the centrifu- 
gal stiffening. Additionally, a lower centrifugal force 
(mass x velocity 2 / radius) on the rotor-wing-body 
component results in a higher bending moment. A 
higher centrifugal force could be achieved with a 
higher rotational velocity of the rotor-wing, but this 
is not desirable as the required power would increase. 
Therefore, the design objective positions the major- 
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ity of the vehicle mass as far outboard along the 
radius of rotation as possible, while maintaining a 
low rotor-wing speed. 

Having more than three rotor- wing satellites also 
increases the difficulty of maintaining control of each 
element, because each component can independently 
vary the lift, thrust, and radial force. More elements 
may provide increased control robustness (through 
having the rotor tip path plane over constrained), 
but at the cost of both the vehicle performance and 
hub complexity. 

When wind conditions are present, each satel- 
lite will operate at different speeds. One element 
is essentially acting in a retreating blade condition, 
while another is acting in an advancing blade condi- 
tion. Therefore, the thrust of each element in wind 
conditions varies throughout each cycle, and the az- 
imuthal locations of each satellite element is non- 
equally spaced. Proving that the entire multibody 
vehicle system can maintain control is a key uncer- 
tainty of this concept that needs to be answered. 

Controls Overview 

Each rotor- wing is covered with thin film solar cells 
to collect solar energy during daylight hours. This 
provides power to the propulsion system for thrust, 
which permits the rotor-wings to rotate, and to keep 
the vehicle aloft. Each satellite includes an elec- 
tric motor and propeller in the nose of the body for 
propulsion. Control surfaces are oriented vertically 
and horizontally to maintain vehicle control. 

Pitch control of each wing changes the angle of 
attack of each individual wing. This is accomplished 
by a surface extending outward from the wing, and in 
the same plane as the wing. This pitch trim surface 
is the equivalent of a tail on a conventional aircraft; 
however, it operates in the wing wake vortex up- 
wash, so the required angle of attack decreases and 
the drag is further reduced by providing increased 
spanwise efficiency. When all wing pitch trim sur- 
faces change together, they provide the equivalent of 
helicopter collective pitch which increase or decrease 
the altitude of the vehicle. 

Roll control of the entire vehicle system is ac- 
complished by changing the pitch of each rotor-wing 
differentially. This achieves the equivalence of cyclic 


pitch change of a helicopter rotor, and enables the 
tip path plane to either tilt side-to-side or back-to- 
front. Yaw control is unnecessary for this vehicle 
since there it is axisymmetric about the center of ro- 
tation, and because the rotation is induced at by the 
wingtip propulsors. When winds are present, the en- 
tire vehicle system is angled through differential lift 
to provide an angular incidence between the forward 
wind direction and the effective tip path plane of the 
rotor- wings to provide forward thrust to keep the ve- 
hicle in place. 

A vertical surface is located on the body at the 
end of each rotor-wing in a similar fashion to a 
winglet, but extending both upwards and downwards 
from the wing. This surface provides the ability 
to augment or decrement the centrifugal force to 
provide stability across the different rotor-wing el- 
ements through a force that is independent of the 
rotation rate or thrust. 

The method of control for the entire multi-body 
vehicle system is envisioned where the azimuth lo- 
cation of each element is known to a controller 
through wireless signals communicated between each 
element. However, this may not be necessary and 
this is discussed in detail in the control development 
section of this report. 

1.5.2 Central Hub Payload 

While the center hub is described as stationary, it 
is permitted to rotate at the extremely slow angu- 
lar rates of the vehicle. Even for small Unmanned 
Aerial System (sUAS) scales of 55 pounds the rota- 
tional rate is less than 20 rpm and for larger satellite 
systems the rates are less than 5 rpm. 

The hub is located at the center of rotation and 
provides a payload pod independent of the wing 
cruise speed (i.e. large volume antenna pods that 
only result in drag during wind conditions). Since 
the tethers are operated with a large radius of ac- 
tion spanning hundreds of feet, the rotor rotation 
rate is quite low, with the payload rotating at that 
same low rate. 

The mass of the center pod is around 20% of the 
total vehicle mass, so the tether angle is relatively 
small (<15°) due to the centrifugal force balancing. 
Solar cells on the center hub provide the payload 
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energy, along with the potential to transmit electric- 
ity through the tether to the props depending on 
whether a more complex tether is warranted in case 
the optimal wing-rotor sizing does not permit suffi- 
cient solar cell area for propulsion purposes. How- 
ever, a conductive tether is significantly larger diam- 
eter and results in higher tether drag and weight. 

The hub has structural arms where each tether 
extends from the center section to provide rotation 
rate stability across each rotor-wing-body compo- 
nent. Specifically if one rotor- wing gets ahead of 
the others, a restoring angular force will be induced 
that slows that wing, while an angular force is in- 
duced on the other rotor-wings that accelerate the 
other rotor-wings that are lagging behind. These 
structural arms hold the tether, and act as ’whips’ 
to induce these stabilizing rotational forces across 
the entire wing/rotor system. While independent 
thrust at each rotor- wing permits differential spacing 
of the elements, the tether structural whips provide 
a mechanism for incurring a stabilizing force across 
the system of elements. 

1.5.3 Tethers 

The tethers connect the rotor-wings to the center 
hub. Because the vehicle is always in rotation, it only 
experience tensile loading. This is ideally suited for 
advanced composite and braided materials designed 
specifically for high tensile strength, like Spectra. 
The tether can be reeled in to achieve a compact 
configuration for takeoff and landing. But this con- 
figuration requires significantly more power due to 
the higher velocities needed to maintain lift and from 
the non-optimal triangular lift loading distribution 
across each wing. 

In flight, the tether is reeled out to provide a 
dramatically larger radius, which results in a toroid 
lift area and an optimal elliptical lift loading distri- 
bution across each wing. The amount of area the 
rotor-wing traverses is the effective disc loading of 
the vehicle and determines the power required. The 
tether comprises the inner 70 to 80% of the total di- 
ameter when in flight, depending of the optimized 
wing sizing for different mission needs. The tether is 
at a slower velocity than the rotor-wing, where the 
average tether velocity is only 35-40% of the rotor- 


wing speed. This means that only about 5 to 10% of 
the total drag of the vehicle comes from the tether. 
For example, a 0.03” diameter compressed aerospace 
grade Spectra has a tensile strength of 200 pounds. 

1.6 CSR Feasibility Analysis 

Having described the general concept of the CSR ve- 
hicle, the next task describes the feasibility of this 
system achieving eternal flight. This section de- 
scribes the benefits of centrifugal stiffening, and how 
the CSR achieves high L/D ratios while minimizing 
parasitic drag. It compares this concept with the 
leading competitors, and illustrates the distinct ad- 
vantages this concept offers. 

1.6.1 Benefits of Centrifugal Stiffening 

The first goal of the CSR concept is to achieve a 
large reduction in the wing spar mass fraction of the 
total vehicle mass, which increases the structural ef- 
ficiency. Another objective is to increase the Lift 
to Drag ratio, which increases the aerodynamic effi- 
ciency of the HALE aircraft. These goals are primar- 
ily required for station keeping missions. The CSR 
system achieves these goals through centrifugal stiff- 
ening, just as a helicopter rotor uses CF which allows 
extremely high aspect ratio blades while using very 
lightweight rotors. 

Most rotorcraft blades have aspect ratios of 30 
to 50, but the majority aircraft rarely exceed an as- 
pect ratio of 10, due to the high resulting weight of 
the wing and aeroelastic constraints. Rotor blades, 
which develop all the lift of a helicopter, typically 
weigh only 2 to 3% of the gross helicopter weight. 
Compare that to aircraft wings which are typically 
10 to 20% of the gross aircraft weight. HALE air- 
craft experience a similar wing weight mass fraction 
of 40 to 50% of the total aircraft mass. Clearly cen- 
trifugal stiffening of high aspect ratio rotors-wings 
is a highly effective method of taking lifting surfaces 
to significantly greater aspect ratios, while achieving 
very low lifting surface weight. 

While rotor blades use high tip speeds (600 to 
850 ft/sec) to achieve high radial accelerations, cen- 
trifugal forces can also be achieved at relatively low 
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CSR rotor-wing is free to rotate at tether joint instead of imparting large bending moments, 
with centrifugal restoring forces to achieve an efficient tensegrity structure 
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Figure 1.5: Lift load distribution on CSR concept versus conventional wing configuration. 


tip speeds by having a large mass at the rotor outer 
most tip location. Slowed rotor concepts such as 
the CarterCopter use depleted uranium at the rotor 
tips to maintain rotor blade stiffness at low rota- 
tional speeds. The centrifugal force in rotor blades 
results in the wing spar material being very small 
to counteract the rotor lift root bending moment. 
Helicopter rotors experience 1000’s of g’s of acceler- 
ation, while the CSR rotor-wing needs only about 2 
g’s in combination with the large tip mass to provide 
sufficient stiffening for minimal bending. 

The CSR provides a tensegrity structure ap- 
proach to achieve high aspect ratios instead of HPA 
structures which depend on large wing chords to pro- 
vide a thick structural spar depth. Spreadsheet and 
Matlab analysis methods were developed to provide 
a summation of forces and determination of tensile 
forces across the wing and tether. This method was 
first applied to existing HALE conventional wing 
structures to validate the ability of the method to 
capture accurate mass trends, and then applied to 
the CSR structure. 

Figure 1.5 shows the lift load distribution on the 
CSR concept versus a conventional wing configura- 
tion showing the distribution of forces and the abil- 
ity of the rotor- wing elements to pivot at the tether 
connection to alleviate a root bending moment. The 


NASA/AeroVironment Helios is shown as an exam- 
ple of HALE aircraft experiencing very high bending 
moments over large wing spans. The Helios failed 
during a test flight because the wing fluttered due 
to a lack of structural stiffness. 

1.6.2 Achieving High Lift To Drag Ratios 

An objective of the CSR design uses the high aspect 
ratio rotor- wing to achieve high L/D ratios. Heli- 
copter rotors achieve poor lift efficiency because their 
lift is developed over the entire diameter. This re- 
sults in a triangular load distribution with almost 
the entire lift generated by the outer 25% of the 
rotor blade. Achieving low induced drag depends 
on achieving an elliptic load distribution across the 
lifting surface. Helicopter rotors have high induced 
drag, because they do not achieve an elliptic loading. 

Helicopter rotors also experience high induced 
drag because each rotor acts in close proximity to 
the downwash of the blade rotating in front of it, 
which increases the angle of attack required for lift 
generation. Instead, the CSR concept has each satel- 
lite operating in a downwash field independent of the 
other rotor- wings due to the large spatial separation. 
An effective way of visualizing the resulting reduc- 
tion in downwash induced drag is to consider things 
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from a rotor disc loading perspective, as seen in Fig- 
ure 1.6. The CSR concept achieves a downwash field 
in the shape of a toroid of the traversed rotational 
rotor-wings based on the extremely large tether di- 
ameters. A conventional rotor with 20 foot radius 
achieves a swept disc area of 1256 square feet, while 
the same 20 foot radius with an 80 foot tether tra- 
verses a swept area of 11,304 square feet. This is 
far greater than that of the span of the rotor-wing 
alone. The total rotor disc area relates directly to 
the induced drag. 

In addition, a conventional helicopter rotor has 
a triangular lift distribution with poor aerodynamic 
efficiency. However, the tethered rotor- wing achieves 
an ideal elliptical load distribution because the en- 
tire satellite has a relatively constant velocity field. 
There is approximately 20% variation in airflow ve- 
locity, which can be made elliptic through a modest 
twist and taper distribution across the rotor-wing. 
The final method for minimizing induced drag uses 
the same Outboard Horizontal Stabilizer (OHS) lay- 
out used by the Boeing Vulture concept, which was 
originally developed by the WW II Blohm and Voss 
aircraft concept, and with application as recently as 
by Scaled Composites Spaceship One. 



Figure 1.6: Simplified traversed rotor disc area. 


1.6.3 Minimizing Parasitic Drag 

While the resulting low induced drag of the CSR 
rotor-wing is critical, operation at the maximum 
L/D requires that equal portions of induced and par- 
asite drag are achieved in balance. Therefore, it is 
also important to provide a mechanism for lower 
parasite (skin friction) drag. This is accomplished 
in two ways with the CSR concept, which are illus- 
trated in Figure 1.7. 

Conventional HALE concepts require that pay- 
load volumes move at the same speed as the rest of 
the vehicle. For bulky antenna or optic payloads this 
can result in substantial parasite drag. The payload 
for the CSR concept remains stationary at the hub, 
which is independent of the required lift generation. 
The size of the satellite is minimal due to the high 
density of the enclosed components (electric motor 
and batteries). 

Another mechanism of achieving lower parasite 
drag comes from the rotor-wings having low thick- 
ness to chord (t/c) ratios, because of the low result- 
ing bending moments of the centrifugally stiffened 
structure. The use of low t/c airfoils is important 
for HALE aircraft operating at low Reynolds num- 
ber. This condition maintains laminar flow for much 
of the airfoil while insuring that laminar separation 
does not occur. 

At first glance the tether drag may appear to 
significantly increase the parasite drag. The non- 
dimensional cross-sectional Coefficient of Drag (CD) 
of a cylindrical cable is extremely high (around 1.15 
compared to a typical airfoil of .01). As the tether di- 
ameter increases with larger scale vehicles, the CD of 
the tether improves through braided cross-sections. 
These have improved separation characteristics that 
achieve CD’s of about 0.75. 

The tether cable diameter is able to be small be- 
cause it is extremely structurally efficient at carrying 
tensile loads. Combined with the fact that the av- 
erage tether speed is only 40% of the rotor-wing tip 
speed, the resulting tether drag is only 5 to 10% of 
the total vehicle drag. Of course, these values de- 
pend on the amount of centrifugal tensile load that 
sizes the tether diameter and whether the tether is 
designed to conduct electricity for energy sharing 
across the payload hub and the rotor-wing elements. 
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Half Span 



Wing tip vortex provides upwash to pitch trim tail, resulting in thrust 

Why is the CSR concept more aerodvnamicallv efficient? 
• High span with ideal loading, without weight penalty 



if* 


Original Blohm & Voss 
Concept 


Scaled Composites 
Spaceship One 


OHS configuration layout, without weight penalty 

Use thinner airfoils to maximize low Reynolds 
number laminar flow, without weigh penalty 
Payload volume independent of cruise speed 
Tether drag is <10% of total drag 


Figure 1.7: CSR aerodynamic loading versus conventional helicopter rotor. 


A factor of safety of 6 was applied to the tether 
tensile load to insure reasonable tether cable diam- 
eter. Spectra offers a good match to CSR require- 
ments due to low stretching of about 1%. Spread- 
sheet and Matlab methods were assembled for cal- 
culating the parasite and induced drag of the CSR 
concept using a blade element aerodynamic analysis. 

1.6.4 CSR Comparisons 

The prior mass, force balance, and drag methods 
were combined along with methods to determine the 
power required for maintaining flight. The results of 
these analyses are presented within Table 1.2 across 
several payload sizes ranging from sUAS to full satel- 
lite equivalent payloads. The vehicle characteristics 
shown are not optimized, but merely a result of sim- 
ple sizing sensitivities to yield reasonable configura- 
tions. A useful tracking metric to compare the CSR 
concept results to other approaches is the power re- 
quired per vehicle pound to maintain flight. 

An interesting comparison reference to the CSR 


concept is the recent successful winner of the Siko- 
rsky Human Powered Helicopter prize (the Aerovelo 
Atlas), which lifted a total mass of 281 lbf (121 lbf ve- 
hicle plus 160 lbf pilot) for 64 seconds with a human 
power output of almost 750 watts (2.67 watts/lbf). 
This was accomplished with a total rotor lift area of 
13,700 ft, that had 4 rotors comprising a total of 156 
ft diameter, and an incredibly lightweight and flimsy 
structure that provided no design to ultimate load 
or safety margins. 

A total safety margin over the design load of 
3 was utilized for sizing the CSR structure, with 
the resulting specific power required being somewhat 
higher across all vehicle sizes than that required by 
the Atlas. At a similar vehicle weight the CSR con- 
cept is able to achieve around 3.5 watts/lbf. As an- 
other comparison, the Zephyr 7 is calculated to use 
approximately 550 watts at cruise for the 110 lbf 
vehicle, yielding a power required per pound of 5 
watts/lbf, which is 45% more power than required 
for a CSR concept at a similar mass. 
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sUAS 

Camera 

Small AT- Sat 

AT-Sat 

Geometry 

Diameter - Outer (ft) 

70 

200 

400 

800 

Diameter - Inner (ft) 

56 

160 

320 

640 

Span (ft) 

7 

20 

40 

80 

Chord (ft) 

0.5 

1.0 

1.0 

1.5 

Aspect Ratio 

18.7 

26.7 

53.3 

71.1 

Swept Lifting Surface (ft "2) 

1385 

11304 

45216 

180864 

Forces 

Rotational Rate (rpm) 

20 

7 

4 

2.25 

Rotor- wing Tip Speed - Outer (ft/sec) 

73 

73 

84 

94 

Rotor- wing Tip Speed - Inner (ft/sec) 

59 

59 

67 

75 

Reynolds Number (avg) 

233023 

466045 

532623 

898802 

Centrifugal Acceleration (g’s) 

4.52 

1.58 

1.03 

0.65 

Rotor-wing Angle vs Horizon (deg) 

-15 

-12 

-12 

-18 

Tether Tension (lbf) 

62 

122 

256 

823 

Weights 

Gross Weight (lbf) 

39.2 

226.8 

592.5 

2187.7 

Rotor- wing Weight /Element (lbf) 

13.3 

60.5 

158.0 

583.3 

Payload Weight (lbf) 

3.9 

22.7 

59.2 

218.7 

Swept Disk Loading (lbf/ ft ^2) 

0.028 

0.020 

0.013 

0.012 

Drag and Power 

Induced Drag (lbf) 

133 

646 

1363 

4835 

Parasite Drag (lbf) 

18 

101 

303 

1294 

Tether Drag (lbf) 

9 

35 

180 

602 

Power Required (watts) 

159 

783 

1847 

6734 

Specific Power (watts/lbf) 

4.07 

3.45 

3.12 

3.08 


Table 1.2: Comparison of CSR concepts across several payload weights. 


1.6.5 Analysis Results 

The results of this analysis indicate a 35% improve- 
ment in the L/D ratio, while at the same time achiev- 
ing a reduction in structural mass by 50% may be 
possible. Further gains may be feasible after opti- 
mization is applied to the current non-optimum re- 
sults. Reduction in power by 35% appears possible, 
which will dramatically increase HALE mission fea- 
sibility. The improvement in drag and weight are 
the critical concept performance ingredients for this 
mission concept to achieve feasibility. 

Applying these improvements potentially en- 
abled by the CSR concept provide a basis for achiev- 
ing a day/night energy balance while using near- 


term rechargeable batteries with specific energy of 
400 Whr/kg, and thin film solar cells with thin- 
film efficiency of 35%. Use the same tool method- 
ology developed in NASA/TP-2007-214861, but re- 
evaluate that analysis at the aerodynamic and struc- 
tural mass ratios indicated by the CSR concept. 
This result is visualized in Figure 1.8 where the CSR 
concept aerodynamic and structural factors are spec- 
ified for the same Pareto sensitivity frontier analysis 
of the energy required compared to energy available 
for completing a day/night cycle during the worst 
winter day conditions, as a function of the specific 
energy of the energy storage system and the solar 
cell efficiency (analysis performed by Mark Guynn 
of NASA Langley Research Center). 
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Figure 1.8: CSR concept specified for same Pareto sensitivity frontier. 


1.6.6 Remaining Questions 

The ability to achieve feasibility with near-term en- 
ergy solutions provides a compelling reason to fur- 
ther investigate the assumptions buried within this 
initial analysis. Two assumptions not evaluated 
in this vehicle performance analysis are the ability 
to achieve a 10% mass fraction for the center hub 
structure (including the tether reel-in/out mechani- 
cal system) , and the ability to maintain control with 
complex multibody dynamics present in a system 
that will experience nonuniform disturbances in real 
world operations. 

The remainder of the Phase I investigation was 
directed towards answering the control feasibility, 
with detailed development of a multibody analysis 
method. While not within the scope of Phase I, the 
researchers were also highly motivated to test the 


resulting control algorithms while also verifying the 
assumptions present in the hub mechanics. A sub- 
scale test vehicle has been initiated for this purpose; 
however, it is still in development and will be com- 
pleted prior to a Phase II proposal submission. 

1.7 Unique New Missions 

While the primary mission for the CSR concept 
would be for HALE communication platforms, a 
number of compelling missions have been identified 
that map directly into the unique capabilities of the 
CSR advanced concept. These unique advantages in- 
clude achieving an extremely high L/D ratio, a sta- 
tionary center hub, a large swept rotor- wing area, 
extremely low effective disc loading, and the ability 
to change the effective disc loading by reeling the 
tether in/out. 
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Dual-Use Airborne Wind Energy Platform and Coastal Ocean Observation 

Ground tethered wind energy extraction (conductive CNT tether key enabler) 
Rotor-wing tip propellers act as turbines to extract energy 



Cross-wind flight velocity scales with vehicle L/D 
What if... 

Offshore wind didn't need huge platforms; but just a buoy? 


7*2:18 PM ■ -0 
3,8/2012 


Conventional HAWT 
80 m 


r 600 


200 - 


7-27:57 PM 
3/5/2012 


72757 AM 

3/6/2012 


727:57 PM 
3/6/2012 


727:57 AM 
3/7/2012 


7:2757 PM 
3/7/2012 


727:57 AM 
3/8,2012 


32 - 

94400 AM 
3/5/201 2 


Makani/Google 

Airborne Wind Energy 
200 - 600 m 


What if... 

We could tap into faster wind speeds at higher altitudes 


Figure 1.9: Airborne Wind Energy (AWE) and coastal ocean observatory platforms. 


1.7.1 Airborne Wind Energy 

A number of companies are currently investigating 
Airborne Wind Energy (AWE) production through 
aerial platforms. These systems harness the wind at 
altitudes where the winds are faster, the turbulence 
is lower, and the winds are more consistently avail- 
able. Such systems require a tethered connection to 
the ground to harness wind energy, otherwise they 
would simply be blown downwind. 

Companies such as Makani (purchased by Google 
this past year) are developing concepts as shown in 
Figure 1.9. Small aircraft are flown in circles, acting 
as the tip of a ground-based horizontal wind turbine, 
since the tips are performing >75% of the energy 
capture. The FA A ground obstruction rules enable 
such tethered platforms to be operated at altitudes 
as high as 2000 feet. With the Makani approach, the 
entire tether length is moving at the aircraft veloc- 
ity, which is typically 200 to 300 ft/sec, resulting in 
high tether drag. Energy is extracted by propellers 
on the aircraft that act as turbines turning electric 


alternators, which then transmit the energy down 
the tether to the ground. 

Because the tether connects the aircraft to the 
ground, the tether orientation to the aircraft is 90 
degrees to each other so that the tether imparts a 
significant bending load on the wing. If the CSR con- 
cept were used in place of the conventional Makani 
aircraft concept, the ground tether would connect to 
the stationary CSR center hub and not experience 
a velocity or drag. The resulting drag of the sys- 
tem would be much less, and permit smaller land 
use as decreased tether angles to the horizon are ac- 
complished. In addition, the CSR approach would 
maintain the ability to keep the rotating tether be- 
tween the hub and rotor-wings along the span path 
to achieve a lighter weight solution. 

Makani has already accomplished testing at the 
10 to 50 kW scale in Hawaii, which has very high 
costs of electricity due to this coming from oil based 
energy production. In time, it may be proven that 
lower cost power can be achieved, but currently with 
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these smaller early adopter systems, the cost of en- 
ergy is quite high, especially when land use costs are 
considered. If these AWE platforms were able to be 
positioned in locations where the platforms could act 
as multi-functional surveillance/communications as 
well as AWE energy harvesters, there would be the 
ability to lower the effective costs. 

If a CSR platform was operated at 2000 feet, 
the communication transmission coverage diameter 
would be nearly 50 miles. The CSR platform could 
act as a cell phone tower or WiFi provider with a 
communication payload in the center hub. The visi- 
bility of the CSR platform would be far less than the 
alternatives of a tower or ground-based wind turbine. 

The CSR concept achieves a higher Lift/Drag ra- 
tio compared to other airborne wind concepts, which 
results in higher tether angles and less land area 
under the vehicle radius of operation that depends 
on incoming wind direction. Off-shore application 
is particularly appealing since a CSR aerial vehicle 
would eliminate large and expensive mooring plat- 
form, which is required for ground-based wind tur- 


bines. Ground-based wind turbines experience huge 
bending moments in the tower structure and require 
massive concrete pads, while an AWE transmits no 
bending moment down the tether. 

Off-shore also provides significant dual-use mar- 
ket opportunities. The first long-term civil use of a 
UAV was recently approved in Australia to fly along 
the coastline as a shark spotting platform. This 
was performed previously by small aircraft operat- 
ing with two people on board, the potential cost sav- 
ings are considerable. A series of off-shore CSR plat- 
forms could be performing a combination of energy 
harvesting, distressed swimmer surveillance and Wifi 
from each platform with very lightweight payloads. 

While such applications may be considered niche 
markets, hundreds of thousands of lifeguards in the 
world monitor coastal waters. Potentially there are 
even larger markets for this mission concept as it is 
applied to distributed residential power production, 
with a sUAS sized version likely capable of provid- 
ing 2 to 5 kW of power at altitudes of 500 feet with 
average wind speeds of less than 20 mph. 


Distributed Aperture Atmospheric Satellite Observatory 

High resolution linked sensors rotate around azimuth, at 450 and 600 ft radius locations 


Low resolution sensor in center payload section provides filling 



Result: Imaging resolution never before possible 


CSR Configuration 


What if... 

We combined new imaging techs with distributed rotating arrays? 


C* 

Azimuth 

Rotational 

Axis 


Low-Mass Planar Photonic Imaging Sensor 


Very Large Array (VLA) Radio Telescopes 


Figure 1.10: Distributed aperture atmospheric satellite observatory mission concept. 
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What if... 

Instead of attempting a 
Mars flyer only works for 
an hour, at transonic 
speeds 



Instead... 

We could fly along Mars 
canyon walls at slow 
speeds or even hover to 
look at all the layers of 
Martian geological strata 


Ultra Low Power. Compact VTOL Mars Explorer for Low Density Atmospheres 

Fully retracting version of the CSR rotor-wing to fit within a transit aeroshell 
Operation from a recharging station (with potential for CSR concept to capture wind energy) 



NASA Langley Rapid Prototyped 
Centrifugally Stiffened Rotor 


Initial 
prototype 
fully retracting 
CSR rotor 
(tested) 



Rotor-wing Retraction Mechanism 
CAD Design 


Figure 1.11: Ultra low power, compact VTOL Mars Explorer mission concept. 


1.7.2 AT- Sat Observatory 

The CSR concept provides a large rotating struc- 
ture, as well as large volume at the center hub, with 
the capability of operation above the majority of the 
atmosphere. Integrating high resolution, compact, 
linked optic sensors (such as the Low Mass Planar 
Photonic Imaging Sensor currently being developed 
as Phase II NIAC research) at the tip and root of 
each rotor- wing will permit a rotating imaging array 
across a full 600 diameter azimuth. Combining this 
with a lower resolution conventional imaging system 
in the center hub for image filling could provide a 
resolution never before possible. Such an observa- 
tory could be designed as a dual-purpose system that 
provides imaging both upwards for space investiga- 
tions, as well as for Earth imaging. A depiction of 
the device is shown in Figure 1.10. 

1.7.3 Compact VTOL Vehicle 

A version of the CSR concept has been identified 
that enables full retraction of the rotor-wing, shown 
in Figure 1.11 to achieve a highly compact VTOL 
vehicle that can operate at extremely low discload- 
ing (<0.03 lb/ft2) on atmospheres such as Mars 
that have atmospheric densities similar to Earth at 


100,000 feet altitude. Because it was highly uncer- 
tain whether the centrifugal stiffening loads would 
be sufficient to maintain a rigid rotor-wing for this 
highly faceted spanwise design, a sub-scale version 
of the fully retracting rotor-wing was fabricated and 
tested with success. The rotor-wing experiences 
an increase in weight and complexity with a decre- 
ment in performance to achieve compact packag- 
ing, however was shown to be a feasible approach 
that could facilitate very low discloading renewable 
VTOL aerial vehicles which could be packaged effec- 
tively in an aeroshell. 

1.7.4 Non- Commercial Applications 

A large number of other missions exists for the CSR 
concept with most relating to HALE missions on 
Earth. A comprehensive list is provided in Figure 
1.12. While clearly Facebook and Google will be- 
come major investors in this technology area due to 
the civil mission application potential, it is interest- 
ing to note how many Earth Science and Govern- 
ment Service mission are also of high importance to 
our societal needs. Because of this, it seems impor- 
tant that government research institutions such as 
NASA help to establish feasible HALE platforms, 
which can focus on non-commercial services. 
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High Altitude Long Endurance (@ 60,000 ft+) 


Earth Science Missions 

Extreme Weather Observation 
Hurricane Genesis and Evolution 
Tornado Monitoring 
Cloud and Aerosol Measurements 




Stratospheric Ozone Chemistry 
Tropospheric Pollution/Air Quality 
Water Pollution Monitoring 
Water Vapor Measurements 
02 and C02 Flux Measurements 
Vegetation Composition 
Aerosol and Precipitation Distribution 
Glacier and Ice Sheet Dynamics 
Antarctic Exploration Surveyor 
Imaging Spectroscopy 
Topographic LIDAR Mapping 
Magnetic Fields Measurement 
Oceanographic research 
Volcano Monitoring 
Surface Deformation Interferometry 


Government Services 

Active Plume Assessments 
Forest Fire Damage Assessment 
Forest Fire Communications 
Wildlife Census / Animal Tracking 
Land Resource Management 
Crop Disease Tracking 
Invasive Plant Assessment 
Pirate / Coastal Patrol 
Emergency On-Demand Comms 
Disaster Assessment 
Border Patrol 
Search and Rescue 
Mapping 


Civil Missions 

Omnipresent Wireless Broadband 
Sat-Cell Phone Service 
High Resolution Imagery 
Agriculture Yield Maximizing 


Figure 1.12: HALE missions that relate to CSR concept use. 




Chapter 2 

Plant Development 


This chapter expands the preliminary analysis performed on the CSR vehicle, and develops the 
plant model required for the controls development. It presents how the system is modeled and 
derives the equations of motion of the vehicle. These dynamics are implemented in a multi- 
body dynamics structure which governs the CSR simulations. Finally, the nonlinear system is 
trimmed and linearized which is the starting point for the controller development. 


Motivation 

The CSR concept vehicle is a radical departure 
from conventional fixed-wing and rotor craft design 
methodology, but this concept offers greater flight 
capabilities by combining the best features of both 
fixed-wing and rotor craft designs. However, conven- 
tional aircraft dynamics are not applicable to this 
concept vehicle, so a complete mathematical model 
for the CSR dynamics needs to be established. Sim- 
ilarly, existing control techniques are not appropri- 
ate for this concept vehicle. The primary objective 
with this initial analysis is to demonstrate that a 
suitable control law exists and can be successfully 
implemented on this vehicle. 

Proof of concept is demonstrated through scaled 
prototypes. Before developing the prototypes and 
full scale vehicle, simulation and analysis are re- 
quired to gain a better insight into how the system 
behaves. Before starting the controls analysis the 
system needs to be linearized so that control theory 
can be applied, and prior to linearizing the model, 
a set of initial conditions must be identified which 
put the system in a steady-state equilibrium condi- 
tion. Developing the equations of motion (EOM), 


finding trimmed control inputs, and obtaining a lin- 
earized model of the system is a crucial step towards 
beginning the controls analysis. 

Objectives 

The goal of this chapter is to develop a complete 
and accurate description of the CSR system. The de- 
sign methodology for this vehicle follows an iterative 
process, so basic models and prototypes are used as 
starting points, and more refined models and proto- 
types are developed over time. The models must pro- 
vide sufficient accuracy to evaluate the system, but 
also be modular enough to accept more sophisticated 
models throughout the development process. Differ- 
ent vehicle geometries at various scales are evalu- 
ated, so the models must be general enough to ac- 
commodate alternate configurations. Parameter files 
are used to store the vehicle geometry and physical 
attributes, so various designs can be evaluated by 
simply changing the parameter files. 

The first objective defines the system, develops 
the equations of motion (EOM), and implements 
the simulation of the CSR system. The system is 
reduced to a collection individual subcomponents. 
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Each component has a description, some physical 
parameters, and a coordinate frame. As the vehicle 
moves through space, various forces and moments 
act on these components. These external forces and 
moments are calculated for each individual entity at 
their current state, and then summed together to 
provide the total force and moment acting on that 
rigid body. All the relationships between compo- 
nents are compiled into a multibody dynamics struc- 
ture, which implements the numerical routine re- 
sponsible for the nonlinear simulation. 

The next goal of this chapter takes the developed 
CSR dynamics, and expresses those relationships in 
a form suitable for control analysis. Before applying 
controls, the system must operate in a steady-state 
equilibrium point, such as hover. Trimming is the 
process of finding states and control inputs which 
attain a particular equilibrium point. The trimming 
routine takes an educated guess value for the system 
initial conditions, and then optimizes a force balance 
equality for the lift generated by the system. Once 
the vehicle is trimmed, a linearization routine takes 
the nonlinear equations of motion and determines 
the state space matrices used for control analysis. 

Outline 

Before undertaking any of the controls development 
for the CSR system, a complete and accurate de- 
scription of the vehicle is required. This chapter de- 
velops all aspects of the physical plant used during 
the design and analysis process. 

Section 2.1 addresses Modeling , which illustrates 
the main features of the CSR vehicle, then decom- 
poses the system into individual components for the 
hub, tether, and satellite. It describes the physical 
properties of all the components, and assigns coor- 


dinate frames used during the modeling process. 

Section 2.2 derives the Dynamics of the CSR sys- 
tem. These are the mathematical relationships that 
describe the equations of motion of the vehicle. It 
identifies all the pertinent forces and moments act- 
ing on each component, and uses first principles to 
obtain those force and moment values. Finally, it 
presents a parallel CFD analysis which models the 
aerodynamics of the satellite body. 

Section 2.3 implements the Simulation that mod- 
els the vehicle within the Matlab environment. It 
formulates a multibody dynamics structure based on 
a numerical routine, presents a generalized approach 
to implement the routine, and applies the theory to 
the CSR models. 

Section 2.4 address Trimming the vehicle. This 
process seeks to obtain the set of states and con- 
trol inputs that achieve hover. While the CSR is in 
hover, the system is in a steady-state static equilib- 
rium condition, which is the starting point for lin- 
earization and controls development. 

Section 2.5 performs Linearization on the sys- 
tem. It takes the nonlinear system dynamics and 
yields linear relationships around the hover condi- 
tion. The values of the partial derivatives populate 
the state space matrices used for controls analysis. 

2.1 Modeling 

Modeling the CSR system is the first step towards 
deriving the system dynamics, building the simu- 
lations, and developing the control analysis. This 
section presents the overall vehicle structure, de- 
scribes the individual components, introduces coor- 
dinate frames, and illustrates the forces present on 
each component. 



Figure 2.1: Complete Centrifugally Stiffened Rotor (CSR) system. 
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Figure 2.2: Hub model. 


Figure 2.3: Tether model. 


2.1.1 System Overview 

Three tethered satellite vehicles rotate around a cen- 
tral hub, illustrated in Figure 2.1. Each satellite has 
a propeller mounted on a fuselage, a horizontal sta- 
bilizer with an elevator, a vertical stabilizer with a 
rudder, and a wing section. Centrifugal forces keep 
the tethers taught and lifting forces from the wings 
support each satellite and the central hub. Individ- 
ual satellites are controlled through their propeller, 
elevator, and rudder, and the central hub is con- 
trolled by manipulating tether forces coming from 
the three satellite vehicles. The axes for the inertial 
frame {N} are designated as: 

+X: North; +Y : East; + Z: Down. 

2.1.2 Hub Description 

The central hub rigid body anchors the three tether 
systems, illustrated in Figure 2.2. It has mass uih 
and moment inertia Jh = diag{ Jh x J Hy Jhz } ex- 
pressed in the hub frame {H}. The three tethers 
are arbitrarily labeled as A, B, and C, to distinguish 
them from each other. Position vectors ( ) 
relate the tether anchor points to the CG of the hub. 
Forces that act on the hub include gravity, tension 
from the tethers, and aerodynamic forces as the hub 
moves through space. The frame associated with the 
hub is located at the CG of the rigid body, with 
+X : Tether A; +Y : Right; +Z: Down. 

2.1.3 Tether Description 

Each tether is modeled as a chain of rigid body 
elements connected together with spherical joints, 
depicted in Figure 2.3. Each element has a cylin- 
drical shape with mass tut and moment inertia 
Jt — d±ag{ Jtx Jry Jtz} expressed in the tether 


frame {T}. Two position vectors, ri and r G , point 
from the CG oi the element to the inboard and out- 
board endpoints, respectively. Each tether element 
has four forces acting on it. Two forces, F{ and F OJ 
come from the joint constraints and are located on 
each end point of the rigid body. Gravity and aero- 
dynamic drag make up the other two forces, F g and 
Fy. The frame of each tether element is located at 
the CG of the rigid body, with 
+X: Outward; +Y : Right; +Z : Down. 

2.1.4 Satellite Description 

The satellites are attached to the tethers, and are 
the rigid bodies which contain all the control inputs 
to manipulate the CSR system. A satellite layout is 
provided in Figure 2.4. Each vehicle moves through 
the inertial frame as a rigid body, with mass ms 
and moment inertia Js = diag { Jsx Jsy Jsz} ex- 
pressed in the satellite frame {S}. External forces 
and moments that act on the vehicle come from grav- 
ity, the propulsion system, and aerodynamic effects. 
The aerodynamic forces and moments are generated 
from the wing section, horizontal stabilizer, vertical 
stabilizer, elevator, and rudder. The frame origin of 
the satellite is located at the CG of the vehicle, and 
its coordinate frame is oriented as 
+X: Nose; +Y : Right; + Z : Down. 



Figure 2.4: Satellite model. 
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Figure 2.5: COTS Prototype satellite vehicle. 


Figure 2.6: Matlab satellite vehicle. 


2.1.5 COTS Prototype Vehicle 

The first prototype in the spiral development life cy- 
cle is a commercial off the shelf (COTS) control line 
aircraft from Brodak, shown in Figure 2.5. This ve- 
hicle was selected as the first prototype platform be- 
cause it is small, light, inexpensive (< $100), and 
designed to fly out of the box. The simplicity of this 
prototype facilitates bench testing to configure new 
avionics hardware, test system identification tech- 
niques, and implement new control laws; all while 
reducing the risk of a crash from a custom proto- 
type with unknown dynamics. The vehicle weighs 

2.5 lbf (11.12 N), has a 38 in (0.97 m) wingspan, 
with a 24 in (0.61 m) fuselage length. 

2.1.6 Matlab Satellite Model 

The control line aircraft was modeled in Matlab 
to store vehicle geometries for simulation purposes. 
The Matlab satellite model is depicted in Figure 2.6. 
All the geometries are stored in a parameter file, and 


the satellite characteristics are calculated from those 
input dimensions. The satellite structure was as- 
sembled with modularity in mind, so parameter files 
can be changed to reflect different satellite layouts 
(different prototype versions, and various sized full 
scale systems) . The lifting force is located at quarter 
chord of the wing, and at the geometric center of all 
control surfaces. The motor vector includes tilt in 
pitch and yaw directions, if needed. 

2.1.7 Scaled Prototype Satellite 

The COTS control line aircraft tests new compo- 
nents and control algorithms. After successful initial 
demonstrations, testing advances to a scaled proto- 
type. This vehicle more closely resembles the full 
scale system, and collects data to determine the aero- 
dynamic characteristics of the CSR vehicle. The fol- 
lowing sections describe the geometry of the scaled 
and full sized satellites, and provides dimensions of 
the prototype vehicle. 



Figure 2.7: Isometric view of the prototype satellite. Figure 2.8: Side view of the prototype satellite. 
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Component 

Length/ 
Span (ft) 

Diameter/ 
Chord (ft) 

Shape/ 

Airfoil 

Fuselage 

5.5 

0.1667 

Circular 

Wing 

7.0 

0.5 

GOE596 

Horz Tail 

1.0 

0.5 

GE0596 

Vert Tail 

1.0 

0.5 

GE0596 


Table 2.1: Prototype parameters. 


Description 

The satellite has an asymmetric mono-wing design, 
inboard wing section, cylindrical fuselage, outboard 
horizontal stabilizer, and vertical stabilizer below the 
vehicle. The prototype layout is modeled in Vehicle 
Sketch Pad (VSP), where Figure 2.7 shows the iso- 
metric view of the satellite, and Figure 2.8 illustrates 
the side profile viewed from inboard to outboard. 
The figures depict the inboard wing in grey, the ver- 
tical stabilizer in blue, the outboard horizontal sta- 
bilizer in red, and the fuselage in tan. Dimensions 
of the prototype satellite are provided in Table 2.1. 

Wing Surface 

The prototype and full scale systems implement a 
unique twist distribution on the wing. Because the 
system is in a constant state of rotation, there is a 
nonuniform freestream velocity over the wing. The 
tapered twit geometry was designed to provide an 
optimized elliptical lift distribution over the lifting 
surface. The wing has a 7.0 ft (2.13 m) wing span, 
a constant chord length of 0.5 ft (0.15 m), and is 
located 0.8 ft (0.24 m) from the nose of the fuselage. 
Both the root and wingtip have zero incidence angle 
with a 6.5° linear washout that meet 3 ft (0.91 m) 
from the inboard wingtip. The wing uses a GOE596 
airfoil profile, displayed in Figure 2.9. 

Control Surfaces 

There is no twist or taper on the horizontal or ver- 
tical tail surfaces. The horizontal stabilizer is posi- 
tioned outboard of the fuselage. This location takes 
advantage of the wing tip vortices, which generates 
additional lift across the horizontal surface. The 
leading edge is located 1.75 ft (0.53 m) behind the 



Figure 2.9: GOE596 airfoil profile. 


nose of the fuselage, and 100% of the control surface 
is dedicated to pitch control. The vertical control 
surface extends below the fuselage. Camber is de- 
signed into this airfoil, so even with a zero angle of 
attack the surface generates lift in the outboard di- 
rection. The vertical stabilizer houses the landing 
gear, its leading edge is located 1.0 ft (0.30 m) from 
the nose of the fuselage, and 25% of the surface is 
used for yaw control. 

2.2 Dynamics 

Modeling defined the individual components that 
make up the CSR system, outlined relationships be- 
tween those components, and identified the forces 
acting on them. This section describes the dynam- 
ics of the system, which outlines the mathematical 
relationships for all the forces and moments acting 
on each individual component. This establishes the 
equations of motion that are used for both simula- 
tion and controls analysis. 

2.2.1 Overview 

The following sections outline the various forces and 
moments that act on each individual component of 
the CSR system. Then, the total force and moment 
on the rigid body element is the summation of all 
the forces and moments acting on that body. Dis- 
placement vectors define the locations where forces 
act on the rigid body, typically located at the CG, 


CHAPTER 2. PLANT DEVELOPMENT 


28 


or at the geometric center of a control surface. The 
forces and displacement vectors determine the mo- 
ments with the relationship 

M = d x F (2.1) 

where the cross product of the displacement vector 
d and the force vector F yield the resultant moment 
M for that subcomponent. 

2.2.2 Common Forces 

Many forces are common to all the CSR compo- 
nents. Such forces include gravity, centrifugal force, 
and aerodynamic lift/drag forces. Gravity is always 
downward (positive z in the inertial frame), and is 
expressed by 

F g = mg[0 0 1] T (2.2) 

where F g is the force of gravity, m is the mass of the 
element, and g = 9.81m/s 2 is the acceleration due 
to gravity. Nearly all CSR components are in a state 
of rotation which are subject to a centrifugal force. 
This force is expressed as 

F c = m(cj x (uo x r)) (2.3) 

where the centrifugal force vector F c is determined 
by the mass m, angular rate vector c j, and position 
vector r, of the element. Similarly, all rigid bod- 
ies experience aerodynamic drag, and many provide 
aerodynamic lift. Lift and drag forces are governed 
by the relationships 

F* = (l/2)pv 2 Ac*, * = l,d (2.4) 

where the lift Fi and drag Fj forces are dictated by 
the density of air p, velocity v, cross sectional area 
A , and the lift q and drag c d coefficients. 

2.2.3 Hub Forces 

Primary forces acting on the hub include: grav- 

ity, tension from the tethers, and aerodynamic drag. 
Gravity was previously described, and tether tension 
is accounted for within the simulation. Aerodynamic 
drag is a direct extension of (2.4), modeled as 

Fh = [ Fhx Fh v F Hz ] t , Fh* = c Hf *v 2 H * (2.5) 


for * = rr,y, z, which describes a linear relationship 
between hub drag force T#* and the square of the 
hub velocity vh *• A similar relationship exists be- 
tween hub drag moments Mh* and the square of the 
hub angular rates a;#*, described by 

Mh = [ M Hx Muy M Hz ] T , Mh* = ChtuMh* ( 2 - 6 ) 

for * = x, p, z. The coefficients, cuf * and c# m *, ma P 
the aerodynamic drag relationships for the hub. Ini- 
tial values for the coefficients are obtained by esti- 
mating terminal velocity and maximum torque on 
the hub. Aerodynamic drag force and moment are 
expressed in the hub frame and act at the hub CG. 

2.2.4 Tether Forces 

Each tether element feels five forces acting on it. 
Two forces, Fn and EV 0 , come from the inner and 
outer joint constraints, and are accounted for in the 
simulation structure. The force of gravity Ft 9 , and 
the centrifugal force Ft a have been previously de- 
scribed. The final force is aerodynamic drag Ftt 
which is modeled with (2.4). The tether moves 
through the inertial frame with velocity vt- Ac- 
counting for wind acting on the tether section, yields 
the relative airflow in the tether frame. The airflow 
vector component acts at the CG of the tether, and 
the component perpendicular to the tether gives the 
velocity v used in 

Frd = (l/ 2 )P V<2 Ac Td (2.7) 

where cxd — 0.47 is the drag coefficient around a 
cylindrical body. 

2.2.5 Satellite Forces 

The satellites contain all the available control in- 
puts, so controlling the CSR is accomplished en- 
tirely through these entities. Each satellite applies 
its own control locally, and then all the satellites are 
used collectively to control the complete CSR sys- 
tem. Therefore, an accurate assessment of the satel- 
lite behavior is needed for both vehicle simulations 
and controller design. 

Each satellite is modeled as subcomponents, 
which include: propeller, fuselage, wing, horizontal 
stabilizer with elevator, and vertical stabilizer with 
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a rudder. Each component has a designated frame 
and a position vector pointing from the frame origin 
to the CG of the satellite. Rotating between frames 
is done with an Euler rotation matrix R given as 


R = 


Cq C^jj 
CqS^ 

-so 


CfiS'ip H - $(j)S6C"ip 
C(f)C"ip H - S(f)SoS'ip 
S<f>C0 


— S^C^Pc^SqS^ 
C(f)C0 


(2.8) 


where cpf and cp m are the propeller coefficients. Fu- 
ture prototypes may include a canted motor, so 

F m = R( 0 , 9m, i>M )Fp ( 2 . 12 ) 

maps the force from the propeller frame to the motor 
frame. Similarly, the moment in the motor frame is 
described by 

Mm = dM x Fm + R( 0, 9m, ^ m )Mp (2.13) 


where the attitude (roll, pitch, and yaw) is given by 
(</>, 0, / 0), and s* and c* are sine and cosine functions. 

The force of gravity and the centrifugal force 
have already been presented, so the following sec- 
tions describe the remaining forces acting on the 
satellite. These include the propulsion from the pro- 
peller, skin friction drag from the fuselage, and the 
lift and drag profiles and how they are applied to the 
aerodynamic surfaces. Any forces that do not act 
through the CG of the satellite will have a moment 
described by (2.1), and the total force and moment 
on the satellite is the summation of all individual 
forces and moments. 


Propulsion Model 

The motor spins at angular rate £lp, and generates 
forces and moments described by 

F P = [Fp x Fp y Fp z ]T Mp = [M Px Mp y Mp z ]T 

(2.9) 

The inner workings of a brushless DC motor, and 
the aerodynamic intricacies of a fixed pitch propeller 
may be handled in future developments. For the ini- 
tial analysis, a simplified propulsion model is utilized 
as a first approximation for the propulsion dynamics. 
Blade flapping and side airflow are neglected, so the 
force and moment of the propeller only act parallel 
to the z-axis of rotation; therefore, (2.9) reduces to 

Fp = [0 0F Pz ]T M p = [0 0M Pz } t (2.10) 

The simplified model contains a linear relationship 
between the squared rotor angular rate Lip and the 
force and moment (Fp z , Mp z ) delivered from the 
propulsion system. This relationship is expressed as 

( 2 . 11 ) 


where dM is the position vector from the CG of the 
satellite to the motor frame. 


Skin Friction 


As the satellite moves through space, the fuselage is 
subjected to skin friction; a type of drag force from 
the airflow grabbing the fuselage surface. The skin 
friction drag coefficient is defined as 


CFd = 


hpul 


(2.14) 


where r w is the wall shear stress, p is the density 
of air, and Uqq is the free stream velocity. Treating 
Uoo = v yields 

FFd = (l / 2 )P v2 Ac F d 


— (1/2 )pv 2 A 




(2.15) 


= At,, 


where r w is defined as 


(2.16) 
y = o 

and where p is the dynamic viscosity air, v is the 
fluid velocity along the boundary, and y is the dis- 
tance above the boundary layer. 


r w = V 


dv 

dy 


Lift and Drag Profiles 

Several prototypes are planned as part of the spiral 
development design process, and different wing ge- 
ometries will be evaluated during these design itera- 
tions. To easily accommodate these revisions in the 
simulation, a Matlab function models the lift and 
drag profiles for each of the control surfaces. The 
simulation stores lift/drag attributes including: zero 
angle q and Cd, linear range and slope, and max/min 
cl values. Then the program uses these parameters, 
to assemble the lift and drag coefficients as functions 


Fp z = cpfCtp , M Pz = cp m Dp, 
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of the relative angle of attack of the surface. The lift 
and drag functions are depicted in Figure 2.10. 

Wing Lift and Drag Profiles 




Figure 2.10: Lift and drag profile curves. 


Aerodynamic Surfaces 

The five surfaces on the satellite include the wing, 
horizontal stabilizer with an elevator, and vertical 
stabilizer with a rudder. Each surface has its own lift 
and drag profile to accurately represent each compo- 
nent. The satellite moves through the inertial frame 
with a known velocity. Accounting for wind within 
the simulation, yields an airspeed expressed in the 
inertial frame. The attitude of the satellite forms 
a rotation matrix that converts the airspeed vector 
to the local satellite frame. For each surface, the 
the airspeed vector is projected onto a plane parallel 
with the cross section of the surface. This projection 
yields the relative angle of attack and the magnitude 
of the airflow over the surface. The relative angle of 
attack is fed into the lift and drag profile function, 
which returns the lift and drag coefficients (cl,cd) 
for that particular angle and surface. The magni- 
tude of the airspeed v, together with cl or cd, are 
passed into (2.4) which yields the magnitudes of the 
lift force and drag force. These magnitudes are con- 
verted into vectors where the drag force acts parallel 
to the airflow, and the lift force acts perpendicular 
to the airflow. The force vectors act at the quar- 
ter chord of the wing surface, and at the geometric 
center for all other surfaces. Each surface has a po- 
sition vector pointing from the CG of the satellite to 
the point where the aerodynamic force acts. Then 
individual moments are calculated with (2.1). 


2.2.6 CFD Aerodynamic Model 

Satellite aerodynamics account for the largest con- 
tribution towards the CSR equations of motion, so a 
high fidelity aerodynamic model is desired to accu- 
rately simulate the vehicle dynamics. The first ap- 
proach for modeling the satellite relied on first prin- 
ciples that govern rigid bodies and lifting surfaces. A 
parallel effort models the satellite vehicle using com- 
putational fluid dynamic (CFD) techniques. The fol- 
lowing sections describe the software, how the satel- 
lite is modeled, the operation point, and the look-up 
table used within the Matlab simulation. 

Athena Vortex Lattice 

Satellite aerodynamics are obtained with Athena 
Vortex Lattice (AVL) software. The AVL program 
performs aerodynamic analysis on rigid aircraft with 
arbitrary configurations. It employs an extended lat- 
tice vortex model for the lifting surfaces and a slen- 
der body model for the fuselage. The software starts 
with the nonlinear system and generates a full lin- 
earization around any flight condition. The satel- 
lite is parsed into surfaces and fuselage components, 
which are input into the AVL program. Surfaces are 
modeled as a series of attached sections, where each 
section is described by a chord, leading edge loca- 
tion, airfoil, and pitch angle. Control surfaces have 
additional parameters which include a gain, hinge 
location, hinge vector, and sign of deflection. Co- 
ordinates of the top and bottom edges define the 
shape of the fuselage. Finally, all the user-defined 
sections are joined together through a linear interpo- 
lation process which generates the complete vehicle 
model. Figure 2.11 shows the AVL representation of 
the CSR satellite. 
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Modeling Elliptical Lift Distribution 

The most important goal of any modeling effort is 
to generate a model that accurately represents the 
physical system in the real world. The CSR satel- 
lite is designed to have an elliptical lift distribution 
across the wing surface. This aerodynamic effect is 
an important contributing factor for eternal flight, 
and should be reflected in the aerodynamic mod- 
els. Because the satellite rotates around the cen- 
tral hub, outboard wing sections have faster velocity 
profiles than inboard sections. The tapered wing 
is designed to account for this differential velocity 
while attaining the desired elliptical loading. Unfor- 
tunately, AVL does not permit differential velocities 
within the program. To work around this issue, the 
twist in the main wing was removed. This adjusted 
model in a uniform flow field, has a similar elliptical 
lift distribution as the CSR satellite operating in a 
state of rotation. 

Operating Conditions 

Differing scales of the CSR vehicle will have differ- 
ent operating points, which are dependent on the 
vehicle rotational velocity and the distance between 
the satellite and the hub. Similarly, each system will 
have two different operating points for hover and for 
takeoffs/ landings. During start up, the tether is re- 
tracted so the satellite must spin more quickly to 
generate sufficient lift. With a retracted diameter of 
22 ft (6.71 m), the vehicle must rotate at 74 rpm 
(7.75 rad/sec) for takeoffs and landings. Once the 
vehicle is airborne, the tethers are let out and the 
satellites generate the necessary lift at lower rota- 
tional velocities. When the tether is extended to 
100 ft (30.48 m), the hub must spin at 11 rpm (1.15 
rad/sec) to maintain the satellite cruise speed. Dur- 
ing hover, the linear differential velocity across satel- 
lite wing, from inboard to outboard, is 29 mph (12.96 
m/s) to 34 mph (15.20 m/s), respectively. 

Look-up Table 

When the AVL software is run, a unique operat- 
ing point must be selected, then the software pro- 
duces a linearized model around that point. How- 
ever, linearization is only accurate for a small re- 
gion around the point. To increase the fidelity of 
the model, a lookup table is populated with trends 


obtained by running the software at a variety of op- 
erating points. These operating conditions span the 
entire range of the flight envelope and actuator lim- 
its, which provide a complete dynamic picture of the 
satellite. Five parameters are adjusted 

• angle of attack, 

• side-slip angle, 

• velocity magnitude, 

• rudder deflection, and 

• elevator deflection 

and forces and moments are stored for each oper- 
ating point. Then, a 5D linear interpolation pro- 
cess approximates the forces and moments between 
points. The AVL program is automated to generate 
the aerodynamic lookup table, where the boundaries 
of the variables are shown in Table 2.2. 


Variable 

Min 

Max 

Points 

AoA (<a) [deg] 

-12 

12 

7 

Side-slip (/?) [deg] 

-12 

12 

7 

Rudder (Sr) [deg] 

-25 

25 

7 

Elevator (Se) [deg] 

-25 

25 

7 

Velocity [m/s] 

30 

15 

8 


Table 2.2: Aerodynamic boundaries. 


Trends 

This section presents a series of figures that show- 
case the AVL generated results. For all plots, angles 
are in degrees, forces are in Newtons, and moments 
are in Newton-meters. Lift and drag are primarily 
a function of angle of attack a and vehicle airspeed 
V, so the first four plots demonstrate these relation- 
ships in Figure 2.12. The next six plots present the 
change in force and moment with respect to a change 
in a single independent variable. The nominal air- 
speed is 22.5 m/s, and the angle of attack, side-slip 
angle, rudder deflection, and elevator deflection are 
all set to zero degrees. Each of the five subplots 
changes a single variable, while all other variables 
are held at their nominal values. Figures 2.13, 2.14, 
and 2.15 showcase the forces, and Figures 2.16, 2.17 
and 2.18 present the moments. All the plots demon- 
strate well-behaved trends which are smooth func- 
tions with no discontinuities. 
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Figure 2.12: Lift and drag forces as functions of angle of attack and airspeed. 


2.2.7 Future Revisions 

The CSR development is an iterative process, and 
the presented work only represents the first model 
development. Several simplifying assumptions are 
present which will be addressed in future work. 

Current derivations treat each satellite vehicle as 
a rigid body. In reality, the long and slender wing 
section may experience bend and/or twist during 
flight. Subsequent developments can model these ad- 
ditional dynamics through a Finite Element Analysis 
(FEA) or a more refined multibody dynamics model. 

The existing model only evaluates individual 
satellite components, so it does not account for 
aerodynamic interactions between the components. 
Items such as wing tip vortices, aerodynamics 
around the wing and fuselage connection, and inter- 
actions between the stabilizers and control surfaces, 
are not accounted for. Future work will continue 
to develop the CFD model of the satellite to popu- 


late the Matlab look-up table. This data will return 
forces and moments acting on the satellite for given 
airspeeds and control inputs. 

Finally, a basic linear relationship models the 
propeller dynamics. More refined models should ac- 
count for off-axis forces and moments, side airflow 
into the propeller, blade flapping, and time constants 
within the motor system dynamics. 

2.3 Simulation 

With the CSR dynamics established, the final step 
toward implementing the equations of motion, es- 
tablishes the relationships between all the rigid body 
elements. This section formulates the multibody dy- 
namics problem through an illustrative generic ex- 
ample, expands the theory to the Round the Pole 
(RTP) and CSR systems, and explains how the nu- 
merical routine is implemented in Matlab. 
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Figure 2.13: Longitudinal force trend. Figure 2.16: Longitudinal moment trend. 



Figure 2.14: Lateral force trend. 


Figure 2.17: Lateral moment trend. 



Figure 2.15: Vertical force trend. 


Figure 2.18: Vertical moment trend. 
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2.3.1 Multibody Dynamics Motivation 

Simple elements, like springs and dampers, have an 
exact representation between forces and states, so a 
closed form solution is available. With the CSR dy- 
namics, airflow over the tether generates drag forces 
which are dependent on the current state of the 
tether, so no closed form solution exists. Therefore, 
a numerical approach is required to represent the 
tether in the CSR model. 

Several approaches are available to model a 
multibody system. Hand derivations like Kane’s and 
Lagrange’s methods, are laborious and changes to 
the model or number of elements require new deriva- 
tions. The following approach is similar to an Finite 
Element Analysis (FEA) model where each section 
is modeled as a smaller building block, but FEA lin- 
earizes early in the process, and is only suitable for 
small perturbations. 

The CSR simulation implements a multibody dy- 
namics formulation. This approach keeps the non- 
linearity in the model, is appropriate for numeri- 
cal integration, and can be assembled by inspection. 
This approach offers the most accurate model, and 
is modular enough to accommodate design revisions 
and various prototypes during the design process. 

2.3.2 P reliminar ies 

The system is broken down into rigid body subcom- 
ponents connected together with spherical joints. 
The simulation structure represents some states in 
the inertial frame and other states in a local body 
frame. This section describes the dynamics of a 
rigid body, the rotation matrix between coordinate 
frames, how to express derivatives in different coor- 
dinate frames, and a common notation used within 
the multibody development. 

Rigid Body Dynamics 

The dynamics of a rigid body under external forces 
and moments is governed by 

mv — EF, Jd) + uo x Joj = EM (2-17) 

where EF G M 3 and EM G M 3 represent external 
forces and moments acting on the body, v is the ve- 
locity vector in the inertial frame, uo is the angular 


rate vector, and m and J are the mass and moment 
inertia of the rigid body. 


Rotation Matrix 


During the development of the equations of motion, 
some elements need to be transformed between the 
inertial frame and its local body frame. When work- 
ing with Euler angles, the rotation from body frame 
to inertial frame is given by the rotation matrix F, 
described by 


R = 


CqS^Jj S^SQS^p -pC^SQS^p 

-So S<j)CQ C(f)CQ 


(2.18) 


where sine and cosine ( 5 *, c*) act on the attitude 
(0,0,^) of the element. Rotation matrices possess 
a special property where the inverse is equal to the 
transpose, F _1 = F T , so mapping from the inertial 
frame to the body frame is accomplished with F T . 


Derivatives 

Differentiation is done with respect to a local frame. 
The relationship 

nL v ) = b Lv) + n ui b xv = v + n u b xv (2.19) 

U/l/ U/l/ 

equates time derivatives in the body frame B to the 
inertial frame N. 


Notation 


Introduce two operators to keep the notation clean. 
The cross product of two vectors is equivalent to the 
matrix multiplication of a skew symmetric matrix 
and a vector, so that xy = x x y, where 


0 

~%3 

X 2 

£3 

0 

-xi 

-X 2 

Xi 

0 


(2.20) 


is the skew symmetric matrix of x. Similarly, a dou- 
ble cross product can be represented in matrix form, 
so that xy — x x (x x y), where 

I" — x\ — x\ X\X 2 X\X 3 | 


x — 


X\X 2 ~x\ — x\ X 2 X 3 


(2.21) 


|_ XlX 3 X2X3 -x\ -x\ \ 

is the matrix equivalent double cross product of x. 
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' Ji 
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Jo 
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Vo 


Fo 

— R* r T 
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R° r T 
n N' o 

~h 

h 

0 


_ Fq _ 


R l N O iri - R° N LV 0 r 0 


(2.22) 



2.3.3 Two Rigid Bodies 

The simplest case to consider is a two body system 
connected with a spherical joint, which is illustrated 
in Figure 2.19. There is an inner body Bi and outer 
body B 0 connected at joint G, which has a joint 
force Fq and joint torque Tq. The joint force and 
torque are arbitrarily designated as acting in a posi- 
tive sense on the inner body. Position vectors and 
r Q extend from each CG to the connection point, and 
both bodies have a force, F\ and F 0 , and torque, Ti 
and T 0 , acting on them. 

Euler’s equation describes the two relationships 
for rotation as 

— Hi W. i X Hi + Tq + T{ X Fq , ^ 

J o L0 o = To - U Q x Ho - Tq - r Q X Fq 

where J* is the moment of inertia, is the angular 
rate vector, and H* is the angular momentum of the 
rigid body. When there are no momentum storage 
devices, the angular momentum term is expressed as 
H — Juo. Newton’s Law gives the relationships for 
translational as 

rriiVi = Fi + Fq, m 0 v 0 = F 0 - F G (2.24) 



where m* is the mass, and v* is the velocity vector 
of the rigid body. A fifth relationship models the ve- 
locity of the joint constraint between the two bodies 

v G = Vi + Ui x Vi = Vo + LOo x r Q (2.25) 

and after differentiating within the inertial frame N 
yields the relationship 

Vi + iOi x n + LViTi = Vo + LVo x r Q + W 0 r 0 . (2.26) 

The five equations are put into matrix notation, and 
rotation matrices are added to convert between co- 
ordinate frames. The result is shown in (2.22). 

2.3.4 N Rigid Bodies 

The previous section laid the ground work for form- 
ing the equations of motion. This section expands 
upon that foundation and generalizes the equations 
of motion for N rigid bodies, where a sample four 
body system is depicted in Figure 2.20. A connec- 
tivity table identifies the inner and outer bodies of 
each joint, presented in Table 2.3. Each joint con- 
nects two bodies, where i and o designate the inner 
and outer body, respectively. 
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G1 

G2 

G3 

B1 

i 



B2 

0 

i 

i 

B3 


0 


B4 



0 


Table 2.3: Four body connectivity table. 


The equations of motion are formed through in- 
spection using the connectivity table. The matrix 
relationship for an N rigid body system is given by 


J OR' 


(jj 


T 

0 MI 


V 

= 

7 

1 

O 

b 

h 

0$ 


_ F _ 


_ C _ 


(2.27) 


The states of the system are stored in the first col- 
umn vector, where uo G R 3Ar contains the angular 
velocities, v G R 3Ar contains the linear velocities, 
and F G R 3 ^ -1 ) contains the joint constraint forces. 
The right column vector stores the external torques 
in r G R 3iV , the external forces in T G R 3Ar , and the 
constraint relationships in C G M 3 ^ -1 ). The left ma- 
trix is populated with sub matrices. Matrices J and 
M are each M 3iVx3iv diagonal matrices, that contain 
the masses and moment inertias for each rigid body. 
Matrices R and / are both ]j 3Ar x3(7V-i)^ and follow 
the same structure as the connectivity table. 

Matrix / has negative 3x3 identity matrices for 
inner body connections, and positive 3x3 identity 
matrices for outer body connections, formed as 


I 3 0 0 

h ~h ~h 
0 h 0 
0 0/ 3 


(2.28) 


Matrix R follows an identical sign convention, but is 
populated with position vectors and rotation matri- 
ces, described by 


' -hiRi 0 0 

hi R2 -rnR? 

0 f 32J Rf 0 

0 0 r 43 i?f 


(2.29) 


This concludes the generalized procedure for de- 
termining the equations of motion of a multibody 
system. The next sections apply this procedure to 
generate the dynamics for two different CSR models. 


2.3.5 Round the Pole Model 

Part of the CSR development will build a scaled 
model prototype of the system. Before the complete 
system is flown in free flight, several development 
stages will collect data from a single satellite vehicle 
flying around a fixed pole. This collected data will 
characterize the satellite and tether model, and flight 
data compared to simulation data will validate and 
refine the simulation model. Therefore, a model of 
the single satellite system is needed for this future 
work. This section describes the multibody struc- 
ture of the Round the Pole (RTP) system. 

The system layout is shown in Figure 2.21. The 
multibody structure requires one more rigid body 
than the number of constraints, so the hub H must 
be included in the model. The hub should not move 
with respect to the inertial frame, so it is modeled 
as a rigid body with an extremely large mass and 
moment of inertia, and has no gravity acting on the 
body. This provides a very good model for a pole 
fixed in the inertial frame. 

The model contains a single chain of rigid bod- 
ies that start with tether link Tf, work outward to 
tether link T n , and finish with the satellite S as the 
final chain link. The connectivity table has a very 
simple structure, provided in Table 2.4, which is used 
to form matrices R and / as described previously. 

With n tether elements, there are N = n+2 rigid 
bodies, and n + 1 constraints. Each tether element 
is identical, so they all have the same mass tut and 
moment inertia Jt . Sub matrices J and M are both 
jg)3Wx3W ^agonal matrices, where the first three en- 
tries contain the hub properties (ra# and Jh), the 
last three entries contain the satellite properties (ms 
and Js), and the middle 3 n entries contain the tether 
properties (my and Jt) repeated n times. 



Figure 2.21: Round the Pole (RTP) diagram. 
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Table 2.4: RTP connectivity table. Figure 2.22: Free flight CSR model. 


2.3.6 Full System Model 

The other model describes the full CSR system. It 
includes the central hub, three tether arms with 
satellites, and no external constraints on the system. 
The system layout is depicted in Figure 2.22. 

Several connectivity tables are possible. The 
most desirable configuration was selected so simula- 
tion code could easily scale to accommodate different 
numbers of tether elements, or additional rigid bod- 
ies to model flexibility in the wing section. The con- 
nectivity table for the free flight system is shown in 
Table 2.5, and the other sub matrices are populated 
in a similar way to the RTP flight model, where the 
first elements store the hub parameters, the last el- 
ements store the satellite parameters, and the inner 
elements store the tether parameters. 
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Table 2.5: CSR connectivity table. 


2.3.7 Future Work 

The multibody structure offers a nonlinear model 
that scales easily and is assembled through inspec- 
tion. The priority was to develop the dynamics as 
easily as possible, so controls analysis could begin as 
soon as possible. However, the current multibody 
structure has some potential areas of improvement 
which may be addressed in future work. 

The current structure completely models each 
rigid body and all constraints among them. Many 
relationships can be expressed in terms of states of 
neighboring elements. Thus, the current structure 
models redundant states and does not offer a mini- 
mal set. Because the simulation is based on a numer- 
ical routine, small numerical discrepancies between 
these redundant states may accumulate over time. 

The simulation runs in Matlab with an ordinary 
differential equation (ODE) solver. Solving for the 
derivatives of the states requires the inverse of the 
square matrix given in (2.27). This is not a concern 
for a smaller number of rigid bodies, but the size of 
the matrix grows rapidly with higher fidelity mod- 
els (more elements are added). The structure of the 
matrix is predominantly diagonal, so a more efficient 
inverse procedure may be devised. 

2.4 Trimming 

Dynamics establishes how the system moves, and 
the multibody structure implements those equations 
of motion in the simulation. This section addresses 
trimming, which seeks the set of states and control 
inputs that achieve an equilibrium point. In other 
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words, the set of values that keep the system states 
and control inputs constant. 

Trimming is an important milestone for several 
reasons. First, it obtains the values for specific flight 
conditions, like hover and forward flight. These 
steady state values are also needed to perform lin- 
earization, which is required to begin the controls 
analysis. Finally, trimmed states are used as system 
initial conditions before running a simulation. 

The CSR has a complicated nonlinear model 
with intertwined states and control inputs, so an an- 
alytic approach to determining equilibrium points is 
not possible. Matlab has built-in functions which 
can trim certain classes of problems; however, the 
special structure of the CSR multibody problem pre- 
cludes the use of these existing functions. This sec- 
tion presents the procedure for trimming the CSR 
system, addresses limitations with traditional trim- 
ming routines applied to multibody dynamics prob- 
lems, and presents an alternative approach using 
nonlinear optimization techniques. 

2.4.1 Problem Formulation 

The objective is to find the system states x(t) and 
control inputs u(t) to attain an equilibrium point, 
which means the system is not accelerating. Achiev- 
ing zero acceleration requires a force balance within 
all the elements. Thus, the goal is to find values 
where the forces and moments sum to zero for every 
rigid body in the system. Following sections illus- 
trate deficiencies with current trimming techniques, 
and outlines a solution to overcome the problem. 



Figure 2.23: Typical model to trim. 



Figure 2.24: Multibody model to trim. 


Multibody Trimming Limitation 

Many types of modeling problems are formulated as 
decoupled rigid bodies which are connected together 
through springs and dampers. Figure 2.23 illustrates 
this common modeling form. For this arrangement, 
the masses (and their states) are independent of one 
another, so a single state can be perturbed while 
the others remain static. In the multibody dynam- 
ics structure, the rigid bodies must maintain joint 
constraints, so the states of the rigid bodies are not 
independent. Perturbing a single rigid body element 
leads to a disjoined chain structure, depicted in Fig- 
ure 2.24. This is not a viable solution, so an alter- 
native trimming procedure alleviates this concern. 

Optimization Function 

Trimming is a type of optimization problem. The 
Matlab function, lsqnonlin, solves nonlinear 
least-squares curve fitting problems. The input is 
a vector x with initial condition xo, and a vector as 
a function of x is the output 

fix) = [ fi(x) f 2 (x) ■ ■ ■ f n (x) } T . (2.30) 

The optimization routine solves 

min-ll/CO 111 =min(/i(a;) 2 +f 2 (xf 4 (2.31) 

X X X 7 

which minimizes the square of the Euclidean norm 
of the vector function. 

Process Overview 

The trimming process is accomplished in two loops: 
an inner loop and an outer loop. The inner loop 
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starts at a tether anchor point and works outwards 
toward the satellite. Optimization runs on each rigid 
body to achieve a force balance within that element. 
After the satellite has been optimized, the entire 
tether arm achieves a force balance condition. How- 
ever, the lift generated by the satellite depends en- 
tirely on the initial values used at the hub, so there 
is no assurance that the generated lift equals the 
required lift. The outer optimization loop incre- 
mentally adjusts the hub rotation and anchor point 
forces, runs the inner loop again, and reevaluates the 
amount of lift generated. The outer loop optimiza- 
tion continues until there is a force balance within 
the central hub. 

2.4.2 Optimization 

The trimming routine uses several different opti- 
mizations to to obtain a force balance within the 
elements. The first optimization balances a single 
tether element. Then a loop applies this routine 
to all the links in the tether. Another optimization 
finds the satellite states and control inputs to achieve 
a force balance within the outermost segment of the 
tether chain. This process completes the inner loop 
routine. The final optimization seeks a force balance 
within the central hub. This outer loop optimization 
adjusts the hub angular velocity and tether anchor 
force until the lift generated by the satellite equals 
the weight it needs to support. The following sec- 
tions describe each process in detail. 

Tether Force Balance 

Consider a single tether element. The goal is to find 
the state that achieves a force balance between all 
the forces acting on the element. These forces in- 
clude the inner and outer joint constraint forces, F{ 
and F 0) a centrifugal force F c , aerodynamic drag i^, 
and gravity F g . Assume that the hub angular rate 
and the tether’s inboard position, velocity, and force 
are all known. The attitude of the element dictates 
the CG location and the direction of relative airflow 
over the tether section. These states determine the 
drag force, centrifugal force, and outboard endpoint 
force. The optimization adjusts the attitude of the 


tether until the force balance equality 

Fi + F c + F d + F g + F 0 = 0 (2.32) 

is satisfied. A similar relationship balances the mo- 
ments acting on the tether element. 

Tether Projection 

The previous section outlined a process that takes 
an inboard tether state, balances the forces within 
the element, and returns the outboard tether state. 
Tether projection describes a program loop that bal- 
ances all the elements of the tether. The loop starts 
with the hub anchor force and the angular rate of 
the hub, and applies the optimization to the inner 
most tether element. The outboard states of this 
first element become the inboard states for the next 
element. The loop proceeds outward through each 
section until it reaches the final tether element. The 
outboard states of the final tether are used as the 
anchor position, velocity and force for the satellite. 

Satellite Force Balance 

The satellite optimization is very similar to the 
tether optimization. The tether end point states 
were determined from the tether projection loop, but 
the satellite attitude and control inputs are still un- 
known. As a starting point, the satellite is initial- 
ized with a zero attitude, zero control surface de- 
flection, and a reasonable guess value for throttle 
input. The optimization takes the tether end point 
states, and the current values for satellite attitude 
and control inputs. For this particular state, the 
program evaluates the centrifugal force and aerody- 
namic forces/moments, calculates the summation of 
forces and moments on the satellite, and compares 
this value to zero. The optimization incrementally 
adjusts the satellite states and control inputs until 
the net forces and moments are reduced to zero. This 
completes the inner loop. 

Hub Force Balance 

The previous sections describe the inner loop which 
balances the entire tether arm with a satellite. That 
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process began with an initial guess for the hub an- 
chor force and system angular rate. Using these ini- 
tial guess values, the lift generated from the satellite 
may not equal the lift force needed for the system. 
Therefore, an outer loop optimization is needed to 
balance the hub forces. The optimization compares 
the current amount of lift to the force required to 
support the weight of the entire system. It also bal- 
ances the moment generated by the tether force with 
the rotational aerodynamic drag from the spinning 
hub. If the net forces and moments are non zero, the 
optimization routine adjusts the hub angular rate 
and the hub anchor force, and then runs another in- 
ner loop optimization. This iterative process contin- 
ues until the hub and satellite achieve static equilib- 
rium. Once the CSR system is completely balanced, 
the states and control inputs are stored as initial con- 
ditions for the simulation and as operating points for 
linearization. 

2.4.3 Initial Guess Values 

All the minimization routines require initial con- 
ditions to begin the optimization process. Values 
for the inner loop optimizations (tether and satel- 
lite) were initialized locally. However, the outer 
loop initial conditions (hub angular rate and an- 
chor force) are still unaccounted for. This section 
presents derivations that yield approximate values 
suitable for initializing the outer loop optimization. 
It addresses the satellite lift force, linear velocity, 
and position, which is used to estimate the system 
angular rate and anchor forces. 

Lift Force 

The masses of the elements are known, so the first 
step finds the total weight of the system. There are 
three tether arms, so the lifting force required of each 
satellite must carry one third of the total weight of 
the system. The program stores the required lift 
force Fg which is used in the next section. 

Satellite Velocity 

The lift equation is given by 

F S = (1/2 )pv 2 s Ac S L (2.33) 


where Fg is the satellite lift force previously calcu- 
lated, p is the density of air, v§ is the forward linear 
velocity of the satellite, A is the planform area of the 
wing, and cgL is the lift coefficient of the satellite 
wing. The only unknown value in this relationship 
is the satellite velocity. Rearranging the equation 
leads to the equality 

vs = \fM~L {2M) 

which yields the forward velocity vg required to 
achieve the desired lift force Fg. This result as- 
sumes a level flight condition. The true velocity will 
be slightly greater to because the trimmed satellite 
attitude will not be perfectly level. 

Satellite Position 

During flight the tether will curve from aerodynamic 
drag and from the force of gravity. The amount of 
curvature is unknown until after the optimization 
routine. As a starting point, assume the tether has 
no curvature and projects straight out towards the 
satellite. The hub anchor position, the length of the 
tether, and the satellite anchor position have known 
dimensions, which can be added together to obtain 
rg, an approximate starting position for the satellite. 

Satellite Angular Rate 

With an approximate satellite position and velocity, 
the next step estimates the angular rate ujh of the 
central hub. Relating the satellite position rg and 
velocity vg with 

uj H = — • (2.35) 

r s 

yields the estimated angular velocity uh- 

Anchor Force 

The final set of calculations estimate the force vec- 
tor acting at the hub anchor points. The estimated 
satellite position rg and the hub angular rate ujh 
determines the centrifugal force F x described by 

F x = msrsojn 


(2.36) 
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Figure 2.25: Various views of the trimmed RTP system. 


where ms is the mass of the satellite. Rotational 
aerodynamic drag dictates F y with 

F y = - (2.37) 

3 T 

where Mh z — CHmz^H Z the ^~ ax i s moment on the 
hub for a particular angular rate, r is the distance 
from the hub CG to the tether anchor point, and the 
one third term equally distributes the force to each 
of the satellite anchor points. The vertical force F z 
is influenced by gravity g and the mass of the hub 
rriH described by 

F z = -(1/3 )m H g (2.38) 

where the one third distributes the total force to each 
tether arm, and the negative sign remains consistent 
with the coordinate frame. 

2.4.4 Plots of Trimmed Systems 

The previous sections described the process for trim- 
ming the RTP and CSR systems. This section 


presents the results of that effort by displaying two 
figures of the trimmed systems. 

Figure 2.25 depicts several perspective views of 
the trimmed RTP system. Tether parameters (den- 
sity, diameter, mass) were set to larger values to ex- 
aggerate the effects of drag and gravity, and to show- 
case the trimming routine. Notice, that the tether 
bends uniformly in the vertical direction from grav- 
ity, and the horizontal curvature becomes more pro- 
nounced outward as drag increases. Furthermore, 
the CG of the satellite aligns with the tether end- 
point force vector, as expected. 

Figure 2.26 displays the states of a trimmed CSR 
satellite in the inertial frame. The satellite flies a 
circle around the hub, so the horizontal position and 
velocity are sinusoidal. The central hub rotates at 
a constant velocity, as seen in the upper right cor- 
ner. All other states are in static equilibrium, so 
they maintain very nearly zero values. 
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Figure 2.26: Trimmed CSR satellite states in the inertial frame. 


2.5 Linearization 

Thus far, the CSR plant has trimmed operating 
points and nonlinear dynamics implemented in the 
simulation. The last step toward developing the CRS 
plant, takes the nonlinear system and creates a linear 
model around those trimmed operating points. 

Linearization is an important topic because the 
vast majority of the controls analysis is based on the 
linear system model. Linearizing multibody dynam- 
ics problems have many of the same hurdles as with 
trimming, so this research utilizes the alternative 
trimming technique and adapts it to the lineariza- 
tion process for multibody dynamics problems. 

This section presents the Taylor series expansion 
of a nonlinear model, and describes how the partial 
derivatives form state space matrices of the linear 
system. It addresses issues with linearizing multi- 
body dynamics problems, proposes an alternative 


approach, and presents the matrices that represent 
the Matlab simulated COTS control line aircraft. 


2.5.1 Taylor Series Expansion 


The equations of motion of the plant are expressed 
in compact form as x(t) — f(x{t),u{t )) where x{t) G 
W 1 are the n states of the plant and u{t) G M m are 
the m control inputs. Denote an operating point 
as x and u , so that perturbations are described by 
Ax — x — x and A u — u — u. 

Taylor series expansion for the nonlinear system 
is described by 


df 


X = f(x,u ) = f(x,u ) + — 


du 


(x,u) 


(x — x) 

(x,u) (2.39) 

(u — u) + HOT 


and disregarding the higher order terms (HOT), 
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yields the form 

Aa; = f(x, u) + AAx + BAu (2.40) 


where the state matrix A is described by 


A = 


- dh dh dh - 

dx\ dx2 ' ’ * dx n 

d% 

dx± 


dfn dfn 

_ dxi * * * dx n _ 


and the input matrix B is given by 


B = 


dh dh dh ~ 

dui du2 ' ' * dum 

dh 

du± 


dim dim 

du\ ’ ’ ’ dum - 


(2.41) 


(2.42) 


When an operating point is also an equilibrium 
point (derivatives are equal zero), then f(x,u) = 0 
and the system reduces to Ax = AAx + BAu where 
future equations omit the A symbol for readability. 


2.5.2 Process 


set of states that are as close as possible to the ideal 
perturbed state vector, and call this new set the ac- 
tual perturbed state vector. Run lsqnonlin opti- 
mization to find the “closest match”, which mini- 
mizes the Euclidean norm of the difference between 
the ideal perturbed state vector and the actual per- 
turbed state vector. 

Using x+ and xj % two optimizations return two 
actual perturbed state vectors x + and xj. These 
state vectors are passed into the system dynamics, 
which return two derivative function vectors 

h = f(xp u), f~ = f(xj , u) (2.43) 

Working through each row of the derivative vector 
yields the partial derivatives 

dfi f ji f ji f ji f ji rn a a\ 

aT 05 wj = < 2 - 44) 

JL j JL j LtO 

The code runs through a loop and performs this pro- 
cess for each state xj for j = 1, . . . , n which popu- 
lates the A matrix. An identical process is used with 
the control inputs u to derive the B matrix. 


For simple systems with closed form solutions, the 
partial derivative terms that form A and B matri- 
ces may be obtained by hand calculations. Because 
the CSR system requires a numerical routine (does 
not have a closed form solution) this is not a viable 
option. Matlab has built-in functions that linearize 
a wide range of systems. However, the Matlab pro- 
cess is expecting a model that resembles Figure 2.23, 
which is not suitable for multibody problems. 

An alternative approach is used. For each entry 
in the A matrix select a single state (xj) at equi- 
librium. Adjust this value by a small increment 
6 in both directions of its nominal trimmed value 
while all other states are kept at their equi- 
librium values. For reference, call this the ideal per- 
turbed state vector. Traditional trimming routines 
run the system dynamics at this point, but for a 
multibody problem the perturbed system resembles 
Figure 2.24, which is not a valid configuration. 

The ideal perturbed state vector does not main- 
tain the physical constraints of the problem, so an 
alternative vector is sought. The goal is to find the 


2.5.3 RTP State Space Form 

This section applies the linearization process to the 
Round the Pole system, where the state Matrix A 
and input matrix B are provided in (2.45). The first 
step defines the states and inputs within the system. 

There are three control inputs for each satellite; 
throttle, elevator, and rudder. These control signals 
influence the relative positions between the hub and 
the satellite. The throttle impacts relative horizontal 
states, the elevator impacts relative vertical states, 
and the rudder impacts relative heading states. 

The RTP system has ten total states. The first 
two are the roll and pitch of the satellite. The next 
three express relative positions between the satellite 
and hub. They are heading, horizontal position, and 
vertical position. The final five states are the deriva- 
tives of the first five states. 

Initially, it is assumed that all states are available 
for measurement, which allows for full state feed- 
back control. However, standard GPS is not accu- 
rate enough to measure the relative positions of the 
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satellites. More elaborate GPS techniques, like dif- will not be available, and introduces controllers with 
ferential GPS, RFID triangulation, landmark track- state estimation and reduced order state estimation 
ing, and INS navigation might offer a localization so- which can control the system without measuring the 
lution after further review. As a worst case scenario, exact location of the satellite, 
this research assumes that relative satellite positions 



Chapter 3 

Controls Theory 


The system dynamics are implemented within the simulation environment, and the vehicle 
model is trimmed and linearized. The current task reviews the controllers developed for the CSR 
vehicle. Fundamental concepts are presented which describe the system to be controlled and 
the properties used during control design. Controllers with progressive levels of sophistication 
are introduced and related back to their purpose for the CSR system. 


Motivation 

The CSR concept vehicle is a radical departure from 
customary aircraft design, so there is no existing re- 
search to indicate that a particular control strategy 
is best suited for this type of system. The general 
design methodology for this research starts with ba- 
sic models and controllers, learns more about the 
particular system, and then refines and expands the 
models and controller capabilities over time. 

The full CSR system does not have any direct 
control over the central hub. Stabilizing and ma- 
neuvering the hub, as well as the rest of the system, 
is accomplished indirectly through the tension in the 
tethers, which is manipulated by the control inputs 
applied to each individual satellite vehicle. Inves- 
tigating these interactions is the first step towards 
developing the system controllers. 

Understanding the background theory used dur- 
ing the controller development is of utmost impor- 
tance in understanding how the controller is imple- 
mented. The purpose of this chapter introduces con- 
cepts of state feedback and state estimation to a 
technical audience not well versed in control design. 
Entire text books address the underlying theory of 


these controllers, so this content strives to deliver 
the most pertinent concepts needed to understand 
basic principles used for the CSR controller. 

Objectives 

The primary goal of this initial research is to demon- 
strate that the CSR vehicle is a viable design con- 
cept, and that a controller exists which can success- 
fully navigate the system. Because no control theory 
is established for this system, the research on the 
CSR concept vehicle includes several different con- 
troller methodologies. Prerequisite knowledge is re- 
quired before applying these techniques, so the goal 
is to present fundamental concepts and then gradu- 
ate to more advanced strategies. 

Before introducing the controllers, it is impor- 
tant to understand the concept of a plant, which is 
the system to be controlled. Describing how to ex- 
press a system in state space form, how to determine 
the eigenvalues, and understanding how those eigen- 
values contribute to the overall response of the sys- 
tem, is required to understand the controller theory. 
Other topics are directly related to the controller 
development. Theory on controllability and observ- 
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ability, and canonical transformations are needed to 
understand the controller formulation. 

After developing a firm foundation, the various 
controllers are presented. The analysis works from 
simple techniques to more advanced concepts, so the 
starting point for the controls analysis begins with 
full state feedback applied to single-input single- 
output (SISO) subsystems, which provides a foun- 
dation for multiple-input multiple-output (MIMO) 
control. Full state feedback assumes that all the 
satellite states are available for measurement, which 
is not possible due to the limitations with GPS. State 
estimation addresses this issue, where a reduced 
order estimator is required to estimate unknown 
states. Finally, parameter identification and miti- 
gating model uncertainty is only achieved through 
adaptive control techniques. 

Outline 

Before assessing any eternal flight capabilities, The 
CSR vehicle must demonstrate an ability to be con- 
trolled during flight. This chapter focuses on the 
controls theory and analysis to address this issue. 

Section 3.1 presents the Controller Methodol- 
ogy which makes the CSR control problem more 
tractable. It describes the available control modes, 
presents an alternative perspective, and defines the 
roles of the inner and outer loop feedback systems. 

Section 3.2 outlines the System Analysis per- 
formed on the CSR model. It describes its state 
space representation, transfer matrix, eigenvalues, 
and system poles and zeros. Topics pertaining to 
controllability, observability, and decentralized con- 
trol are related back to the system model. 

Section 3.3 describes SISO State Feedback which 
is the first control implemented on the vehicle. A 
custom designed system response is put into a canon- 
ical transformation, and implemented through the 
state feedback controller. 

Section 3.4 expands to MIMO State Feedback 
which accommodates coupling between input and 
output channels. The SISO structure is general- 
ized to handle multiple inputs and multiple outputs, 
while still achieving a custom closed-loop response. 

Section 3.5 addresses State Estimation which is 
used when some of the CSR states are not directly 


measured. Control inputs and measured outputs are 
used to estimate the states and reduce the error to 
zero over time. Then a reduced order estimator is 
implemented to only estimate unknown states. 

Section 3.6 introduces Adaptive Control which 
accomplishes system identification and mitigates 
modeling uncertainties. It presents a parameter 
identification technique where the adaptive laws en- 
sure stability and parameter convergence. Finally, it 
discusses model reference adaptive control (MRAC) 
and L\ adaptive control strategies. 

3.1 Controller Methodology 

At first glance, controlling the CSR vehicle appears 
very complicated, which requires full cooperation 
from all entities involved. However, this research 
will show that a slight shift in perspective makes the 
control problem much more tractable. 

This section outlines the different control modes 
available for the CSR system, which leads into an al- 
ternative perspective describing the big picture ap- 
proach developing the system controller. It describes 
the roles of the inner and outer loop controllers, and 
draws a parallel between the controller development 
and the prototype hardware development. 

3.1.1 Control Modes 

Before embarking on the specifics of the CSR con- 
troller development; first recognize that the novel 
CSR design operates in a perpetual state of rota- 
tion. Thus, it is important to redefine some common 
conventions used within general aviation. 

Roll and pitch are synonymous, because there is 
no “nose” or “side” to differentiate between them. 
Roll and pitch describe the same flight maneuver for 
any given direction in the inertial frame. Likewise, 
yaw is arbitrary, because the system rotates con- 
stantly. Yaw only has meaning referring to an instru- 
ment (like a camera system) that maintains a static 
heading. Accounting for the specifics of the CSR 
geometry, only two inertial descriptions are needed: 
vertical translation and horizontal translation. 

Borrowing from rotorcraft terminology, vertical 
translation utilizes collective control inputs, and hor- 
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izontal translation uses cyclic control inputs. Hori- 
zontal and vertical translation both have two control 
modes available to induce each type of motion. 

For vertical translation, the first mode applies 
collective throttle commands. This changes the 
satellite velocity which alters the lift force. Chang- 
ing the lift force in unison results in a climb or a 
descent. The other mode applies collective elevator 
commands. This changes the satellite pitch angle, 
which points the nose of the satellite up or down. 
With all satellites pointing upward or downward, the 
system enters into a climb or descent. 

Similarly, two cyclic control modes exist for hor- 
izontal translation. The first mode applies cyclic 
rudder commands. This changes the relative head- 
ing of the satellites, which adjusts the radial tension 
in the tethers acting on the hub. The command 
is cyclic, so one point in the rotation has maximum 
outward tension, and the opposite side has minimum 
radial tension. Because the forces are imbalanced, 
the CSR system translates horizontally in the direc- 
tion of maximum tension. Cyclic elevator commands 
can also achieve horizontal translation. This com- 
mand influences the relative vertical position of the 
satellites, so each satellite flies through a low point 
on one side and a high point directly opposite. This 
coordinated flight has the effect of tilting the rotor 
plane, which tilts the overall force vector and induces 
horizontal translation. 

3.1.2 Alternative Perspective 

Controlling the CSR vehicle appears to be a very 
complicated task. With many different states, mul- 
tiple tethers and satellites, and force interactions 
pulling on the central hub, the problem appears to 
require full cooperation from all entities involved. 

In many ways, the coordinated effort resembles 
a highway traffic system. When viewed from above, 
the number of vehicles and potential routes make the 
traffic system appear unmanageable. However, driv- 
ing is a very simple task as as seen from the driver’s 
seat. As long as each driver adheres to basic rules 
(stay centered within a lane, keep a safe following 
distance, obey traffic signals), the complex driving 
system becomes much more tractable. 

The same approach can be applied to the CSR 


vehicle. Imagine sitting on the hub while peering at 
an individual satellite, operating in a trimmed state 
with no external disturbances. From this perspective 
the satellite will maintain a static position relative 
to the hub. It will appear to be “fixed” in space. 

The control task identifies meaningful relative 
positions of the satellite with respect to the hub, 
then applies control signals to manipulate the rel- 
ative states of the satellite. Altering the relative 
states changes the tether force vector anchored at 
the hub. Manipulating the tether anchor forces con- 
trols the central hub and the entire CSR system. 

3.1.3 Inner and Outer Loop Control 

Previous sections described the motion of the CSR 
system in the inertial frame, and how to control 
satellites with respect to their local frame. These 
two tasks indicate that two control loops are required 
for the control system. An inner loop controller is 
responsible for stabilization, and an outer loop con- 
troller executes the waypoint navigation. 

Inner loop control works locally on each satel- 
lite. The controller adjusts the throttle, elevator, 
and rudder of each satellite. These inputs maintain 
the desired relative positions with respect to the cen- 
tral hub. This subsystem is mirrored as a Round the 
Pole (RTP) flight model. 

Outer loop control works globally on the central 
hub. This moves the CSR vertically and horizon- 
tally within the inertial frame. The controller ac- 
complishes this task by manipulating tether forces 
acting on the hub anchor points. It achieves this 
because tether anchor forces are mirrored by the rel- 
ative positions of the satellites. 

3.1.4 Parallel with Hardware Testing 

The CSR project follows a spiral development life cy- 
cle. Because of the novel design, there is very little 
existing research to assist the development process. 
To minimize risk during early development, small in- 
cremental steps are taken to learn about the system, 
collect data, and refine simulation models. More ad- 
vanced techniques are introduced over time, building 
upon previous lessons learned. 
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The first prototype system implements RTP 
flight. For this model, a single satellite flies around a 
fixed pole. The main task identifies trimmed control 
inputs and performs system identification on the ve- 
hicle. Collected data is compared to simulation data, 
and the computer model is revised to more closely 
match the physical system. 

The next milestone is a fixed flight. This model 
has three satellites constrained to a pole. The pri- 
mary objective maintains equal spacing between the 
satellites, and tests coordinated flight maneuvers. 
Equal spacing between satellites is a requirement be- 
fore graduating to a free flight test. Similarly, coordi- 
nated flights demonstrate collective and cyclic com- 
mands executed by the satellites. These capabilities 
are accomplished within a controlled environment, 
before advancing to a free flight. 

The final stage implements the free flight model. 
This is the fully unconstrained CSR system. The 
objective is to demonstrate hover abilities and to 
achieve waypoint navigation through coordinated 
flight maneuvers. Once a flying prototype is avail- 
able, the research will evaluate disturbance rejection 
and the effectiveness of various adaptive controllers. 

3.2 System Analysis 

Previous work described the modeling and dynamics 
that describe the RTP and CSR systems. This sec- 
tion presents concepts that analyze the plant from 
a controls perspective. It introduces the idea of a 
plant, and describes the state space and transfer 
function forms. It covers the characteristic polyno- 
mial of the system, the process for determining the 
system eigenvalues, and how those relate to the poles 
and zeros of the system. Finally, it addresses con- 
trollability, observability, and decentralized control. 

3.2.1 System Description 

The plant is the physical system to be controlled. 
It is a mathematical model that relates the control 
inputs u(t) to the system outputs y(t). A linear 
time- invariant (LTI) model of a physical system is 
expressed in state space form as 

x(t) — Ax(t ) + Bu(t ), y(t) = Cx(t ) (3.1) 


where x(t) G is the state vector, u(t) G M m is 
the control input, y(t) G R q is the system output, 
A G R nxn is the state matrix, B G R nXm is the in- 
put matrix, and C G R qXn is the output matrix. The 
most general form of a state space system includes a 
feedforward matrix, D G R qXm . However, this term 
does not appear in physical systems, and is dropped 
from subsequent derivations. 

The system (3.1) describes a multiple-input 
multiple-output (MIMO) system. For single-input 
single-output (SISO) systems, the input u(t) and 
output y(t) become scalar values, and B and C ma- 
trices reduce to vectors with appropriate dimensions. 

3.2.2 Transfer Function/Matrix 

The transfer function G(s) of a SISO LTI state space 
model expresses the relationship between the control 
input and the system output of the system with zero 
initial condition and zero point equilibrium. The 
plant transfer function G(s) can be obtained directly 
from the state space model. 

Convert the time domain representation of (3.1) 
to the frequency domain with the Laplace transform 

sX(s) = AX(s) + BU(s), Y(s) = CX(s) (3.2) 
solve for X(s) as 

sX(s)-AX(s) = BU(s) 

(si — A)X(s) = BU(s ) (3.3) 

X(s) = (si - A)~ l BU(s) 

substitute X(s) into the output equation 

Y(s) = C(sl - A)~ 1 BU(s) (3.4) 

and define the transfer function G(s) as 

G(s) = C(sl — A)~ l B (3.5) 

which satisfies the relationship 

°< s > = W) (3 ' 6) 

as desired. This relationship is expressed in block 
diagram notation depicted in Figure 3.1. 
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Figure 3.1: Block diagram of the plant model. 
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Figure 3.2: System poles in the complex plane. 


For a SISO system, the transfer function G(s) 
is single a rational polynomial that relates a single 
input to a single output. For MIMO systems, G(s) 
becomes a transfer matrix which is populated with 
rational polynomials within the matrix elements. 


3.2.3 System Eigenvalues 


The transfer function G(s) has an equivalent repre- 
sentation formed by replacing the inverse term 

C[Adj(sI - A)]B = n(s) 

det (si -A) d(s)' y ’ 

The denominator d(s) = det(sI-A), called the char- 
acteristic polynomial, describes the response charac- 
teristics of the system. It is calculated as 


si- A I 
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and solving for d(s) = 0 yields the system eigenval- 
ues, A i for i — 1,2, ..., n. 

Eigenvalues (also called system poles) have real 
and imaginary parts, and their location in the com- 
plex plane dictates the response characteristics of the 
system. The real part of a pole, Re[A^], determines 
its horizontal position. Poles in the Left Hand Plane 
(LHP) are stable, poles on the vertical axis produce 
a bounded oscillation, and poles in the Right Hand 
Plane (RHP) are unstable. The imaginary part of 
a pole, Im[Ai], form complex conjugate pairs above 
and below the horizontal axis. The greater the mag- 
nitude away from the axis, the greater the frequency 
of oscillation. Figure 3.2 illustrates poles and their 
respective system response. 

3.2.4 Poles and Zeros 

The numerator and denominator polynomials of a 
SISO transfer function determine the poles and ze- 
ros of the system. The roots of the numerator n(s) 
are system zeroes, and the roots of the denominator 
d(s) (or characteristic equation) are poles of the sys- 
tem. Poles of SISO systems are also the eigenvalues 
of the state matrix A. 

Transfer matrix poles and zeros carry slightly dif- 
ferent meaning from their transfer function coun- 
terparts. For nearly all s, a MIMO transfer ma- 
trix G(s) = N(s)/d(s ), has rank[G(s)] = min {m,g}. 
System zeros are defined such that N(s) drops rank, 
and system poles are d(s) = det (si — A) = 0. 

The Smith-McMillan Form transforms the trans- 
fer matrix G(s) into a useful form for evaluating 
poles and zeros. Given e^(s) and Si(s) are a pair 
of monic and coprime polynomials, e^(s) is a factor 
of ei+i(s), and Si(s) is a factor of 5i-i(s), the Smith- 
McMillan form of G(s) is described by 

Gsm(s) = diag ( Yj-\, • • • , jj~li 0, • • • , o) (3.9) 
l<h(s) d r {s) ) 

where r — rank (N(s)). 

3.2.5 Controllability 

Controllability describes whether the control inputs 
can influence the states of the system. It is an im- 
portant tool for system transformations and imple- 
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meriting MIMO state feedback. The controllability 
matrix C is given by 

C = [B, AB , A 2 B , • • • , A n ~ 1 B] e R nxnm (3.10) 

where the system (A, B ) is controllable, if C is full 
rank (rank(C) = n), which signifies that all rows and 
columns are linearly independent. 

Controllability indexes perform a bookkeeping 
function while working with MIMO systems. Define 
column vectors of B as 

B = [ bi, b 2 , ■ ■ ■ , b m ], h e M n (3.11) 


Significant coupling between inputs and outputs 
means SISO techniques are not appropriate. How- 
ever, when there is strong correlation between in- 
put/output pairs, and little cross coupling, then de- 
centralized control is a valid technique. 

The objective is to identify one-to-one mappings 
between control inputs and system outputs, with 
little coupling between other inputs and outputs. 
Strong correlations permit SISO design techniques. 
For example, for the system Y(s) = G(s)U(s) 


’ yi(s) ‘ 


9u{s) 

912(a) 

ui(s) 

. 2/2 (s) _ 


_ 921 (s) 

922(s) 

_ u 2 {s ) 


then, select the first n linearly independent columns 
of C, and reorder them as 

[6i, Ab i,---, A^~%, • • • , 6 m , • • • , A^bm] (3.12) 

where /i^, i — 1, . . . , m are controllability indexes. 


decentralized control is appropriate when < 712 ( 5 ) and 
c /21 ( 5 ) are “small”, because it indicates that 1 / 1 ( 5 ) 
maps to 1/1(5) and 1/2(5) maps to 1/2(5). The Rela- 
tive Gain Array (RGA) is a matrix A that quantifies 
the subjective term “small”. It is described by 

f An • • • ^lm I 


3.2.6 Observability 

Observability describes whether initial states can be 
observed from the system outputs. It is an important 
tool for state estimation when states are not directly 
measured. The observability matrix O is given by 

0= [ 1 C,CA,CA V- - .CT-fer- (3.13) 


A = 



mm 


(3.17) 


where A ij = [G( 0)]^- [G _1 (0)]^. Large positive 
terms on the matrix diagonal indicate strong corre- 
lations. Nearly zero values for off-diagonal elements 
imply little coupling between I/O pairs. 


where the system (A, C ) is observable, if O is full 
rank (rank(O) = n), which signifies that all rows 
and columns are linearly independent. 

Observability indexes serve the same function as 
controllability indexes. Define row vectors of C as 

c =[c 1, c 2 , • • • , Cq, ] T , cf e M n (3.14) 

then, select the first n linearly independent rows of 
O, and reorder as 

[ci, Cl A, • • • , ClA 1 ' 1-1 , • • • , Cq, ■ ■ ■ , CgA Uq ~ 1 ] T (3.15) 
where i = 1 , . . . , q are observability indexes. 

3.2.7 Decentralized Control 

Developing MIMO control theory is more involved 
than SISO control theory, so applying SISO de- 
sign techniques to a MIMO system is attractive. 


3.3 SISO State Feedback 

SISO systems have one control input ?/(£), any num- 
ber of states x(t), and one system output y(t). These 
systems are less complex than MIMO systems, so un- 
derstanding the behavior of the SISO subsystems is a 
natural progression toward a more complex analysis. 

State feedback control is a well established con- 
trol technique, straight forward to implement, and 
allows for arbitrary selection of the desired system 
response. Although, state feedback only utilizes a 
state gain K and a reference gain this control tech- 
nique achieves any desired system response (within 
the physical limits of the control inputs). The sim- 
plicity of the controller, and the flexibility to design 
a prescribed system response, makes this controller 
an ideal starting point. This section introduces fun- 
damental concepts needed to understand the theory 
implementing state feedback control. 
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3.3.1 System Response 

Consider a generic system represented by a second 
order dynamic model, described by 


which forms an alternative state space realization 
with transformed matrices of the form 

r o i oi 


G( s ) 




s 2 + 2(uj n s + (J 2 


(3.18) 


where ( and uo n are the desired damping ratio and 
natural frequency, respectively. The characteristic 
polynomial is described by 


d(s) — (s + a + ooj)(s + a — c oj) (3.19) 
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Tl n —1 


(3.25) 


where Re [A] — cr, and Im[A] — c u. 

When designing a feedback controller, it is more 
convenient to specify the desired dynamics in terms 
of settling time t s and percent overshoot M p . How- 
ever, implementing the controller requires knowledge 
of the system poles. Understanding the relationships 
between the system dynamics and the pole locations 
aids the design process. 

The settling time t s dictates cr, the percent over- 
shoot M p yields £, and a and ( are used to determine 
uo and uo ni from the following relationships 

cr 4.6 jt s 

hmYE 

7 r 2 + In {Mp) 2 
LO n = cr/C 
U = U n \Jl ~ C 2 . 

Starting with the desired system characteristics 
(£ s ,Mp,£), the previous equations are used to find 
the system poles (cr, cj), needed for control design. 

3.3.2 Canonical Transformation 

Before developing the state feedback controller, the 
system must be transformed into a canonical form. 
This alternative description is equivalent to the orig- 
inal system, but has a special simplified structure 
used during the controls development. 

The plant transfer function G(s) is expressed as 
rational polynomials of n(s) and d(s), given by 

g (.)=^= n ^ r 1+ ;" +ni ; i+n ° (3.24) 

d(s) s n +d n -is n 1 H \-diSi +d 0 



(3.20) 

(3.21) 

(3.22) 

(3.23) 


The control input can manipulate every state within 
the model, so the special structure of this realization 
guarantees the system (H, B ) is controllable. An im- 
portant point to note, is that similarity transforms 
do not alter the original system eigenvalues. 

The mapping between the controller canonical 
form and the original representation, is handled by 
a transformation matrix T G R nXn , where 

x = Tx , A = TAT~ 1 , B = TB, C = CT~ 1 (3.26) 

and the the transformed system is described by 

x{t) — Ax(t ) + Bu[t ), y(t) — Cx(t). (3.27) 

Transformation matrix T is obtained with the con- 
trollability matrix C. Starting with (3.10), compute 
the last row of the controllability matrix inverse as 

t n = [ 0 0 • • • 1 ]C _1 (3.28) 


construct the transformation matrix inverse as 


T -i 


t n A 
t n A n ~ 2 

tnA n ~ l 


(3.29) 


and use the transformation matrix T to obtain the 
transformed state space realization (3.27). The 
derivation began with a fully controllable realization, 
then developed the transformation. Therefore, when 
starting with the transformation to find the canoni- 
cal form, the system (A, B ) must be controllable. 
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and the closed-loop characteristic polynomial is 


d(s) = s n + (d n - 1 - k n -i)s n 1 . , 

+ • • • + ( d\ — k\)s + (g?o — fco). 


Figure 3.3: State feedback block diagram. 

3.3.3 State Feedback Development 

The goal of state feedback is to generate a control 
input u(t) using direct measurements of the states 
x(t) to arbitrarily place the poles of the closed-loop 
system in prescribed locations. With a SISO LTI 
plant described by (3.1), define control law u(t) as 

u(t) = Kx(t ) + kr{t) (3.30) 


This simplified structure allows for direct matching 
between coefficients in the closed-loop characteris- 
tic polynomial d(s), and coefficients in the desired 
characteristic polynomial a(s) 

cy(s) = s n T OL n — \s n -F • • • T ot\s T cro* (3.35) 

Select the transformed feedback gain K as 
K—[ do — o^o, d\ — ot\, ••*, d n -\ — a n -i ] (3.36) 


with r(t) G M as a reference input, constant state 
gain vector K T G R n , and constant reference gain 
k G R. The block diagram structure of the feedback 
controller is displayed in Figure 3.3. 

Substituting the control law (3.30) into the plant 
(3.1), produces the following closed- loop system 

x{t) = (A + BK) x(t) + Bkr(t) , , 

= Ax(t) + Br(t) [ } 

where 

A = A + BK, B = Bk , (3.32) 

describes the desired closed-loop system response. 
The design task is to find K such that the eigen- 
values of the closed-loop system A i(A + BK) match 
the eigenvalues of some desired characteristic poly- 
nomial a(s), and to find k so the output y(t) tracks a 
constant reference set point r(t) = r. The state gain 
K determines transient response, the reference gain 
k dictates the steady state response, and the resul- 
tant system follows the matching condition (3.32). 

When the system is expressed in controller 
canonical form (A,B,C), the closed-loop state ma- 
trix A is described by 


which places the closed-loop poles of A in the desired 
locations. Because K is designed for the canonical 
form, it must be transformed with K — KT to work 
for the original system. 

The final step determines the reference gain k 
which achieves steady-state set point tracking. The 
steady-state of a system is described by the equilib- 
rium point once the state derivatives settle to zero. 
The closed-loop state equation is expressed as 

0 = Ax e + Br , , 

x e = —A~ l Br 

and the output equation becomes 

y e Cx e — —CA~ l Br. (3.38) 

For y e to track r, we need —CA~ X B — 1 , which is 
equivalent to 

-C(A + BK)~ 1 Bk = 1. (3.39) 

Solving for the reference gain k yields 


A = A + BK 

0 1 0 

0 0 1 

_ ~ do+fco —d\+ki ■ ■ ■ —d n -i+k n -\ 


(3.33) 


k = —\C(A + BK)~ l B]~ l 
= -[ CM " 1 #]- 1 


(3.40) 


which forces the steady state system output y e to 
track the reference signal r. 
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3.4 MIMO State Feedback 

MIMO systems have multiple control inputs u(t) and 
multiple system outputs y{t). The theory for MIMO 
controllers is more involved than SISO controllers, 
but MIMO control accounts for coupling between 
inputs and outputs, which offers a more unified re- 
sponse. Furthermore, most adaptive control and sys- 
tem identification theory is based on MIMO systems, 
so this controller is the next logical progression. 

State feedback for MIMO systems maintains the 
same structure as SISO state feedback. It still only 
utilizes two gains, and achieves any desired system 
response. The two gains become matrices K x and 
iF r , accounting for additional inputs and outputs. 

Although the control structure is identical, ob- 
taining the gain matrices requires some additional 
derivations. This section presents the theory for de- 
veloping MIMO state feedback, selecting a desired 
system response, and obtaining the gain matrices. 


locations, and to have the output signal y(t) track a 
constant reference set point, r(t) — r. To implement 
this control architecture, the system (A, B ) must be 
controllable and rank(F>) = m. 

3.4.2 Canonical Transformations 

Identical to the SISO development, this process be- 
gins by putting the system into a canonical form 
through a transformation T G W nXn . The trans- 
formed arrays have the same form as (3.26) and the 
transformed system is still described by (3.27). 

With (A, B ) controllable, a transformation T ex- 
ists for the MIMO system, such that 

A = A 0 + B 0 UT~\ B = B 0 R (3.44) 

with 

Aq = diag{ Aqi , Aq 2 i • • • , Ao m } , > 

B 0 =diag{£oi, B 0 2 , ••• , B 0m } ^ °> 


3.4.1 Theory Overview 


where Aqi is given by 


Consider a MIMO LTI system described by (3.1), 
where system matrices are known and constant, and 
all states are available for measurement. 

Define the control law as 

u(t) — K x x{t) + K r r(t) (3-41) 


' 0 1 0 ••• O' 

0 0 1 0 


Am = 


0 


G R^* x ^* 


1 

0 0 


(3.46) 


with control input u(t) G R m , reference signal r(t) G 
R m , state gain matrix K x G R mxn and reference gain 
matrix K r G R mxm . Substitute (3.41) into (3.1) to 
obtain the closed-loop system 


and Bm is described by 

B 0i = [ 0 0 • • • 0 1 ] T € (3.47) 

for i = 1, . . . , m, where hi are controllability indexes. 


x(t) = (A + BK x )x(t ) + BK r r(t) 
— Ax(t ) + Br(t) 


(3.42) 343 Matrix Formulation 


where A G R nXn and B G R nxm are the reference 
state and input matrices which achieve the desired 
closed-loop system response. These relationships 
form the following matching conditions 

A = A + BK x , B = BK r (3.43) 

which mirrors the SISO relationships in (3.32). 

The design task is to place the closed-loop eigen- 
values, A i(A) = A i(A+BK x ) = Si at some prescribed 


Parallel to the SISO derivations, the MIMO theory 
utilizes a modified controllability matrix as an inter- 
mediate step toward finding the gain matrices. Form 
matrix M -1 as 

[bi,Abi, • • • , (3.48) 

and define row vectors as 

M = [M 1 , M 2 , ••• , M n ] T eK^, 

T = [T U T 2 , ■■■ , T m ] T ER nxn . 1 ’ 
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Matrix M determines the sub matrices of T with 


so the closed-loop eigenvalues are given by 


-VnH 1 “M* 

M^+.-.+^A 


G x n 


L 




and populates matrices f/ and f? by 




e R 


mXn 


L J 


(3.50) 


(3.51) 


' M^A^~ l B 
M^A^B 


€ R mxm . (3.52) 




This yields all the matrices needed to put the system 
into the canonical form. 


A i(A) = + BK X ) = A,(A + M,) (3.57) 

where K x = K X T, which places the system poles in 
the desired locations. 

The final step calculates the reference gain K r , 
so the system output y(t) tracks a constant refer- 
ence set point r(t) — r. Steady-state equilibrium is 
dictated by 


y e = —C {A + BK x )~ l BK r r (3.58) 


so, define the reference gain as 


K r 


C '(A + BKxY 1 B 
: CA~ l B]~ l 


(3.59) 


3.4.4 Closed-Loop Response 


to achieve y(t) = r. 


The closed-loop response is dictated by A where 
d(s) — det (si — A) 

= (s - si) • • • (s - s n ) (3.53) 

= s n + d n —is n 1 + • • • + d\s + o?o 

and the desired system response is described by 
(3.35). The closed-loop state matrix A indirectly 
determines the state feedback gain K x , defined as 

K x = -R~ l UT~ l + R~ l P (3.54) 


where P is given by 


e AU + l 
e AU+M 2 + 1 


G» mxn (3.55) 


m — 1 + 1 

[ — <ao, aq, • • • , —OL n ~ i ] 


with ej being the j th row of an n x n identity matrix. 
This canonical form leads to 


A — A^ + BqP 

= 4) + BAUT- 1 + RK X ) 


3.5 State Estimation 

Up to this point, state feedback theory assumes all 
states x(t) are available for measurement. In prac- 
tice, many systems have states that are not directly 
measured, so full state feedback cannot be directly 
implemented. Later sections will show that the CSR 
vehicle is one such system, because the location of 
the satellites is not precisely known. GPS provides 
sufficient accuracy for waypoint navigation, but does 
not have enough fidelity to provide the satellite lo- 
cation with respect to the central hub. 

State estimation is a powerful control technique 
that can overcome these short comings. It is de- 
signed to make use of all the available measurements 
and then estimate states that are not measured. The 
theory is a direct extension of state feedback, so it 
builds heavily upon the previous foundation. This 
section introduces the general theory for full state es- 
timation, shows the error reduction process, and ex- 
pands the theory to a reduced order estimator which 
only estimates the unknown states. 
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3.5.1 Theory Overview 

State estimation uses information from y(t), u(t ) and 
(A, £>, C) to find the estimate x(t) of x(£), such that 
lim^oo (x(t) — x{t)) — 0, exponentially. The MIMO 
LTI system is still described by (3.1), but now the 
control law becomes 

u{t) — K x x(t ) + K r r(t) (3.60) 

which uses the state estimate x(t) to calculate the 
control input u(t). Substituting (3.60) into (3.1) 
forms the closed-loop system as 

x[t) = Ax[t ) + BK x x(t ) + BK r r(t). (3.61) 

Now, define the estimation system as 

x(t) — Ax(t ) + Bu(t) + K e ( y(t ) — Cx(t)) (3.62) 

with initial condition x(0) = xo, and an estimation 
gain of K e G R nXq . 

The estimation system (3.62) has a very spe- 
cific form. Notice y(t) = Cx{t) is the estimate of 
y(t) = Cx(t). When the estimate x(t) equals the 
true state value x(t), then y(t) — y(t ), the final term 
becomes zero, and K e has no influence. Similarly, 
when the estimate has no error, the other terms in 
(3.62) reduce to (3.1). 

3.5.2 Error Reduction 

The goal is to select K e so the state estimate x(t) 
tracks the actual state x(t), and maintains stability 
within the system. Define the estimation error as 

x{t) — x{t) — x(t) (3.63) 

and solve for its derivative 

x(t) — x(t) — x(t) 

— [Ax+Bu] — [Ax+Bu+K e (y — Cx)] 

— Ax — Ax + K e (y — C(x — x)) (3.64) 

= A(x — x) + K e (y — Cx + Cx) 

= (A + K e C ) x. 

The estimation system dynamics are described 
by A i(A+K e C), so the task is to select K e to achieve 
a desired response (stable, fast, etc). Notice the sim- 
ilarity between state estimation ( A + K e C ) and state 


feedback (A + BK X ). The duality between controlla- 
bility and observability, means selecting K e follows 
the same process as selecting K x . With (A,C) ob- 
servable, it is known that (A T , C T ) controllable. The 
estimation structure is identical to state feedback 
formulation, so the previous process finding the state 
feedback gain (K x from A and B) can be used to find 
the state estimator gain [K e from A T and C T ). 


3.5.3 Reduced Order Estimator 

All feedback controlled systems have some states 
available for measurement, so full order estimators 
have redundancy estimating states that are directly 
measured. A reduced order estimator is more effi- 
cient because it uses the direct measurements and 
only estimates the unknown states. 

Apply a transformation described by (3.26) and 
(3.27) where the transformation matrix T is selected 
to achieve a special form of C, such that 

C=[I q , 0 gx(n _ 9) ]€R« xn (3.65) 

which leads to x = [ x ^ x^ ] T where x m G M. q are 
measured states and x e G MA~ q are estimated states. 
The measured states directly map to the system out- 
puts, so y = x m . The expanded transformed system 
is described by 

x(t) — Ax(t) + Bu(t) 

Xm Amm ^me X m B rn (3.66) 

= . . + „ U 

x e A em A ee x e B e 

where ( A ee , A me ) must be observable. 

The reduced order estimator has a special struc- 
ture of (7, so the previous estimation system (3.62) 
cannot be used. Notice that y — Cx — x m — y\ 
therefore, y — Cx = 0, so an estimator gain has 
no influence on this system. To alleviate this con- 
cern, select a substitution variable w = x e — K 0 x m 
where K 0 G and w G MT~ q . Now the es- 

timated states are determined from the relationship 
x e — w T K 0 x m . 
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The substitution variable dynamics are given by 

w = x e - K 0 x m 

— \A em x m T A ee x e T B e ti\ 

K 0 \A mm x m + A me x e + B m u] 

= [^em!/ T A ee (w + K 0 y) + B e u\ (3.67) 

-K 0 [A 

mm y+A 

me (' w+K a y)+B m'A 
(A ee T(-R e -^-o-^m) ^ 

~b ((^ee K 0 A m e) K 0 -\~A ern K 0 A rnrn ) y 

with u>(0) = wq. Define the estimation system as 


— ^em^m A ee X e ~\~ B e U 
+K 0 [ %m ^A mm X m 4" ^me^e "4 S m , u)] 


(3.68) 


where 4 e is the current value of the estimated state 
x e . Just like (3.62), this estimation system has a very 
special structure. Notice that when the estimate x e 
equals the true value x e , the last term drops to zero 
and K 0 has no influence on the system. Similarly, 
when the estimate has zero error, the other terms in 
(3.68) become the bottom row of (3.66). 

The estimation error is described by 


flight data, meaning the prototypes must fly based 
on imperfect models. 

Adaptive control serves two main purposes for 
the CSR system. First, certain adaptive control ar- 
chitectures perform system identification, which uses 
measured signals to determine system characteris- 
tics. Values like polynomial coefficients and elements 
of state space matrices, can validate computer mod- 
els and aid with future prototype designs. Second, 
adaptive controllers are deigned to adjust controller 
parameters in real time, so the controller starts with 
an imperfect model, and adapts to the true physical 
plant. Early flights, with little aerodynamic valida- 
tion, will not have high-fidelity models. Adaptive 
control mitigates some of this risk by adapting to 
the imperfect model. 

Running system identification on the first set of 
prototypes is the highest priority, so this section fo- 
cuses on adaptive control which performs parameter 
estimation. It also presents a stability analysis, and 
the conditions that must be present for the parame- 
ters to converge to their true values. 


x e (t) — x e (t) - x e (t) (3.69) 3.6.1 Parameter Estimation 


and the estimation error dynamics are obtained by 
taking the derivative of (3.69), and substituting in 
both rows of (3.66) and (3.68). This yields 

x e (t) = x e {t ) - x e (t ) /„ 7n x 

= (A ee - K 0 A me ) x e (t) 

which also has the same structure as the full state 
feedback and full state estimator designs. When 
(A ee , A me ) observable, finding the value of K 0 from 
A^ e and A^ e follows the same procedure as selecting 
K x given (A,B), or K e given ( A T ,C T ). 

3.6 Adaptive Control 

Thus far, all the controllers assume perfect knowl- 
edge of the plant. In reality, the dynamics of the 
CSR system are largely unknown due to its novel de- 
sign. Modeling provides a good starting point, but 
there is very little research and no wind tunnel data 
to compare against the simulations. The primary 
form of model validation will come from collected 


The purpose of parameter estimation is to obtain un- 
known values of the plant. This theory builds upon 
state feedback control, and will yield the values of 
the elements within the A and B matrices. These 
values represent the partial derivatives of the sys- 
tem, so their physical significance is especially useful 
developing prototypes and validating simulations. 

Let A and B be the estimates of A and B , re- 
spectively. Substitute the matrix estimates into the 
matching condition (3.43), solve for A and B as 

A = A- BK~ l K x , B = BK~ X (3.71) 

and solve for the control law (3.41) in terms of the 
reference signal 

r(t) = - K~ l K x x{t ) + K~ l u(t). (3.72) 

The reference system 

x(t) — Ax{t)Br{t) (3.73) 

is identical to (3.42) except the actual states are re- 
placed by reference states x{t). Eliminate the state 
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and reference gains by substituting the reference sig- 
nal (3.72) into reference system (3.73), so that 

1(f) = Ax^+Bl-K-^xit) + K~ l u(t )] 

= Ax\t)_ + (. Ax(t ) - Ax{t)) 

—BK~ l K x x{t) + BR- l u{t) y ’ 
— Ax{t ) + (-A — A)xit) + Bu(t). 

When A — A, B — B, and x — x the system re- 
duces down to (3.1). This suggests the adaptive laws 
should update time varying estimates A(t) and B(t). 

Designing the adaptive observer requires that A 
is stable (all eigenvalues in Re[s] < 0), and the con- 
trol input u(t) and state x(t) are bounded. State 
feedback ensures these conditions are satisfied. The 
closed-loop desired characteristics are reflected in A, 
which is designed to be stable. The boundedness of 
the control input u(t) and the state x(t) were demon- 
strated in the MIMO state feedback development. 


3.6.2 Estimation Errors 

Define the parameter estimate errors as 

A(t) = A(t) - A , B(t) = B{t) - B 

with the following column vectors 

A[t) = [ «i(f), ■ ■ ■ ,a n (t ) ], 
B(t) = [ 61 (f), • • • ,b m (t) ]. 


x(t) — x(t) — x(t) 
with dynamics described by 
x(t) = x(t) — x{t) 

= \Ax[t) +(A(t) —A)x(t) +B(t)u(t)] 
~[Ax(t) + Bu(t)\ 

— Ax(t ) + A(t)x(t) + B(t)u(t). 


(3.78) 


3.6.3 Adaptive Laws 

Select an arbitrary square matrix Q £ w j iere 

Q — Q t > 0. Find P £ R nxn , which is the solution 
to the Lyapunov equation PA + A T P — —Q, where 
P = P T > 0 . Adaptive laws for A(t) and B{t) are 


A(t) = A(t) — —T APx(t)x T (t) 

B(t ) = Bit) = -T B Pxit)u T it) 


(3.80) 


(3.75) 


(3.76) 


(3.77) 


where Ta^b £ W nXn are arbitrary constant gain 
matrices, such that > 0 and > 0 . 

3.6.4 Stability Analysis 

Define a positive definite function V in terms of the 
estimation errors e(£) ((£) is omitted for space) 


V (e) — xA Px + 'Pia T S b 


where 


Sa = = tr 


2=1 


and 


X B = Y,~ b J T B%= tr 

3 = 1 


A t tAA 


B t tAB 


and solve for the derivative of V as 


Because A and B are constant, the parameter esti- 
mate error derivatives are given by A(t) — Ait) and 
Bit) = Bit). The state estimate error is defined as 


V — 2x 1 Px + 2 tr 


A t T- a 1 A 


+ 2 tr 


B t T~ 1 B 


(3.81) 


(3.82) 


(3.83) 


. (3.84) 


Substituting (3.78) into the first term of (3.84) yields 


2 x T Px — 2 x T P 


Ax + Ax + Bu 


— 2 x T PAx + 2 x T PAx + 2 x T PBu 

which yields three new terms. The first term is 


(3.85) 


2 x T PAx — —x t Qx 
the second term yields 


(3.86) 


The estimate errors form the closed-loop state vector 

e(t)= x T , aj, ■■■An, bj, e (3.79) 

which is used for stability analysis. 


2 x T PAx = 2 tr 
= 2 tr 
= 2 tr 
= 2 tr 


x T PAx 
xx t PA 
A T Pxx T 
A t Y- a 1 Y a Pxx t 


(3.87) 
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and the third term follows the same derivation 


3.6.6 Future Work 


2x 1 PBu = 2tr 


B T T B 1 r B Pxu T 


(3.88) 


Substitute all the manipulated terms back into V 
V —x t Qx 


+2tr 

'a t T~ a 1 {T a Pxx t ) 

+ 2tr 

Try(i)l 

+2tr 

B t T~ b \T b Pxu t ) 

+ 2tr 

B T L b 1 (B) 


where the adaptive laws (3.80), cancel all of the trace 
terms leaving 

V = —x t Qx (3.90) 

which shows that V is negative semi-definite. This 
means the equilibrium state e = 0 of the closed-loop 
system is uniformly stable, and the solution eft) is 
uniformly bounded. It also proves that Aft), Bft ), 
and xft) are all bounded. Finally, (3.90) implies that 
xft) G L 2 , so x(t) — 0. 


3.6.5 Parameter Convergence 

The stability analysis shows that all the closed-loop 
signals are bounded; however, Barbalat Lemma does 
not ensure lim^oo Aft) — 0 and lim t^ooBfb) — 0. 
Ensuring the parameter estimates converge to their 
true values requires persistent excitation to make the 
signal “rich” enough. 

Multisine signal injection accomplishes this task. 
The process applies sinusoidal control inputs to the 
system. This excites the system dynamics, so the 
adaptation continues with new measurements. 

The control input is persistently exciting if it con- 
tains at least f — n + m + 1 frequencies. A constant 
control input contains a single frequency at zero. Si- 
nusoidal inputs with frequency uo and amplitude cr, 
contribute two frequencies at ±c o. Therefore, the 
minimum number of sine inputs /o is the smallest 
integer that satisfies /o > (n + m + l)/2. 

Specify the control input as 

fo 

uft) = cq sin (3.91) 

i=i 

where ^ 0, and uoi are all distinct, for i = 
1,2,--- , fo. Once the persistently exciting signal is 
introduced into the system, lim^oo Aft) — 0 and 
lim^oo Bft) = 0, exponentially. 


Collecting data for the first prototype vehicle is the 
highest priority for the hardware development plan. 
Thus, the first adaptive control architecture focused 
on parameter estimation and system identification. 
The next steps with adaptive control will develop 
controllers that adapt to uncertainties within the dy- 
namic model. Two architectures, model reference 
adaptive control (MRAC) and L\ adaptive control, 
will be developed, simulated, and implemented on 
the prototype vehicles. 

Model Reference Adaptive Control (MRAC) 

When it comes to adaptive control, MRAC is the de 
facto standard. It is the most widely researched and 
utilized adaptive control architecture available in the 
literature. An established reference model contains 
the desired system dynamics, then the controller 
feedback gains adapt so the physical system matches 
the desired model. The type of MRAC presented in 
Figure 3.4, yields state estimates. Thus, the CSR 
system can use this architecture even though relative 
positions of the satellites are not measured. Despite 
its widespread use, MRAC has some notable limita- 
tions. High adaptation gains are desirable to accel- 
erate the adaptation process; however, this leads to 
high frequency control inputs which may exceed the 
actuator bandwidth. Selecting appropriate gains is 
something of a black art, where improper selection 
may lead to robustness issues. Finally, MRAC is a 
nonlinear controller, so it is not scalable for alterna- 
tive prototypes and full scale vehicles. 

LI Adaptive Control 

Despite being a relatively new adaptive control ar- 
chitecture, L\ adaptive control has received a great 
deal of attention in recent years, because it alleviates 
many of the issues found with traditional MRAC. 
Similar to MRAC, L\ defines a reference model, 
and the controller adapts so the closed-loop response 
matches the reference model. The key difference, il- 
lustrated in Figure 3.5, is L\ adaptive control intro- 
duces a low pass filter C(s) prior to the control in- 
put u(t). This keeps the actuators within appropri- 
ate bandwidths, which eliminates the high frequency 
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control inputs. Because of this feature, the archi- 
tecture demonstrates guaranteed robustness. There- 
fore, the adaptation gain can be set as high as pos- 
sible, which induces fast estimation, and is only lim- 
ited by the hardware processing speed within the 


controller. Finally, it shifts emphasis away from a 
nonlinear stability analysis, and towards the design 
of a linear filter. This allows for traditional linear de- 
sign analysis, and ensures inputs and outputs scale 
uniformly to specified reference signals. 



Adaptation Law 

Figure 3.4: Block diagram of MRAC control. 



Figure 3.5: Block diagram of L\ control. 









Chapter 4 

System Analysis and Simulations 


Simulations are the primary means of gaining confidence in the system dynamics and controller 
effectiveness prior to flying an untested prototype vehicle. Each round of prototype testing 
exhibits a natural progression from simple tasks to more complex maneuvers, so the simulations 
follow a similar process which mirrors the prototype development effort. 


Motivation 

Very little research exists to aid in the CSR design. 
Whereas, new aircraft designs can utilize existing 
wind tunnel data and control derivatives from other 
similar aircraft, there is no existing model or research 
to start the CSR development. The dynamics are 
highly nonlinear, so paper analysis and hand deriva- 
tions are not sufficient enough for the first proto- 
types. Early controller development will heavily uti- 
lize simulations to gain confidence in the vehicle dy- 
namics before flying a physical prototype. 

The first research focuses on the RTP flight 
model, because this system is used early on for data 
collection, avionics platform development, and sys- 
tem identification. The CSR system does not have 
any direct control over the central hub. Stabiliz- 
ing and maneuvering the hub, as well as the rest of 
the vehicle, is accomplished indirectly through ten- 
sion in the tethers. These forces are manipulated by 
the control inputs applied to each individual satel- 
lite vehicle, which is exhibited in a RTP flight model. 
Therefore, successfully implementing the RTP con- 
troller is the first step towards controlling the full 
CSR vehicle, and this controller will be directly ap- 
plied towards inner loop stabilization. 


Objectives 

With an established set of control theory, the 
next objective is to put that theory into practice. 
The task is to analyze the first prototype vehicle, 
and identify the linear model to develop the con- 
troller. Several types of prototypes will be developed 
throughout the research process, but this first anal- 
ysis focuses exclusively on the control line aircraft, 
because it is the first prototype vehicle. These vehi- 
cles are ready for flight out of the box, which is an 
ideal platform to test the avionics package and flight 
controls before moving on to a much more expensive 
custom aircraft with unknown flight characteristics. 

The CSR control system is composed of inner 
loop control for stabilization, and outer loop con- 
trol for navigation. The first set of simulations ad- 
dress inner loop control which is built upon the RTP 
satellite configuration. The first objective is to as- 
semble the SISO subsystems that corresponds to the 
three control inputs available on the satellite vehi- 
cles: throttle, elevator, and rudder. The structure 
of the system and the desired response is defined for 
each control input, and a control law is developed for 
each input-output pair. Once SISO control is estab- 
lished, the next objective implements MIMO state 
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feedback on the RTP model, which achieves the in- 
ner loop controller used on the CSR vehicle. 

Outline 

While the previous chapter developed the controller 
theory, this chapter applies that theory to the CSR 
system. It uses the analysis techniques on several 
CSR models, and showcases the controller perfor- 
mance by simulating the nonlinear system. 

Section 4.1 performs the System Analysis on the 
CSR vehicle. It describes the model parameters used 
during the simulations, the three control inputs for 
each satellite, the state space signals, and the rela- 
tive gain array for a decentralized control technique. 

Section 4.2 applies SISO State Feedback to the 
round the pole system. It defines the desired system 
response for each input/output channel, and then 
calculates the feedback gains needed to achieve that 
dynamic response. 

Section 4.3 expands the analysis to MIMO State 
Feedback . It reviews the characteristics of the MIMO 
system, reviews controllability and observability, 
and derives the feedback gains for the controller. 

Section 4.4 presents Simulations for both SISO 
and MIMO models. All the controllers are applied to 
the nonlinear CSR model, and several scenarios are 
illustrated which demonstrate that the controllers 
successfully achieve the desired system performance. 

4.1 System Analysis 

Several different prototypes are developed for this re- 
search. The scaled prototype is a carbon fiber Kevlar 
construction, mono-wing design, that more closely 
resembles the full scale CSR system. This prototype 
has many unknown dynamics and aerodynamics, so 
a commercial off the shelf (COTS) control line air- 
plane serves as the first prototype development, be- 
cause it is designed to fly out of the box. 

4.1.1 Model Description 

Before moving to a free flight test, the RTP model 
will serve as a stepping stone. This tests new avion- 
ics hardware and control algorithms before imple- 
menting them in a more difficult environment. This 


section describes the properties of the hub, tether, 
and satellite that generate the simulation model. 

Hub 

For the RTP system the central hub is constrained 
to a fixed pole but is free to rotate about the ver- 
tical axis. The bearing is lubricated, but still pro- 
vides some nominal amount of rotational drag. A 
magnetometer mounted on the hub rotates with the 
system, and provides heading measurements. 

Tether 

The full scale system uses space grade Spectra line 
for the tether, which is extremely light and strong. 
The COTS control line airplane prototype comes 
with a thin steel cable as the tether. The tether 
is 30 feet (9.14 m) long, with a 3/64 inch (1.19 mm) 
diameter, and is modeled with a density of 500 lb/ft 3 
(8,000 kg/m 3 ). The tether is completely flexible, but 
is modeled as 20 rigid body elements. 

Satellite 

The control line airplane resembles a conventional 
aircraft with symmetrical wing and horizontal sta- 
bilizer, and an upper vertical stabilizer. The plane 
uses electric power with LiPo batteries, a brushless 
motor, and a single fixed pitch propeller. The air- 
craft comes with pitch control already installed, and 
yaw control is implemented with a custom addition 
of a servo controlled rudder. The wingspan is 38 
inches (96.5 cm) and the planform wing area is 356 
in 2 (2297 cm 2 ). The vehicle is 24 inches (61.0 cm) 
long with a weight of 2.5 lbf (11.12 N). The wing 
section has an estimated zero angle of attack coeffi- 
cient of lift of q = 0.3, and the wing is inclined 3° 
(0.052 rad) relative to the water line. 

4.1.2 Control Modes 

With a general description of the modeled system, 
the next task formulates how the system will be 
controlled. All control is accomplished through the 
satellite, which has three different control inputs: 
throttle, elevator, and rudder. The next sections 
look at each mode in detail, describe how the control 
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impacts the system, and how this is used with the 
outer control loop. It also describes how the three 
control inputs are mirrored by the three component 
force vectors acting at the hub anchor point. 

Throttle 

The throttle adjusts the amount of force generated 
by the motor system, which influences the horizontal 
position of the satellite relative to the central hub. 
During operation, a lateral tension pulls at the hub 
anchor point. Adjusting the relative horizontal po- 
sition, changes the magnitude of the lateral force. A 
relative position leading the hub increases the lateral 
tension, and a relative position behind the hub de- 
creases the lateral tension. Throttle control is a very 
stable subsystem, because the satellite drag counter- 
acts the throttle force from the motor. 

Elevator 

The elevator adjusts the deflection of the control sur- 
face on the horizontal stabilizer, which influences the 
vertical position of the satellite relative to the cen- 
tral hub. During operation, a vertical tension pulls 
at the hub anchor point. Adjusting the relative ver- 
tical position, changes the magnitude of the vertical 
force. A relative position above the hub increases 
vertical tension, and a relative position below the 
hub decreases the vertical tension. Elevator control 
is a marginally stable subsystem, because the satel- 
lite pivots around a point with no restoring forces. 

Rudder 

The rudder adjusts the deflection of the control sur- 
face on the vertical stabilizer, which influences the 
heading of the satellite relative to the central hub. 
During operation, a radial tension pulls outward on 
the hub anchor point. Adjusting the relative head- 
ing, changes the magnitude of the radial force. A 
relative heading away from the hub increases the ra- 
dial tension, and a relative heading toward the hub 
decreases the radial tension. Rudder control is the 
most stable subsystem, because the centrifugal force 
acting on the satellite wants to remain aligned with 
the tether anchored on the satellite. 


4.1.3 State Space Signals 

To linearize the RTP system, suitable state space 
signals must be established. The following describes 
the inputs, states, and outputs used during the lin- 
earization process. 

Inputs 

There are three control inputs for each satellite, 
where each control input influences a state of the 
satellite relative to the central hub. The throttle 
command 7 , is in Newtons. It adjusts the force gen- 
erated by the motor system, and influences the rela- 
tive horizontal position of the satellite. The elevator 
command a is in radians. It is a measure of the 
elevator angle of deflection, and it influences the rel- 
ative vertical position of the satellite. The rudder 
command /? is also in radians. It describes the rud- 
der angle of deflection, and it influences the relative 
heading of the satellite. Together they form 

u(t) = [ 7 , a, P] T (4.1) 

which is the control input vector signal. 

States 

The RTP system has ten total states. The first two 
states describe roll </> and pitch 6 of the satellite in 
the inertial frame. The next three states express rel- 
ative positions between the satellite and hub. They 
are heading horizontal position 77, and vertical po- 
sition fi. The final five states are the derivatives of 
the first five states. The vector signal 

x(t) = [<p, e, ip, r\ , n, </>, 0, ip, fj, fi, } T (4.2) 

describes the states of the satellite system. 

Outputs 

The outer loop controller uses relative positions of 
the satellites to impart forces on the central hub. 
Because the outer loop specifies relative positions, it 
is natural to design the inner loop controller so that 
the satellites track these relative positions. The out- 
put vector signal is populated as 

y(t) = [ rj, y,, ip ] T (4.3) 
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which are the relative horizontal, vertical, and head- 
ing positions. Because there is a direct mapping be- 
tween states and outputs, define 


C = 


0 0 
0 0 
0 0 


0 1 
0 0 
1 0 


0 0 
1 0 
0 0 


0 0 
0 0 
0 0 


0 0 
0 0 
0 0 


(4.4) 


which is the state space output matrix. 


4.1.5 SISO Subsystem 

The RGA indicates that decentralized control is ap- 
propriate, so the first controller development centers 
around state feedback applied to SISO subsystems 
within the RTP model. The full A and B matri- 
ces are used to extract the linearized SISO models. 
Linearizing the rudder subsystem yields 


Linearized System 

The selection of states and outputs lead directly to- 
ward the output matrix C. A direct relationship 
between the inputs and the states is not analytically 
determined. Following the linearization procedure 
described earlier, yields the state matrix A and in- 
put matrix B. Both matrices are provided in (2.45). 
Together, (A, B,C) form the RTP transfer matrix 
with G(s) = C(sl - A)~ X B. 


A P = 


0 1 

-62.64 0 ’ 


Bp = 


0 

13.78 ’ 


r T - 


(4.8) 


where xp(t) — [ ip(t), ip(t) ] T is based on the relative 
heading. Linearizing the throttle subsystem leads to 


A 7 = 


0 1 

0 2.45 ’ 


B 7 — 


0 

0.54 


, <4= 


i 

0 


(4.9) 


4.1.4 Relative Gain Array 


The RGA uses the DC gain of the transfer matrix 


G(0) — 


and its inverse 


G~\0) = 


-1.26 52.84 0.34 

-0.04 17.93 -0.02 


-0.04 

0.86 

0.30 

-0.90 

2.60 

1.18 

-0.002 

0.06 

0.006 

-0.11 

0.15 

3.45 


(4.5) 


(4.6) 


The elements of the relative gain array A are popu- 
lated with A ij — [G( 0)]^- [G-^O)] which yields 


where x 7 (t) = [ 77 (t), fj(t) ] T uses the relative hori- 
zontal position for feedback. The elevator subsystem 
is marginally stable, so both the relative vertical po- 
sition fi{t) and the pitch angle 6(t) are used as feed- 
back signals. Linearize the elevator subsystem as 


A 


a — 


B 


0 0 10 
0 0 0 1 
0 187.75 9.40 0.02 

0 -59.95 -5.91 0 


0 


' l ' 

0 

, C T = 

0 

13.70 

7 a 

0 

12.41 


0 


(4.10) 


A 


(-l.26K-0.90) (52.84) (-0.002) (0.34)(-0.11) 
(-0.04)(2.60) (17.93) (0.06) (-0.02)(0.15) 
(-0.04)(1.18) (0.86) (0.006) (0.30)(3.45) 

r 1.138 -0.102 -0.036] 

-0.094 1.097 -0.003 

-0.044 0.005 1.039 


(4.7) 


for the CSR system. The matrix has relatively large 
positive elements on the main diagonal, and nearly 
zero values for all other entries. This indicates a 
strong correlation between the input-output pairs, 
so decentralized control is a viable option. 


where x a (t) = [ n(t), 9(t), 9(t) } T . 

The SISO controller requires direct measure- 
ments of the relative states used as feedback dur- 
ing the control loop. Initially, it is assumed that all 
states are available for measurement, which allows 
for this type of control. However, standard GPS is 
not accurate enough to measure the relative posi- 
tions of the satellites. This issue is resolved through 
MIMO state feedback and state estimators. For now, 
SISO simulations serve as a benchmark to compare 
against the MIMO controllers. 
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4.2 SISO State Feedback 

The previous section described the input channels 
as: ( 1 ) throttle for relative horizontal position, ( 2 ) 
elevator for relative vertical position, and (3) rud- 
der for relative heading. It also linearized the full 
nonlinear RTP system and reduced the system into 
three linear subsystems. This first controls analy- 
sis uses those linear models, describes the process 
for selecting the desired system response, and then 
implements state feedback controllers on the three 
SISO subsystems that make up the RTP system. 

4.2.1 Desired System Response 


values are rounded slightly to a = 10 and Ccj = 10. 
This gives the complex pair as 

s + a ± uo j — s + 10 d= 10j (4.15) 

so the desired characteristic polynomial becomes 

a(s) = (s + 10 + 10j) (a + 10 - 10 j) 

= s 2 + 20s + 200 (4 ' lbj 

which yields polynomial coefficients of oto — 200 and 
ol\ — 20. The rudder and throttle both have two 
states, so this polynomial is used directly. The ele- 
vator has four total states, so a second set of complex 
poles are added at the same location, which leads to 


The desired closed-loop system response can be ar- 
bitrarily selected with state feedback. For illustra- 
tive purposes, the first analysis will implement iden- 
tical response characteristics on each of the subsys- 
tems. In future developments, each system may have 
a uniquely optimized dynamic response to accom- 
modate other design considerations, such as power 
consumption or disturbance rejection. 

As a starting point, the desired settling time 
should be around t s = 0.5 seconds, and the percent 
overshoot should be close to M p = 5%. The goal is 
to obtain the real and imaginary parts of the com- 
plex roots, a and a;, respectively. Obtain the real 
part (a = Re [A]) of the complex root with 


4.6 
a — — 

t s 


0.5 


(4.11) 


and solve for the damping ratio as 


C=' 


ln(M p 


7T 2 +ln(M p ) 2 


/ ln(0.05) 2 
' 7r 2 +ln(0.05) 2 


= 0.7. (4.12) 


The natural frequency is determined by 

to n = — = — = 13.14 (4-13) 

so the damped natural frequency becomes 

CJ = UnVi - C 2 = 13.14 \/ 1 - 0.7 2 = 9.39 (4.14) 


a(s) = (s 2 + 20 s + 200 ) (s 2 + 20 s + 200 ) , , 

= s 4 +40s 3 +800s 2 +8000s +40000 ^ ) 

with polynomial coefficients of a q = 40000, oq = 
8000, oi 2 = 800, and 013 = 40. 


4.2.2 Rudder Controller 

Begin with rudder control, because it is the most 
straight forward to implement. Rudder angle /?(£) 
is the input up(t) to the subsystem, and the output 
yp(t) is designated as the relative heading The 
state space realization of the rudder subsystem is 


x/ 3 (t) = Apxp (t) + Bpup (t ) , yp (t) = Cpxp (t) (4.18) 


where Ap, Bp, and Cp are provided in (4.8). Solve 
for the controllability matrix as 


C p = 


0 

13.78 


13.78 

0 


(4.19) 


which determines the transformation matrix 


Tp = 


13.78 0 

0 13.78 


(4.20) 


The canonical form of the subsystem becomes 



0 1 ' 


'O' 

XT 

'13.78' 

-62.64 0_ 

5 Bp = 

1 _ 

> Cp = 

0 


(4.21) 


which yields the imaginary part (uo — Im[A]) of the 
complex root. For simplicity during the control de- 
velopment and simulation, the real and imaginary 


which has the following characteristic polynomial 

dp(s) = s 2 + 0s + 62.64. (4.22) 



CHAPTER 4. SYSTEM ANALYSIS AND SIMULATIONS 


65 


Direct matching yields the transformed state gain 

Kp = [ -137.36 -20 ] (4.23) 

which determines the system state gain 

Kp = KpTp 1 = [ -9.97 -1.45 ] . (4.24) 

The closed-loop system is described by 

Ap = Ap + BpKp 

0 1 1 (4.25) 

-200 -20 

which has the desired closed-loop eigenvalues. The 
reference gain is given by 

kp = -{CpA^Bp}- 1 = 14.51 (4.26) 

which achieves the desired closed-loop response. 

4.2.3 Throttle Controller 

Throttle control requires one additional step from 
the rudder controller development. It uses the rela- 
tive horizontal position of the satellite with respect 
to the hub, so a rotation matrix is applied to obtain 
this value. The throttle y(t) is the input u^(t) to 
the subsystem, and the output y 1 {t) is designated as 
the relative horizontal position 77 (t). The state space 
realization of the throttle subsystem is given by 

%'y(t)=AjXj(t) + y 1 {t) = C 1 x 1 {t) (4.27) 

where Ay, £> 7 , and C 7 are provided in (4.9). Solve 
for the controllability matrix as 


Direct matching yields the transformed state gain 
iTy = [ -200 -22.45 ] (4.32) 

which determines the system state gain 

if 7 = iTyT" 1 = [ -370.37 -41.57 ] . (4.33) 

The closed-loop system is described by 
A 1 — A 1 T 

0 1 1 (4.34) 

-200 -20 

which has the desired closed-loop eigenvalues. The 
reference gain is given by 

fc 7 = -[C 1 A- l B 1 ]- 1 = 370.37 (4.35) 

which achieves the desired closed-loop response. 

4.2.4 Elevator Controller 

The elevator controller is the most complicated sub- 
system to control. This subsystem is marginally sta- 
ble, so it uses two different states in the feedback 
loop. The elevator a(t) is the input u a (t) to the 
subsystem, and the output y a (t) is designated as the 
relative vertical position y(t). The state space real- 
ization of the elevator subsystem is given by 

x a (t) = A a x a {t) + B a u a {t), y a {t) = C a x a (t ) (4.36) 



0 0.54 

0.54 1.323 _ 


(4.28) 


where A a , B a , and C a are provided in (4.10). Solve 
for the controllability matrix as 
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The canonical form of the subsystem becomes 


4.3.1 MIMO System Characteristics 


A 


a — 


B 


a 


0 1 0 0 ~ 

0 0 10 
0 0 0 1 ’ 

0 -546.1 -60.1 9.4 


' 0 ' 


' 3151 ' 

0 

II 

0.248 

0 

13.70 

_ 1 _ 


0 


(4.39) 


which has the following characteristic polynomial 


The desired system response mirrors the values de- 
veloped for the SISO controllers. The desired closed- 
loop poles are placed at s = — 10 ± lOi. With ten 
total states, that pole location will have a multiplic- 
ity of five, so the desired characteristic polynomial 
for the MIMO system becomes 

a(s) = [(s + 10 + 10i)(« + 10 - 10i)] 5 . (4.45) 


d a (s) = s 4 -9.4s 3 +60.1 s 2 +546. Is +0. (4.40) 

Direct matching yields the transformed state gain 


Using a compact form of notation (used by Matlab), 
let a vector represent the polynomial coefficients 


K a = [-40,000, -7453.9, -739.9, -49.4] (4.41) 
which determines the system state gain 

K a = [ -12.7, -32.0, -4.37, 0.84]. (4.42) 


[ a m *** j a h a o ] a n s n +• - • -\-a\s -\-a$ (4.46) 
so the desired system response is expressed as 


The closed-loop system is described by 
+ B a K a 

0 0 10 
0 0 0 1 
-173.9 -251.1 -50.5 11.6 
-157.5 -457.5 -60.2 10.5 


(4.43) 


which has the desired closed-loop eigenvalues. The 
reference gain is given by 

k a = —[C a A~ l B a \~ l = 12.69 (4.44) 


which achieves the desired closed-loop response. 


a(s) 


1.0e2 

5.0e3 

1.6e5 

3.6e6 

5.9e7 

7.2e8 

6.4e9 

4.0el0 

1.6ell 

3.2ell 


(4.47) 


4.3 MIMO State Feedback 


and the characteristic polynomial is given by 


SISO design techniques selected the exact desired 
system response within all the input/output chan- 
nels. The SISO simulations are an important design 
tool, while analyzing a disturbance for a single chan- 
nel. However, implementing the SISO controllers on 
the CSR is not possible, because they require mea- 
surements of the relative states of the satellites. This 
section presents the MIMO controller development, 
which is able to overcome the SISO limitations. It 
presents the MIMO system characteristics, control- 
lability and observability of the system, the desired 
system response, and the gain matrices needed to 
implement MIMO state feedback. 


d(s) 


det(sl— A) 

1 

14 

140 

1077 

5972 

19369 

73427 

40690 

-53512 

-352074 

-110721 


(4.48) 
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The Smith McMillan Form is given as 


4.3.3 MIMO System Matrices 


Gsm(s) = diag{^,^, 5 (s)} 
g(s) = (s—4.97e— 5±4.79i )(s+ 12. 78)(s — 12.79) 

so the poles are 


The procedure outlined earlier was followed, and 

(4.49) the state feedback matrices were obtained. Matrix 
P, K x , and K r are presented in (4.56), (4.57), and 
(4.58), respectively. 


and the zeros are 


s* = 


1 

io 

00 

± 

3. lli 

0.08 

± 

6.78i 

-0.01 

± 

4.79i 

-1.09 

± 

1.55* 

1.51 



-0.33 



4.97e— 5 

± 

4.79i 

-12.78 



12.79 




4.4 Simulations 

(4.50) Simulations have been developed for two of the con- 
trol architectures presented: RTP SISO state feed- 
back, and RTP MIMO state feedback. The simu- 
lations show very promising results. Due to time 
constraints, simulations for the remaining controllers 
are left as future work. 

(4.51) 

4.4.1 SISO Simulations 


4.3.2 Controllability and Observability 

Controllability and observability are formed directly 
from the MIMO state space matrices (4.4) and 
(2.45). Using (3.10), the first ten columns of the 
MIMO controllability matrix C, are presented in 
(4.54). The matrix is full rank because rank(C) = 
10, so the system is fully controllable. All of the first 
10 columns are linearly independent, so the control- 
lability indexes are 

Mi =4, M2 = 3, M3 = 3. (4.52) 

Using (3.13), the first ten rows of the MIMO observ- 
ability matrix O, are presented in (4.55). The matrix 
is full rank because rank(O) = 10, so the system is 
fully observable. All of the first 10 rows are linearly 
independent, so the observability indexes are 

v\ =4, V 2 = 3, z/ 3 = 3. (4.53) 

This is a very important result for the CSR system. 

Because the system is fully controllable, full state 

feedback can be applied to the system. Because 
the system is fully observable, state estimation and 
reduced order estimation controllers can be imple- 
mented. This result is necessary to move forward 
when GPS is not available to measure the relative 
states of the satellites. 


Several simulations were run to demonstrate the con- 
troller performance on the RTP control line model. 
The SISO controllers are applied to the nonlinear 
model in the simulation to demonstrate the system 
response, disturbance rejection, and reference track- 
ing capabilities. The simulations are shown in Fig- 
ures 4.1, 4.2, 4.3, and 4.4, and the trimmed RTP 
model is displayed in Figure 4.5. 

The first three figures present the system re- 
sponse when only the specified controller is applied, 
and the other two control inputs are left at their 
trimmed values. Figure 4.4 simulates the system 
when all control subsystems are applied and distur- 
bances are applied to each input-output pair. 

Figure 4.1 shows the simulation when throttle 
control is applied, and elevator and rudder inputs 
maintain their trimmed values. At 0.5 seconds, the 
system has a 5 Newton impulse disturbance applied 
at the nose of the satellite. At 1.5 seconds, the refer- 
ence signal changes so the system must track a new 
relative horizontal position of -0.015 meters. 

Figure 4.2 shows the simulation when elevator 
control is applied, and throttle and rudder inputs 
maintain their trimmed values. At 0.5 seconds, the 
system has a 3 Newton impulse disturbance applied 
downward at the CG of the satellite. At 1.5 seconds, 
the reference signal changes so the relative vertical 
position is -0.015 meters of the trimmed value. 
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Figure 4.3 shows the simulation when rudder 
control is applied, and throttle and elevator inputs 
maintain their trimmed values. At 0.5 seconds, the 
system has a 5 Newton-meter torque impulse distur- 
bance applied on the vertical axis at the CG of the 
satellite. At 1.5 seconds, the reference signal changes 
and the system must track a new relative heading of 
-1.625 radians. 

Figure 4.4 shows the simulation when all control 
signals are applied, and the system experiences dis- 
turbances on all three subsystems. At 0.5 seconds, 
the system has the disturbance applied to the rudder 
subsystem, at 1.5 seconds, the system has the dis- 
turbance applied to the throttle subsystem, and at 
2.5 seconds, the system has the disturbance applied 
to the elevator subsystem. 

Each simulation shows that the closed-loop re- 
sponse has a settling time around 0.5 seconds, with 
approximately 5% overshoot, as specified during the 
controller design. Additionally, each individual sys- 
tem demonstrates that the overall system has signif- 
icant decoupling between states, so that each con- 
trol input mostly influences the desired system state. 
The final simulation with all control signals applied, 
shows that each input-output pair does not conflict 
with the other input-out pairs, and the full system 
behaves very similarly to the individual simulations. 

4.4.2 MIMO RTP Simulation 

The SISO controllers are easily implemented within 
the simulation, and serve as a good tool for evaluat- 
ing controller performance on a single input /output 
channel. However, the full CSR satellite will require 
MIMO control because that is the form needed to 
perform state estimation when the relative position 
of the satellite is not known. 

To provide an accurate comparison, the full 
MIMO controller was designed with the same re- 
sponse characteristics as the SISO subsystems. The 
system was evaluated, and the gain matrices were 
obtained in a previous section. The simulation of 
the system was subjected to the same disturbances 
as the previous SISO simulations, but the MIMO 
controller was implemented to mitigate the track- 


ing errors. Because the Relative Gain Array (RGA) 
showed very minimal coupling between the other in- 
put/output channels, the MIMO controller simula- 
tion should look very similar to the SISO controller 
simulations. Figure 4.6 demonstrates that this is the 
case, were the MIMO system response is nearly iden- 
tical to the SISO subsystem controllers. 

4.4.3 Future Simulations 

The current simulations show promise for the CSR 
vehicle. The RGA indicates that there are strong 
one-to-one mappings between the inputs and out- 
puts. This result is very favorable because imple- 
menting the control system is a much more straight 
forward process when the system is decentralized. 
The simulations also show that SISO design tech- 
niques can be used to formulate the closed-loop sys- 
tem dynamics, and the MIMO controller will very 
closely mirror the prescribed dynamics. 

The next round of simulations will implement the 
state estimator and reduced order state estimator. 
Existing control theory states that because the sys- 
tem is both controllable and observable, these con- 
trollers will work for the CSR system. However, sim- 
ulations are still needed to demonstrate this type of 
control. Additionally, the theory states that the er- 
ror can be reduced to zero at an arbitrarily fast de- 
cay rate; however, faster dynamics require more ag- 
gressive control authority. Simulating the estimator 
controllers will indicate what decay rate is achievable 
within the actuator bandwidths. 

The final set of simulations prior to hardware 
testing will develop the parameter estimation adap- 
tive controller. Similar to the reduced order con- 
trollers, theory states that a persistently exciting sig- 
nal ensures that the adaptation will push the param- 
eter estimates toward their true values. However, it 
does not give any indication of what signal frequen- 
cies are appropriate. Simulating the parameter es- 
timation is the best way to make sure the multisine 
signal injection is within the bandwidth of the sys- 
tem dynamics, and that the satellite vehicle has an 
acceptable transient response. 
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(4.56) 


3.6e7 5.4e6 3.2e6 -2.5e5 5.9e6 4.0e7 3.3e5 4.3e6 -3.6e5 7.2e6 

K x = | 2.9e0 9.3e0 -8.6e-l -9.2e-3 l.OeO -8.2e0 1.5e-l -8.9e-l 6.0e-3 -7.6e-l 

-2.0e6 — 3.1e5 -1.8e5 1.4e4 -3.4e5 -2.3e6 -1.8e4 -2.4e5 2.0e4 -4.1e5 


(4.57) 


K r = 


1.7e3 — 4.5e4 -2.2e3 

2.6e— 2 — 6.1e— 1 -2.7e-2 
— 9.8el 2.6e3 1.2e2 


(4.58) 
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Figure 4.1: SISO throttle control. 
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Figure 4.2: SISO elevator control, 
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Figure 4.3: SISO rudder control. 
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Figure 4.5: Trimmed RTP system. 
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Chapter 5 


Prototype Testing Plan 


As the saying goes, “The proof is in the pudding.” Theory, analysis, and simulations are all 
valid endeavors and serve as useful design tools. However, a series of plots will never be as 
convincing as a prototype that successfully accomplishes an objective. The only way to truly 
claim the CSR vehicle is capable of eternal flight, is to demonstrate it on a physical platform. 


Motivation 

A simulation is only as good as its model. Even 
with sophisticated analysis, computational dynam- 
ics, and detailed modeling efforts, models always 
contain some discrepancies with the physical system. 
Some systems, like conventional aircraft, have a suf- 
ficient design history where new designs compared 
to historical data provide some level of validation. 
However, radically new concepts, like the CSR vehi- 
cle, have little basis for comparison, so providing a 
physical “proof-of-concept” is of utmost importance. 

Objectives 

The majority of future work focuses on developing 
hardware prototypes that demonstrate the abilities 
of the CSR system. The main development method- 
ology for this research project follows a spiral de- 
sign process. Rough models and prototypes are im- 
plemented as quickly as possible, data is collected 
for those embodiments, the data aids the model im- 
provements, and more refined models lead into more 
sophisticated prototypes. This section outlines the 
plan of action that guides the prototype develop- 
ment process. It introduces the concept of a spiral 
development cycle, describes the prototypes utilized 


along the way, and illustrates the testing performed 
on each prototype. 

Outline 

Simulations serve as useful design tools, but real 
system validation comes from hardware demonstra- 
tions. The only way to claim eternal flight is pos- 
sible, is to implement the vehicle in the real world. 
This chapter focuses on the prototype testing plan 
that progresses the CSR research and design effort. 

Section 5.1 introduces Spiral Development , which 
is an alternative type of design process. Unlike a lin- 
ear design philosophy, spiral development follows a 
cyclic design progression, which advances in small 
and manageable incremental steps. 

Section 5.2 describes the Prototype Vehicles that 
are part of the spiral development process. Each ve- 
hicle serves a unique purpose during each stage of 
development, and progresses from simple models to 
more sophisticated and higher fidelity designs. 

Section 5.3 outlines the Testing Plan that is per- 
formed on each prototype vehicle. With each new 
model, new hurdles and uncertainties are introduced, 
so the testing plan follows a natural progression from 
simple tasks to more risky maneuvers. 


73 



CHAPTER 5. PROTOTYPE TESTING PLAN 

5.1 Spiral Development 

Large engineering projects, like the design of new 
aircraft, traditionally spend a majority of the time, 
effort, and expense on upfront tasks. Spiral devel- 
opment is an alternative development process that 
offers many advantages over a traditional linear de- 
sign path. This section describes the spiral devel- 
opment process, highlights the advantages and dis- 
advantages, draws a parallel to a Design Build Fly 
methodology, and explains how spiral development 
is applied to the CSR research effort. 

5.1.1 Description 

Spiral development is a cyclic design cycle rather 
than a traditional linear progression. Instead of de- 
veloping a single high fidelity model, each cycle uti- 
lizes prototypes to obtain a greater understanding of 
the conceptual system. Early models are very simple 
and only attempt to capture the most dominant sys- 
tem characteristics. Early prototypes are developed 
quickly and inexpensively, which serve as a proof of 
concept. Simulation models are compared against 
data collected from the prototypes. The validation 
process is the crucial step with spiral development, 
where one of two things will happen. 

In the first outcome, the analysis may prove the 
project is infeasible. Unanticipated issues may be 
uncovered, or seemingly small problems are now rec- 
ognized as major concerns. Quite possibly, the con- 
cept design relies on a key component or feature, 
which did not work out favorably in the prototype. 
When this is the case, a considerable amount of time, 
effort, and money is spared rather than proceeding 
forward with a doomed project. 

For the second outcome, the analysis may indi- 
cate the project is feasible and warrants additional 
consideration. Under this scenario, the collected 
data validates or improves the existing computer 
model. The analysis may show that certain simplify- 
ing assumptions need to be refined, or it might iden- 
tify parameter values that were initially estimated. 

This refined model initiates a new round of proto- 
type development. It serves as a design tool for sizing 
features, indicates where revisions are required, and 
reveals potential design improvements. The spiral 
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development process continues in this cyclic fashion, 
where improved prototypes lead to refined models, 
and refined models create improved prototypes. 

5.1.2 Spiral Development Pros and Cons 

Spiral development is not intended to replace a lin- 
ear design process under all circumstances. For some 
projects a traditional design approach is the best 
course of action. However, spiral development is a 
very useful tool that offers notable advantages under 
certain conditions. 

Disadvantages 

Some projects are not well suited for spiral devel- 
opment. Projects that are sufficiently straight for- 
ward may not need a prototype because pure analy- 
sis will suffice. Building a bridge is such an example, 
where analytic tools are available to calculate the 
loads, and there are minimal dynamics to take into 
account. Other times, building a prototype is not 
feasible. A nuclear power plant meets this criteria, 
because a prototype reactor is not practical or safe. 

Advantages 

Despite notable limitations, spiral development of- 
fers some distinct advantages over a linear design 
process. Whether spiral development is appropriate 
reduces down to a simple cost-benefit analysis. A 
linear design requires upfront costs in the form of 
building elaborate models, running monte carlo sim- 
ulations, and processing data. Spiral development 
cuts many of these costs, but includes the time and 
expense building a series of prototypes. To deter- 
mine which design technique is appropriate, simply 
estimate the total design expense for each. For ex- 
ample, systems with complex dynamics, like large 
ship building, may benefit from spiral development. 
An exorbitant amount of effort could go towards 
modeling the hull hydrodynamics, and the result- 
ing model may still include unaccounted errors. A 
small scale prototype, tested under appropriate cir- 
cumstances, would yield the actual hull dynamics. 
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5.1.3 Design Build Fly Methodology 

The aviation industry traditionally follows the lin- 
ear design process, where a new aircraft is exten- 
sively modeled, high fidelity wind tunnel rigs are 
constructed and tested, and the wind tunnel data 
is post processed and analyzed. The lead time for 
a new aircraft typically takes several years of up- 
front modeling, the wind tunnel models cost around 
a quarter million dollars, running the wind tunnel 
costs nearly $10k/hr, and a team of engineers pour 
over the wind tunnel data to replicate the nonlinear 
dynamics within their simulations. 

The concept of spiral development is catching 
the attention of the aviation industry, where the 
term “Design-Build-Fly” (DFB) expresses the spi- 
ral development process applied to new aircraft de- 
sign. Moving towards a design philosophy that uses 
simplifying assumptions and no wind tunnel data, 
raises concerns about the safety of these new air- 
craft. Reducing the fidelity of the model necessitates 
a comparable improvement in another area. The en- 
abling theory to implement DBF in practice resides 
in adaptive control, which serve two primary roles. 

First, adaptive controllers are built upon esti- 
mation and system identification techniques. When 
applied to a system with uncertain dynamics, these 
controllers take the inputs and measured outputs, 
and adjust the model dynamics until the true sys- 
tem parameters have been obtained. This aids in 
the spiral development process because the proto- 
types are developed on imperfect information, but 
the system identification yields the true parameter 
values which validate or refine the original model. 

Second, adaptive controllers reduce risk while fly- 
ing early prototypes built upon simplified models. 
These controllers operate in real time and tune con- 
troller parameters until the physical system matches 
the dynamics of a desired model. In essence, adap- 
tive controllers start with an imperfect model and 
the prescribed desired response. When there is a 
discrepancy, the controller adapts which adjusts the 
control inputs until the physical system matches the 
desired system response. This greatly reduces the 
risk while flying prototype vehicles because the con- 
troller overcomes model uncertainties. 

The new buzz term for adaptive control within 
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the aviation industry is “Learn-to-Fly” (LTF). This 
describes an exaggerated scenario where the flight 
controller has minimal knowledge of the aircraft dy- 
namics a priori , and the vehicle must “learn to fly” 
before it hits the ground. 

5.1.4 CSR Application 

Of all the potentially new aircraft concepts, the CSR 
system is especially well suited for a DBF approach. 
A lack of aerodynamic data, interactions from the 
tether, flexibility within the slender wing, and the 
multibody structure of the CSR, make the system a 
strong candidate for spiral development. 

Data already exists for many fixed-wing aircraft 
which can ballpark control derivatives for conven- 
tional aircraft design. However, minuscule research 
and aerodynamic data is available for mono- wing de- 
signs, which is needed for CSR modeling. Further- 
more, the satellites spin around the central hub, so 
there is a faster velocity outboard than inboard. It is 
nearly impossible to implement a non-uniform veloc- 
ity profile within a wind tunnel, so no existing data 
or analysis is available to aid the satellite design. 

A vast majority of aircraft are modeled as rigid 
bodies. Some research models flexible wings, but 
these still include a certain amount of rigidity in 
the structure. Very few applications have modeled a 
tether as it moves through space. Simple closed form 
solutions express a catenary model, but they do not 
account for varying drag forces acting throughout 
the length of the tether. 

Finally, the hub has no independent control au- 
thority, so its dynamics are completely dependent on 
the tension within the tethers. A full multibody dy- 
namics formulation is required to capture the equa- 
tions of motion and the tether interactions. Control- 
ling a body in such a manner is extremely novel, and 
not well documented in the literature. 

The CSR project follows a spiral development life 
cycle, because there is very little existing research to 
assist the development process. Small incremental 
steps minimize risk during early development. Each 
cycle helps learn about the system, collects data on 
the current prototype, and refines simulation mod- 
els. More advanced techniques are introduced over 
time, building upon previous lessons learned. 
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5.2 Prototype Vehicles 

True to the spiral development philosophy, several 
prototypes are slated to aid the CSR design and 
analysis. A commercially available control line air- 
craft is the easiest and most economical entry point 
to test new avionics and control theories. A scaled 
mono-wing prototype with carbon fiber construc- 
tion, more closely resembles the full scale vehicle at a 
manageable scale. A full scale prototype is the final 
step assessing eternal flight capabilities of the CSR 
system. At present, the control line aircraft is assem- 
bled and ready for data collection, the scaled model 
is under construction, and no formal plans have been 
established for the full scale prototype. The follow- 
ing sections describe each prototype vehicle and its 
purpose within the spiral development cycle. 

5.2.1 Control Line Aircraft 

The first prototype exhibits the greatest risk. Be- 
cause it is the first vehicle, no data is available to 
contribute towards the design, so all preliminary de- 
sign work is based on unverified simulation models. 
However, the simulation cannot be validated until 
the prototype is flown. Similarly, the avionics hard- 
ware and sensors have never been implemented on a 
real flight platform; only bench tested within a con- 
trolled environment. The first prototype vehicle falls 
into a viscous circle, where it is risky to fly without 
some validation, but the components are only vali- 
dated through flying the vehicle. 

Description 

To minimize the amount of risk with the first pro- 
totype, the selected embodiment is a commercially 
available control line aircraft. Prior to RC radios, 
control line aircraft were the only means for hobby- 
ists to fly. The system has a conventional airplane 
geometry with a gas engine or an electric motor. A 
steel cable tethers the plane to a handle, which is 
held by the operator. The handle has a second line 
that runs out the plane, which manipulates the ele- 
vator control surface. The throttle is run open-loop, 
the plane flies in a circle around the operator, and 
the operator maintains the elevation of he vehicle. 
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Purpose 

This platform was selected as the first prototype 
vehicle, because it is mostly assembled and nearly 
ready to fly out of the box. A few slight modifica- 
tions are needed for the CSR prototype application. 
The control line aircraft is equipped with throttle 
and elevator, so the rudder control surface is the 
only custom addition. Servos and a speed controller 
adjust the control inputs, which eliminates the need 
for a central operator. Finally, wireless serial radios 
communicate between the satellite and the central 
hub, which eliminates the original control line. The 
primary purpose of this vehicle minimizes risk while 
developing the custom electrical system. Utilizing 
a ready made vehicle with stable flight characteris- 
tics is the ideal means to implement the custom and 
untested sensors, avionics, hardware, and software. 

5.2.2 Scaled Prototype 

The first prototype tested avionics and control algo- 
rithms on a stable platform. The second prototype 
takes the next incremental step, and implements a 
system that more closely mirrors the CSR concept. 

Description 

The scaled prototype is a mono-wing design with a 
fuselage at the outboard tip, an elevator outboard of 
the fuselage, and a vertical stabilizer under the fuse- 
lage. The tapered wing starts with zero incidence 
angle at the root and wingtip, and gradually moves 
to a 6.5 ° incidence angle about half way through the 
wing. The wing is carbon fiber Kevlar construction 
which measures seven feet long. The entire horizon- 
tal surface is deflected for pitch control, and 25% of 
the vertical stabilizer is articulated for yaw control. 

Purpose 

The purpose of the second prototype is to learn 
the characteristic aerodynamics of the satellite ve- 
hicle. The earliest models have simplifying assump- 
tions which did not account for all the aerodynam- 
ics. Wing tip vortices and interactions between the 
horizontal stabilizer and elevator were not modeled. 
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The high aspect ratio wing has some degree of flex- 
ibility, which needs to be observed. Other aerody- 
namic features are unique to the CSR system, and 
not completely known. The rotating system has dif- 
ferential velocity over the wing surface, which may 
impact the relative heading. This prototype includes 
the tapered wing which achieves an elliptical lift dis- 
tribution. Upwash from the wingtip produces lift 
on the horizontal stabilizer surface. These are sig- 
nificant aerodynamic effects that contributes to the 
eternal flight capability, and should be accounted for 
in the simulation model. The goal of this prototype 
is to observe the dynamics of the scaled system, and 
revise the simulation model to more accurately re- 
flect these aerodynamic features. 

5.2.3 Full Scale System 

No formal plans have been established for a full scale 
prototype. The scaled prototype will shed valuable 
insight into the vehicle dynamics, which will lead to 
an improved model and a better understanding of 
how the full scale system behaves. However, to com- 
plete the analysis and to fully assess the eternal flight 
capabilities, this prototype will need to be built. 

Description 

The following description is based entirely on prelim- 
inary analysis work, and is subject to change once 
data is collected and evaluated from preliminary pro- 
totype models. The full scale CSR is expected to 
have around a 350 foot diameter. The satellite wing 
has a high aspect ratio greater than 50, with a 0.7 
foot chord and a length of 35 feet. Control surface 
positions and sizes may change once new data is re- 
ceived, but are anticipated to be in nearly the same 
positions with the same relative sizes as the scaled 
prototype. Additional analysis is needed to specify 
the motor, propeller, batteries, and solar cells. 

Purpose 

The most obvious purpose of the full scale proto- 
type is to demonstrate eternal flight capabilities. 
While, the first prototypes validated hardware and 
collected aerodynamic data, this vehicle is primarily 


designed to attain eternal flight. This vehicle follows 
the same testing procedure as the other prototypes 
which yields valuable information about the vehicle 
characteristics. As a worst case scenario, the col- 
lected data will shed insight into what revisions are 
needed to improve the performance of the vehicle; 
and if all goes well, the CSR will be the first aircraft 
in history to achieve eternal flight. 

5.3 Testing Plan 

Each prototype is a natural progression that moves 
from simple vehicles to more complex systems. For 
each prototype, a similar procedure is used which 
minimizes risk by first demonstrating fundamental 
capabilities, and then advances to more complex ma- 
neuvers. This section describes the three general 
types of tests: (1) round the pole, (2) constrained 
flight, and (3) free flight. 

5.3.1 Round the Pole Flight 

This test has the least amount of risk because it only 
involves a single satellite. These tests gather impor- 
tant data that is used for system identification and 
inner-loop controller design. 

Description 

The Round the Pole (RTP) flight most closely re- 
sembles the control line aircraft. A fixed pole an- 
chors the system to the ground, but the pole has a 
bearing that is free to swivel about the vertical axis. 
The rotating hub has the data collection board, the 
transceiver for the RC radio, wireless telemetry to 
communicate with the satellite, and a magnetometer 
and rate gyro to measure heading and angular rate 
of the hub. A tether is attached to the hub and a 
single satellite vehicle. The satellite has throttle, el- 
evator, and rudder control inputs which are obtained 
from the central computer. 

Objectives 

For each prototype, this test marks the maiden voy- 
age, so the emphasis is on reducing risk to the 
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untested vehicle and gaining confidence in the sys- 
tem capabilities. All work has been purely analytical 
and performed in simulations, so these flights serve 
to validate and refine the computer models. The 
main goals of these tests are to focus on the inner 
loop-control mode functionality by: 

• demonstrating the new system in flight, 

• validating all hardware, sensors, and actuators, 

• finding inputs to achieve static equilibrium, 

• characterizing system dynamics through sys- 
tem identification, 

• comparing flight data to simulations, and 

• validating or refine existing model dynamics. 

Approach 

Testing a new prototype is done through incremental 
steps to make sure more risky maneuvers are built 
upon a firm foundation. The first step will bench 
test all the prototype components. All control in- 
puts will be tested to make sure they work properly. 
All sensors will be evaluated to confirm good mea- 
surement, and that the data logger works properly. 

Taxi Test : After the bench test, move on to the 
taxi test. This has the system in a ready to fly con- 
figuration, and throttles up the motor to around 90- 
95% of the value needed to achieve lift. The vehicle 
will spin around the hub at nearly flight speeds, but 
will not be enough to achieve lift. As the satellite 
runs around the pole, data logging is enabled, sen- 
sors acquire data, and the control surfaces are per- 
turbed slightly to evaluate their effectiveness before 
the vehicle ever leaves the safety of the ground. Post 
processing will ensure the sensor data looks good, 
and changes will be made to tune the controller val- 
ues and RC radio settings. 

Trimming: The next series of tests will run the 
system open-loop controlled by a pilot. The vehicle 
will throttle up until lift is attained and the pilot 
will actuate the control inputs. The RC radio is 
equipped to trim the inputs onsite, so adjustments 
will be made to find the trimmed RC input values. 

Doublets: With trimmed RC inputs, the pilot 
will command a series of doublets on each control 
channel. Doublets are a quick maneuver designed 
to perturb a single input and then quickly bring the 


vehicle back to a steady-state. The doublets begin 
with small magnitudes and gradually move up to 
larger command inputs. At some point the doublet 
will be deemed “too aggressive” and the test is com- 
plete. The data is post processed and the values 
of the most aggressive doublets are used to design 
control gains that are appropriate for the system. 

Feedback Control : The open-loop flight data 

sheds light on the system dynamics. This informa- 
tion is used to determine the gains needed for a feed- 
back control algorithm. The system is now flown 
closed-loop with the feedback controller applying in- 
put commands. Data is collected and the feedback 
control gains are adjusted as needed. 

Signal Injection: Once the satellite has a solid 
feedback control law, the final step is to evaluate the 
system through multisine signal injection. This pro- 
cess flies the feedback controlled system and super- 
imposes a summation of sinusoidal signals on top of 
the existing control inputs. These signals excite the 
system dynamics which generate new measurements 
that are used in the system identification process. 
The data is post processed and compared against 
the the existing simulation model, and revisions are 
made as needed, until the computer model matches 
the collected data. This revised computer model is 
used to develop the controller for the next round of 
testing: the constrained flight. 

5.3.2 Constrained Flight 

The RTP tests investigated the behavior of a single 
satellite vehicle. The constrained flight tests increase 
the level of difficulty by working with three satellites 
in unison. However, this series of tests still reduces 
the amount of risk by constraining the central hub. 

Description 

The constrained flight test setup is very similar to 
the RTP setup, except three satellites are tethered to 
the central hub. The satellite system is setup iden- 
tically to the RTP system with all the same com- 
ponents and control inputs. The central hub has all 
the same components but three satellites are teth- 
ered to the hub. The first series of test has the 
hub constrained identically as the RTP system with 
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one single degree of freedom which rotates about the 
vertical axis. Subsequent tests relax the number of 
constraints where first, the hub is free to translate 
vertically, and then the hub is bounded but uncon- 
strained in all six degrees of freedom. 

Objectives 

The RTP tests yielded insight into the dynamics of a 
single satellite which mirrors the inner-loop control. 
To control the complete CSR system also requires 
coordination between all the satellites to achieve the 
outer-loop control. The objective of these tests is 
to evaluate the effectiveness of the controller to im- 
plement coordination between all entities. It begins 
with the most basic task of maintaining equal spac- 
ing between satellites. Then it graduates to coordi- 
nated flights, first with a single degree of freedom, 
then with two degrees of freedom, and finally with 
bounded six degree of freedom tests. The main goals 
of the constrained flight tests are to: 

• demonstrate equal spacing between satellites, 

• single degree of freedom tests, 

— collective throttle and elevator 

— cyclic rudder and elevator 

• two degree of freedom tests, and 

— collective throttle and elevator 

— cyclic rudder and elevator 

• bounded unconstrained tests 

— collective throttle and elevator 

— cyclic rudder and elevator 

Approach 

These tests follow a similar progression to the RTP 
tests. All satellites are tuned using previously col- 
lected data and bench tested to confirm all hardware, 
software, and sensors are working properly. 

Taxi Test : The first test is a taxi test. Just like 
before, all the satellites will throttle up to around 
90-95% of their cruising speed and will taxi around 
the pole. The feedback control will be active on all 
three satellites. While taxiing the satellites will have 
the same elevation, so there will be minimal elevator 
commands. There is no outer-loop command that 
adjusts the relative heading, so there will be min- 
imal rudder commands. The main demonstration 


will focus on the feedback controlled throttle com- 
mands to determine if the control law can maintain 
equal spacing between satellites and track a given 
hub angular rate. 

Equal Spacing: Once the taxi test has been suc- 
cessfully completed, the next incremental step takes 
the vehicles off the ground. The hub is still con- 
strained in all DOF except for vertical rotation. The 
reference signal increases the angular rate of the cen- 
tral hub, and the controller will ramp up the throt- 
tle to track the increased angular rate. Once the 
satellites attain cruising speed they will lift off. The 
feedback controller is responsible for all the inner- 
loop commands for each satellite, so the operator 
is solely responsible for adjusting the hub angular 
rate. The collected data will be analyzed and evalu- 
ate how effective the system maintained equal spac- 
ing. Any deficiencies or errors will be simulated to 
observe how those errors would impact an uncon- 
strained system. If the data reflects positively, then 
the next round of testing begins. 

Collective Commands: With equal spacing ac- 
complished, the next tests perform more elaborate 
coordinated flight maneuvers that will be used for 
outer-loop control. The hub is still constrained to 
a single degree of freedom while the four outer-loop 
modes are evaluated. The tests begin with collec- 
tive modes because they are easier to observe. The 
first test is collective throttle because it most closely 
resembles the previous set of tests. The system is 
ramped up to hover equilibrium state, and a ref- 
erence command increases the angular rate of the 
central hub. When the satellites exceed their cruise 
speed, additional lift is generated and the vehicle will 
increase altitude. The next test runs collective eleva- 
tor commands. When the system is in equilibrium, 
the operator increases the elevator reference signal, 
and the satellites begin to climb together. Collected 
data will evaluate whether the satellites maintain the 
same elevation and attitude during their ascent. 

Cyclic Commands: After demonstrating collec- 
tive commands, attention turns toward cyclic com- 
mands. Cyclic rudder is more stable that the eleva- 
tor, so begin there. The system is brough to hover, 
and the operator commands a cyclic rudder com- 
mand. Because the pole constrains the system, there 
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will be very little to observe onsite. However, the 
collected data should reveal whether the satellites 
change their relative headings in a cyclic manner. 
The final one DOF test is cyclic elevator which tilts 
the entire plane of rotation. Once in hover, the oper- 
ator applies the command and the satellites will al- 
ter their elevations cyclically. The data should show 
that each satellite attains a low and high point at 
the same locations. It should also show that each 
satellite achieves the same elevation changes while 
navigating around the hub. 

Relaxed Degree of Freedom: After successfully 
demonstrating coordinated flight with a single DOF 
hub, the next series of tests relax the number of con- 
straints on the hub. The pole hardware is replaced 
so the hub is free to both rotate and translate about 
the vertical axis. The testing process is repeated for 
each of the four outer-loop control modes to make 
sure the feedback control is working as expected. 
The two DOF hub will shed more useful insight on 
the collective commands because they both influence 
the vertical translation of the system. 

Bounded Free Flight: Demonstrating the hori- 
zontal translation is accomplished with the bounded 
free flight test. This test replaces the central hub 
connection. Rather than a bearing mounted on 
the pole, the central hub is fashioned as an annu- 
lar ring. The ring is larger than the pole diameter 
so the hub is free to roll, pitch, and translate hor- 
izontally. This closely resembles unconstrained free 
flight dynamics, but still has the added security of 
placing a boundary on the central hub. Again, each 
of the four outer-loop control modes is tested un- 
der these new constraints, and data is reviewed to 
make sure the system is well behaved. To further 
reduce the amount of risk, moving into a fully un- 
constrained free flight, several rings of various di- 
ameters can be constructed to adjust the boundary 
distance. Smaller rings offer more security by limit- 
ing the range of motion. Once the system response 
behaves as expected, larger rings can showcase the 
unconstrained system dynamics. 

5.3.3 Free Flight 

Previous testing paved the way for the pinnacle of 
the prototype tests: the completely unconstrained 
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free flight of the CSR system. This series of test 
marks the completion of each prototype, so the in- 
formation gathered here is intended to aid in the 
next prototype development. 

Description 

These tests are run on the complete CSR vehicle. 
The training wheels have been removed, and it is 
time to see the big picture at work. All the avion- 
ics, hardware and sensors are the same for the satel- 
lite vehicles. The annular ring for the central hub 
has been replaced with the preferred embodiment of 
a convex cylindrical housing. The hub components 
are identical to the previous testing. Early tests will 
take off with the satellite arms fully extended, but 
more sophisticated prototypes may include a tether 
retraction mechanism for take off and landing. 

Objectives 

Unlike previous tests which were focused on devel- 
oping capabilities needed for the control system, the 
objectives of the free flight tests are more specific 
to the individual prototype being tested. The focus 
of this testing is to acquire enough insight to aid 
the development of the next prototype, which will 
have different objectives from one prototype to the 
next. However, there are many commonalities that 
are shared among all the prototypes. Testing ro- 
bustness to disturbances, waypoint navigation, and 
controller optimality are all accomplished with the 
free flight tests. The major goals of the free flight 
tests are to: 

• evaluate system robustness to disturbances, 

• demonstrate waypoint navigation, 

• implement optimal control techniques, 

• assess effectiveness of adaptive control, 

• measure dynamics of outer-loop control modes. 

Approach 

Now that the system is flying on its own, this set of 
testing focuses on pushing the boundaries of the ve- 
hicle capbilities. Determining what the current sys- 
tem is capable of is necessary information to progress 
to the next prototype design. 
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Disturbance Rejection: After flying with the 

proverbial safety net from the constrained tests, the 
first order of business will assess how robust the sys- 
tem is to disturbances. Wind gusts are the most 
notable disturbances the system must reject. How- 
ever, the size of the vehicle makes outdoor flying a 
necessity, which is impossible to create a completely 
controlled environment. Furthermore, generating a 
static wind disturbance that acts on the entire vehi- 
cle is infeasible. As an alternative, a series of con- 
trolled experiments is presented. Consider a wind 
gust acting on the nose of a satellite. This creates 
a force that acts longitudinally. Force can also be 
applied longitudinally by adjusting the throttle com- 
mand. To test this type of disturbance, code is added 
which scales up or down the motors outside the con- 
trol loop. The control law is unaware of the motor 
scaling, and only acts on the measurements which 
indicate a tracking error. A similar test may evalu- 
ate the elevator channel by dropping a mass in the 
middle of a flight. The decreased payload simulates a 
force acting in the vertical direction, which the eleva- 
tor must mitigate. These experiments can be easily 
modeled in the simulation environment, so collected 
data may further validate the computer model. 

Waypoint Navigation: Knowing how much dis- 
turbance can be rejected leads into the next round 
of testing: waypoint navigation. Thus far, every- 
thing has focused on maintaining a hover equilib- 
rium state. This test will assess the system dynam- 
ics outside of the hover condition. The test is ex- 
tremely straight forward; start the vehicle in one 
location and command it to go to another point. 
The controller will be designed to have some pre- 
scribed response characteristics, like rise time, over- 
shoot, damping ratio, etc. The goal is to compare 


the physical response to the prescribed dynamics and 
confirm that the controller is implementing the de- 
sired closed-loop dynamics. 

Optimal Control : Once the vehicle is able to nav- 
igate to a new position, the next most pertinent issue 
is how efficiently it gets there. Early controllers use 
pole placement which exactly specifies the charac- 
teristic response. This strategy was used to reduce 
the risk with early flights by giving the controller 
enough authority to overcome any unexpected dis- 
turbances or modeling errors. However, this added 
security comes at the price of increased control en- 
ergy. For an eternal flight vehicle, minimizing the 
amount of control effort is a primary consideration. 
Rather than using a pole placement feedback con- 
troller, the next test will implement an optimized 
control law where the gains are obtained through a 
linear quadratic regulator (LQR) design. The ex- 
periment runs the system with the same waypoint 
navigation, but collects data on the control signals 
used to navigate to that point. 

Adaptive Control Assessment: Adaptive control 
is an important control technique to mitigate model- 
ing uncertainty for both the prototype development 
and for the final full scale CSR vehicle. Understand- 
ing the capabilities and limitations of the adaptive 
controller is important to ensure stability and avoid 
unacceptable transient responses. Previous system 
identification tests yield the actual plant model dy- 
namics. For this test, the controller gains are de- 
signed with slightly erroneous values. The controller 
should recognize that there is modeling error and 
adapt to the correct characteristics. Because the cor- 
rect vehicle dynamics are already known, the adap- 
tation can be evaluated against the true parameters. 



Chapter 6 

Conclusions 


The Centrifugally Stiffened Rotor (CSR) concept 
is a unique configuration approach that has never 
previously been attempted, which could fundamen- 
tally lead to a new approach for High Altitude Long 
Endurance (HALE) or ATmospheric Satellite (AT- 
Sat) missions. AT-Sats map into many missions that 
existing GEO and LEO satellites simply can not per- 
form, and in particular AT-Sats align with missions 
were communication transmission time needs to be 
minimized for such capabilities as omni-present wire- 
less broadband. 

The CSR significantly improved structural and 
aerodynamic efficiency open up the design space to 
become more feasible, with initial analysis indicat- 
ing that feasibility may be possible through this con- 
cept approach in combination with energy storage of 
about 400 Whr/kg and thin- film solar cells around 
35%. Feasibility relates to accomplishing year round 
missions at high U.S. latitudes. 

The Phase I research focused on developing anal- 
ysis tools capable of capturing the unique attributes 
of this multibody approach, including aerodynamic, 
structural, force balance, and power required. Per- 
formance comparisons were made with other HALE 
tool sets used in the analysis of DARPA Vulture con- 
cepts, as well as calibrations of the new tools to the 
most similar concepts. 

The initial performance analysis indicates a re- 
duction in power required on the order of 35% com- 
pared to the prior QinetiQ Zephyr 7 HALE en- 


durance record holder, with the Zephyr carrying less 
than a 5% payload fraction compared to the CSR 
concept with a 10% payload. 

A number of compelling missions have been iden- 
tified that map directly into the unique capabilities 
of this advanced concept. Missions include: 

Multi-Functional Airborne Wind and Surveil- 
lance Commercial Platforms at Lower Altitudes: 
Aerial platforms that operate at altitudes up 
to 2000 ft altitude (without FAA impediment 
of operational feasibility) that can both capture 
wind energy more effectively than ground-based 
wind turbines, and provide close proximity surveil- 
lance/communications with a coverage diameter of 
50 miles. The CSR concept has the potential to 
achieve a higher Lift/Drag ratio compared to other 
airborne wind concepts, which results in higher 
tether angles and less land area underneath the ve- 
hicle radius of operation that depends on incoming 
wind direction. Off-shore application is particularly 
appealing since a CSR aerial vehicle would eliminate 
the need for a large /expensive mooring platform, 
which is required for ground-based wind turbines. 
The ability of this concept to have a non-moving 
tether from the ground to the center hub permits the 
inertial connection to be less complex than current 
airborne wind turbines. Potentially this mission con- 
cept could be applied all the way down to the level 
of distributed residential power production, with a 
sUAS sized version likely capable of providing 2 to 
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5 kW of power at altitudes of 500 feet with average 
wind speeds of less than 20 mph. 

Distributed Aperture AT- Sat Observatory: The 
CSR concept provides a large rotating structure, as 
well as large volume at the center hub, with the ca- 
pability of operation above the majority of the atmo- 
sphere. Integrating high resolution, compact, linked 
optic sensors (such as the Low Mass Planar Photonic 
Imaging Sensor currently being developed as Phase 
II NIAC research) at the tip and root of each rotor- 
wing will permit a rotating imaging array across a 
full 600 diameter azimuth. Combining this with a 
lower resolution conventional imaging system in the 
center hub for image filling could provide a resolution 
never before possible. Such an observatory could 
be designed as a dual-purpose system that provides 
imaging both upwards for space investigations, as 
well as for Earth imaging. 

Extra- Planetary Compact VTOL Exploration 
Platforms: A version of the CSR concept has been 
identified that enables full retraction of the rotor- 
wing to achieve a highly compact VTOL vehicle 
that can operate at extremely low discloading (<0.03 
lb/ft2) on atmospheres such as Mars that have at- 
mospheric densities similar to Earth at 100,000 feet 
altitude. Because it was highly uncertain whether 
the centrifugal stiffening loads would be sufficient 
to maintain a rigid rotor-wing for this highly faceted 
spanwise design, a sub-scale version of the rotor- wing 
was fabricated and tested with success. The rotor- 
wing experiences an increase in weight and complex- 
ity with a decrement in performance to achieve com- 
pact packaging, however was shown to be a feasi- 
ble approach that could facilitate very low discload- 
ing renewable VTOL aerial vehicles which could be 
packaged effectively in an aeroshell. 


Since the performance analysis indicates a com- 
pelling advantage for this concept, the Phase I re- 
search focused on addressing the key risk area for 
this concept, establishing that vehicle control is fea- 
sible for this multibody dynamics problem. A com- 
plete derivation was presented which outlined how 
the CSR is modeled, the dynamics that dictate the 
equations of motion of the system, and a means of 
trimming and linearizing the system. This transi- 
tioned into the theory needed to control this partic- 
ular vehicle, and simulations demonstrated the con- 
troller response implemented on the nonlinear sys- 
tem. The key hurdle to overcome from the controls 
point of view, is maintaining control when the exact 
relative position of the satellite is unknown. Estima- 
tion and reduced order estimators were shown to be 
able to overcome this, because the CSR vehicle was 
determined to be both fully controllable and fully 
observable. Furthermore, the system dynamics can 
be arbitrarily selected to meet any desired response 
within the limits of the actuators. Finally, parame- 
ter estimation and adaptive control techniques can 
be applied to this type of vehicle. 

A sub-scale flight demonstrator at the 30 lbf 
gross weight size has been initiated to validate the 
ability of the control system to accommodate real 
world disturbances as well as winds and maintain 
stable flight. Additionally, the sub-scale prototype 
permits other study assumptions to be addressed, 
such as the mechanical complexity of the tether reel- 
in/out system which has not been analyzed. This 
demonstrator effort will continue as part of the PhD 
dissertation for the lead control researcher, and be 
completed prior to a follow-on proposal for Phase II 
research. 
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Abstract — A comprehensive report outlined the model- 
ing and simulation work performed for the Centrifugally 
Stiffened Rotor (CSR) concept vehicle. Since that time, 
the majority of the research effort focused on developing 
the avionics hardware that serves as a flight controller for 
the CSR prototype and full scale vehicles. This document 
outlines the work performed and the outcomes from that 
research effort. 

This summary describes the progress with the avion- 
ics development. The hardware is run on a BeagleBone 
Black single board computer, which runs the Linux 
operating system. It describes the Device Tree Overlays 
(DTO) that are available to configure the pinmux set- 
tings, the General Purpose Input Output (GPIO) pins, 
the UARTs which allow for wireless communication 
between the satellites and the Ground Control Station 
(GCS), the set of Programmable Realtime Units (PRU) 
which deliver realtime 5 nanosecond clock cycles, and 
a Xenomai cokemel which executes timing commands 
in real time. 

I. Platform 

The controller for the CSR vehicle runs on a single 
board computer called the BeagleBone Black (BBB). 
This board runs linux OS and has two header rails 
for fast prototype development. The Linux kernel on 
the BBB was updated to 3.14 which includes PWM, 
I2C, and SPI interfacing, and makes the device tree 
overlay (DTO) easier to implement. Each pin on the 
board can serve different functions based on how the 
DTO is loaded at start up. The default setting uses many 
pins for HDMI and eMMc chip allocations, which are 
needed for the avionics development. Those pieces of 
hardware were disabled so their respective pins could 
be used for other purposes, like serial communication 
and PWM input and output. 

II. Xenomai Cokernel 

A crucial step for the avionics hardware is ensuring 
the Linux OS runs in realtime. When running an 


avionics controller, reliable timing is essential because 
the radio inputs and motor outputs cannot be delayed. 
This is accomplished through Xenomai; a type of co- 
kernel which runs along-side the actual kernel, and 
properly manages system interrupts to make the Linux 
OS adhere to hard realtime constraints. In the past few 
months, the embedded Linux community has come a 
long way in making Xenomai accessible by integrat- 
ing Xenomai into the general kernel building popu- 
lation. This eliminated a lot of compatibility issues, 
and greatly reduced the number of patches to enable 
Xenomai. The BBB is running Xenomai, so a low level 
OS commands will not interfere with reading a sensor 
or setting a motor value. 

III. Device Tree Overlays (DTO) 

Each pin on the BBB can be assigned different func- 
tionality which makes the board extremely modular. 
A Device Tree Overlay (DTO) is a type of code that 
describes how all the pins are configured on the board. 
The overlay describes what the pin is sued for and 
assigns appropriate settings for the particular hardware 
it is associated with. The Device Tree Compiler (DTC) 
is a special type of compiler that convert the human- 
readable source code into the binary object. For this 
project, a DTO was created which enabled four UARTs, 
two I2C, two PRU, four LEDs, and several GPIOs. 

IV. General Purpose Input/Output (GPIO) 

A general purpose input/output (GPIO) pin does not 
perform any form of communication protocol, and it is 
not associated with any particular piece of hardware. 
Rather, it is a very simple pin that can be configured 
to either read a voltage as an input signal, or to set a 
voltage as an output signal. The avionics board uses 
four LEDs to communicate the status of the system, 
and it has a button input enabled to read user inputs. 

V. UART Serial communication 

The CSR will need to communicate between each 
of the three satellites and a base station, so wireless 



communication is required. The BBB has several com- 
munications protocols, but UART was selected for this 
application because it is the protocol used by XBee RF 
radio systems, which have very reliable signals within 
a small form factor. Four UARTS were enabled and 
communicating to each other on the BBB. Thus, there 
are enough serial ports configured to run the system 
and maintain communication with home base. 

VI. MPU9150 9DOF IMU Sensor 

To implement the controls algorithm, the vehicle 
must be aware of its state. This is accomplished with 
a nine degree of freedom (DOF) inertial measurement 
unit (IMU), called the MPU9150 from InvenSense. It 
has three axis output for magnetometer, rate gyro, and 
accelerometer readings. This board has the added ben- 
efit of a 10th channel that can be used to read altitude 
readings from an external sensor, and an additional 
I2C bus so a second IMU (for redundancy and data 
fusion) can communicate over the same channel. It 
can output attitude in both Euler (nice for debugging) 
and quaternion (better for actual application to prevent 
gimbal lock) formats. The sensor has been connected to 
the BBB where raw data values are processed through 
a data fusion filter algorithm to yield the vehicle state. 

VII. Processing Realtime Units (PRU) 

Another major task requires reading inputs from a ra- 
dio, and sending outputs to electronic speed controllers 
(ESC). Both signals use a form of pulse width modula- 
tion (PWM), which a type of square wave signal with a 
fixed period and a variable duty cycle. Accomplishing 
this required accessing the two Processing Realtime 
Units (PRU) which are dedicated chips that runs in 
realtime (and not subject to delays from the Linux 
userspace OS) with a 5 nanosecond clock cycle. This 
project requires both subsystems to read radio PWM 
inputs and to send high fidelity PWM outputs to the 
ESCs. Working with the PRU requires the assembly 
programming language, which is a set of about 30 
commands that manipulate individual bits directly on 
the chip’s register, hence the nickname ’’bit banging”. 
Following TI’s instruction set, both PRU subsystems 
were enabled, able to load the assemble code, able to 
control both the input and output registers, and able 
to configure the register memory allocation to send 
commands from user-space down into kernel-space. 

VIII. PRU Based Servo Outputs 

Assembly code has been developed to produce PWM 
outputs. One register is used to maintain the total 


number of steps for the duration of a PWM pulse. Eight 
more registers are used as counters for the duration 
of a particular channel. At the start of a new loop, 
all the output pins are set to high, and the code runs 
through a loop which decreases the count each loop. 
When a counter decrements to zero, the output pin of 
that channel is set to low. AT the pins are reset back 
to high at the start of a new loop. 

IX. PRU based Radio Inputs 

Assembly code has been developed to read PWM 
inputs. One register is used to maintain the bits that 
describe the previous state of a channel. Eight more 
registers are used as counters for the duration of a par- 
ticular channel. Assume an input channel is currently 
low, and the counter is set to zero. When the signal goes 
high, the counter increments upwards. At some point 
the signal will drop low. If the previous state was set to 
high, the code passes out the counter value, clears the 
counter value, and resets the bit for the previous state. 
The output value is used to determine the value of the 
input coming in from the radio. 

X. Code Development 

All the individual pieces of hardware need to be 
connected together under a single set of code. The 
avionics code initializes the vehicle, runs a continuous 
loop during the flight, and shuts down the system when 
prompted. The start up code sets all motor values to off, 
initializes the sensors, and clears the buffer for data 
logging. The code looks for several different user input 
commands that specify different flight scenarios (data 
collection, system identification, etc), and commands to 
arm the motors. The main loop has several sub-loops 
that run at different frequencies. The fastest frequency 
is for IMU data, controller commands, and ESC outputs 
to the motors and servos (this is the inner loop control 
for the satellites). The next loop monitors the hub state 
and waypoint navigation (this is the outer control loop 
for the vehicle). The slowest loop is setup to handle 
data logging, performance monitoring, and any other 
non-critical tasks. The code also has an exit command, 
which disables the motors, completes the data logging 
process, and signals the end of the cycle through an 
external LED. The entire project is developed in the 
C programming language to maintain fast execution 
times. The major functions are compartmentalized into 
static libraries so that features can be easily ported to 
other other prototype vehicles. 



XI. Bench Testing 

All the individual capabilities have been transferred 
over onto a single development board. The data log and 
realtime constraints were tested together, where running 
empty code (no controller, data logging only) was able 
to achieve 1000 Hz sample rate, and 30 minutes of 
data logging used less than two percent of the available 
disk space. The radio and speed controllers were tested 
together, where the PRU signals were passed through 
logic level converters to convert between the 3.3V used 
by the BBB and the 5.0V used by the equipment. The 
IMU was tested with the data fusion algorithm, where 
raw data was compared to the fused data, and gain 
values were adjusted to optimal values to overcome 
sensor noise. Finally, a simple PID controller was 
developed to evaluate a complete system which com- 
bines, inputs, outputs, states, and data logging into a 
single demonstration. The controller makes appropriate 
adjustments based on the state and control inputs, so 
the avionics board development is complete. 
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STAR-CCM+ analysis 

• Unstructured mesh 

• Unless otherwise noted: 

- SST (Menter) k-cu turbulence model 

- y-Reg transition model 

• Actuator disk for propellers 

- Radial thrust/torque follows Goldstein optimum 

- Thrust specified; torque and RPM calculated via 
supplied advance ratio table 
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Outline 


• Blown wing analysis (takeoff conditions) 

• Cruise wing analysis (unflapped and 
unblown) 

• Blown wing design space exploration 

• Nacelle design 

• Propeller analysis 
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Introduction 


• "Effective lift" includes the vertical component of the propeller thrust 

• Results labeled "FUN3D" were produced by NASA using the FUN3D CFD code 
with the Spalart-Allmaras turbulence model, and are included for comparison. 

• Unless otherwise noted, these simulations employ the SST (Menter) k-co 
turbulence model and the y- Re 0 turbulence model. 

• Other turbulence model results included for comparison are SST (Menter) k-co 
fully turbulent (labeled as "k-co"), and Spalart-Allmaras with and without 
curvature correction (labeled as "SA" and "SA+CC," respectively). 

• The blown wing mesh size is roughly 12 million cells, except for the results 
labeled "fine," which use a 24-miifion cell variant of this mesh. 

• Unless otherwise noted, these results are with flaps deployed to 40° 

• Unless otherwise noted, "blown" indicates the propellers are powered at 300 hp 
(shaft). 
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Baseline design 




a 

Power 

Lift 

Drag 

Thrust 

Q 

Effective Q 

6° 

300 hp 

2,810 lb 

479 lb 

818 lb 

4.10 

4.15 

7° 

300 hp 

2,890 lb 

505 lb 

818 lb 

4.22 

4.27 

6° 

240 hp 

2,620 lb 

421 lb 

685 lb 

3.83 

3.87 

7° 

240 hp 

2,680 lb 

445 lb 

685 lb 

3.92 

3.97 

9° 

240 hp 

2,520 lb 

510 lb 

685 lb 

3.67 

3.76 
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Baseline design 


Lift distribution, a = 7° 
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Final design 
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61 KEAS, 300 hp, 40° flaps 


Blown (effective) 


•Blown 30° flaps 


Blown (effective) 


Blown 30° flaps 


•Unblown 


•Unblown (no nacelles) Unflapped 


•Unblown 


•Unblown (no nacelles) Unflapped 
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61 KEAS, 40° flaps 



Blown 

Unblown 
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Lift distribution, 1 0°, 300 hp, 40° flaps, 61 KEAS 



Spanwise position (ft) 
— Wing+flap —Wing only 
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Lift distribution, 11°, 300 hp, 40° flaps, 61 KEAS 



Spanwise position (ft) 


— Wing+flap —Wing only 
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Lift distribution, 1 2°, 300 hp, 40° flaps, 61 KEAS 



Spanwise position (ft) 


— Wing+flap —Wing only 
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a=10° 



Pressure Coefficient 

-13 -11 -9 -7 -5 - 3-11 3 
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Pressure Coefficient 
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a-l°, unblown and flaps retracted 
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61 KEAS, unblown, 40° flaps, a- 4° 
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300 hp, 40° flaps, a= 10° 

Comparison to preliminary experimental results 


4.000 
3,500 

3.000 

^ 2,500 

_g 

^ 2,000 
^ 1,500 

1.000 
500 

0 



10 20 


60 70 


Thrust 

-X- Effective lift 
-^-Effective lift (experimental) 
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Cruise wing analysis 
(unflapped and unblown) 
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Cl ; a=0° 


Cl*c/cref ; a=0° 


Cl ; a=8° 


Cl*c/cref ; a=8° 



Spanwise location (ft) 
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Comparison to takeoff conditions 
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Pressure Coefficient 



- 2.2 


- 3.0 



a=8° 








0.02 


0.04 
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a=0° 


Airfoil locations 
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Nacelle locations 
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a=8° 


Airfoil locations 
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Nacelle locations 
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Cruise aircraft component analysis 

Fuselage+VTail only Fuselage+VTail only No nacelles Whole aircraft 
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Specifications 


Wing area 

148 ft 2 

Wingspan 

45.6 ft 

Aspect ratio 

14 

Airfoil 

NACA 23015 

Washout 

3° 

Flap 

30% chord 40° plain flap from y=2.5 ft to y=15 ft 
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Case descriptions 


Baseline 

69 kts 

Taper 

0.6 

Power 

P/qSV m 

Propeller diameter 

Propellers 

T/qS 

T/qA 

Blown span 

A 

69 kts 

0.6 

797 hp 

1.58 

3.35 ft 

6 

0.80 

2.25 

44% 

B 

69 kts 

0.6 

797 hp 

1.58 

4.17 ft 

6 

0.92 

1.67 

55% ◄ — entire flap span 

C 

69 kts 

0.3 

797 hp 

1.58 

3.35 ft 

6 

0.80 

2.25 

44% 

D 

69 kts 

0.3 

1,613 hp 

3.19 

3.35 ft 

6 

1.20 

3.38 

44% 

E 

80 kts 

0.3 

797 hp 

1.01 

3.35 ft 

6 

0.57 

1.60 

44% 

F 

69 kts 

0.3 

531 hp 

1.05 

3.35 ft 

4 

0.54 

2.25 

29% 

HEIST 

(for comparison) 

61 kts 

0.5 

300 hp 

2.31 

1.465 ft 

18 

1.23 

2.23 

85% 
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Case relationships 

more power 

^ — J d 



blowing 

Baseline 1 — — — A 

more taper higher 

1 — — — — ►( C 1 — m E 


larger 

propellers 

fewer 

propellers 

— — >m f 
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Note that E is at a higher V m , so similar Qs will not correspond to similar lift 





Case B 
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Case A 
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Cases C-E 
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Case F 
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Nacelle design 
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Nacelle positioning (Q=1.2) 


Contours show the streamwise (X) component of the product of the surface normal and the 
local surface pressure, in Pascals. 
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Comparison of CFD and experiment 

3- and 5-blade designs 
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ABSTRACT 


Research in Flight Company has conducted an aerodynamic analysis of the D8.5 outer-mold line 
geometry for subsonic flight conditions using the proprietary surface mesh based solver known as 
FlightStream™. Specifically, three major tasks were performed and this report on the project has 
been prepared and presented at the NASA Langley Research Center. 

The analysis scope included the validation of FlightStream™ software for high-lift configurations 
involving flaps. The comparison of aerodynamic results to wind tunnel data for a simple slotted 
flap on a General Aviation Wing model and for the two flapped configurations associated with the 
American Institute of Aeronautics and Astronautics High-Lift Prediction workshop has been 
performed. 

The baseline geometry used for the study involved the addition of single-slotted flaps on the D8.5 
geometry. The configuration is tested at a low Mach number (M=0.2) for an angle of attack sweep 
for 3 different flap deflections, 15°, 25°, and 35°. A motion space for the flap for different 
deflections has been created using variables in the longitudinal and vertical direction and used as 
a design space for the trade study. Aerodynamic results have been evaluated at different positions 
within this motion space for all three deflection angles of the flap. The goal of the trade study was 
to identify near-optimal slotted flap gap size. 
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1. INTRODUCTION 


The D8 aircraft concept is one of several high-performance designs being considered for the next 
generation of transport aircraft. Several studies have been previously commissioned to investigate 
its flight performance in cruise configurations. The present effort is focused on the analysis of the 
D8.5 model in high-lift configurations using flaps. Since a design for the flap to be used on this 
model does not currently exist, this effort encapsulates the design and integration of a simple 
slotted flap to the main wing of the D8.5 design. 

However, considering the geometry of the flap relative to the wing under different deflection 
angles, a translation of the flap is required to match each deflection angle. This translation has a 
direct effect on the aerodynamics of the flap-wing combination and is therefore not arbitrary. An 
aero-mechanical trade study must be conducted to establish the near-optimum positioning of the 
flap. The focus of this report remains solely in the aerodynamic aspects of the gap size and ideal 
mechanical and motion constraints are established. A detailed discussion on this is presented 
further in the report. 

The aerodynamics of the D8.5 flaps configuration is conducted using the FlightStream™ surface 
vorticity numerical solver in conjunction with the open-source NASA Vehicle Sketch Pad (VSP) 
software. FlightStream™ is a proprietary software of the Research in Flight Company and uses the 
surface meshing capabilities of VSP to provide a fast, accurate solution to subsonic flow problems. 
It is capable of using structured and unstructured surface meshes and this capability allows 
effective integration with VSP. 
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2. EXPERIMENTAL VALIDATION 


Surfaced-vorticity based potential-flow methods have not been used extensively for flapped wing 
geometries. For this reason, the validity of FlightStream™ in this genre of aerodynamic analysis 
must be established. Research in Flight has identified some NACA, AIAA and NASA datasets 
believed to be appropriate for the requested validation portion of this study 14 . These test cases are 
summarized in Table-2.1. For each case, the test model is evaluated for an angle of attack sweep 
for both lift and drag conditions. Compressibility models are used for all three cases. However, 
considering the low Mach number regimes involved, the contribution of compressibility is minor. 


The validation process involved the creation of the geometry, surface meshing in either VSP or 
FlightStream™ and the collection and comparison of aerodynamic lift and drag data for an angle 
of attack sweep. Three different geometries were used in the validation process and are shown in 
Table- 1. 


Test case 

Table-2.1: Test cases for the validation effort 

Source 

Flap type 

Slats? 

Advanced General 
Aviation Wing 

NASA TN D-8236 

Semi-span / 
Fowler 

No 

NASA Trap Wing 

1 st AIAA CFD High Lift Prediction Workshop (HLPW1) 

Full- span / 
Fowler 

Yes / Full- 
span 

DLR-F11 

2 nd AIAA CFD High Lift Prediction Workshop (HLPW2) 

Full-span / 
Fowler 

Yes / Full- 
span 


2.1 The Advanced General Aviation Wing (AGAW) 

This geometry was obtained from the experimental work done for the NASA TN D-8236 report. 
The geometry information as well as flight conditions and aerodynamic data are extracted from 
this report and used for comparison with FlightStream™ results. The AGAW wing as tested in 
NASA TN D-8236, consisted of an asymmetric high-lift configuration and included spoilers on 
one side of the wing while employing flaps on both ends. The configuration was tested in the 
Langley V/STOL tunnel in support of the ATLIT aircraft program to determine the general 
characteristics of an ATLIT-type Fowler flap and spoiler lateral-control system and the 
experimental setup is described in NASA TN D-8236. 

The geometry is shown in Figure-2. 1.1 (a-c). The flap model used for this wing was a semi-span 
Fowler design and the airfoil for both the wing and flap are extracted from the basic GA(W)-1 
airfoil. The wing had a span of 4.01 m (13.16 ft.), a chord of 0.45 m (1.46 ft.), and an area of 1.79 
m 2 (19.31 ft 2 ). When the flaps were fully deflected, the wing area was increased by 17 percent to 
2.10 m 2 (22.59 ft 2 ). The wing root was at an incidence of 20 and the wing was linearly twisted to 
a tip incidence of 0°. The experimental data was corrected for tunnel wall effects. No other 
corrections were applied. Most of the data were obtained at a dynamic pressure of 1 .44 kPa (30 
lbf/ft 2 ). Additional data was obtained at a dynamic pressure of 0.48 kPa (10 lbf/ft 2 ). The Reynolds 
numbers corresponding to these dynamic pressures are 1.49xl0 6 and 0.85xl0 6 , respectively. 


4 





4.01 m 



Figure-2.1.1: The Advanced General Aviation Wing (AGAW) geometry and experimental layout 


The geometry and surface mesh used for the flow solution in FlightStream™ has been simplified 
from the experimental setup on account of unavailable experimental data on the mount, walls and 
mechanical joints and testing appendages on the wing. However, this simplification was found not 
to affect the solution in any significant manner. The surface mesh was generated in FlightStream™ 
using the extruder tool, which generates an aligned mesh suited for extrusion surfaces such as 
wings. One major simplification from the setup standpoint was to ignore the stowed spoilers as 
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being a geometric entity. Since the data was collected and compared with the spoilers stowed, this 
was not an egregious simplification and, as results show, did not significantly affect the 
aerodynamic results. However, an advantage of this design choice resulted in the treatment of the 
entire wing as symmetric as the spoilers were the only geometric entities that were causing 
asymmetry as seen Figure-2.1.1. This symmetric geometry takes advantage of the FlightStream™ 
symmetric boundary condition to significantly speed up the solver run time. The FlightStream™ 
surface mesh is shown in Figure-2.1.2 about the symmetry plane. The resulting surface mesh had 
2,470 triangle faces and 1,238 quadrilateral faces after solver initialization using the 
FlightStream™ quad mesher. 



Figure-2.1.2: The AGAW surface mesh as generated by the FlightStream™ extruder tool (left) and after 
solver initialization through the quad mesher (right) for a flap deflection of 10° 

The solver was run in parallel on 2 processors and required -45 seconds for each angle-of-attack 
flow convergence. Aerodynamic data in the deployed flap configuration for 10°, 20° and 30° 
deflections are compared with FlightStream™ results in Figure-2.1.3. The aerodynamic lift and 
drag coefficients for the angle of attack sweep from -4° to 10° are found to be in good agreement 
with the experimental values for all three deflections. 




Angl e of Attack ( D egrets) Angl e of Attack ( D egrees) 

Figure-2.1.3: Lift and Drag coefficients for the AGAW geometry for various flap deflection settings 


2.2 NASA Trap Wing Geometry 

The NASA Trap Wing geometry was obtained from the AIAA 1 st High Lift Prediction Workshop 
(HLPW-1), organized in 2010. The focus of this workshop is on the use of CFD to simulate the 
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aerodynamics of high-lift configurations. As such, the geometry and surface mesh data are 
provided. However, a surface mesh must be generated from the geometry which is available in 
CAD format. The NASA Trap-Wing as tested experimentally is shown in Figure-2.2.1. 



Figure-2.2.1: NASA Trap Wing Geometry in the wind tunnel (Image reference: AIAA HLPW-1) 

There are three configurations of the flap for the HLPW-1 analysis, as shown in Figure-2.2.2. The 
Configuration- 8 setup is used for this effort and modeling is done to include both the slats and 
flaps. For this configuration, the flap deflection angle is 20° and the slat deflection angle is 30°. 



Figure-2.2.2: Configurations for the NASA Trap Wing Geometry 


The free-stream conditions are near incompressible with a very low Mach number (M=0.2). The 
experimental angle-of-attack sweep used is from 0° to 37° at a Reynolds number of 4.3xl0 6 using 
the mean-aerodynamic-chord. The reference temperature is 520R. The experimental data is 
corrected for free-air conditions. The reference area is 22.028 ft 2 and the mean aerodynamic chord 
is 39.634 in. 
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Several simplification assumptions are made from this setup similar to the AGAW geometry 
analysis from section-2.1 of this report. Namely, the symmetry plane condition can be enforced to 
mimic the tunnel configuration as seen in Figure-2.2.1. Additionally, the exposed mechanical 
components of the flap (as seen in Figure-2.2.1) are not modeled in the FlightStream™ geometry 
and certain deviations from the solution are expected as a result. Meshing in FlightStream resulted 
in the use of the extruder for aligned meshes on the slat, wing and flap. An unstructured surface 
mesh was generated on the body pod and the floor of the tunnel was simulated through the use of 
the symmetry plane boundary. The resulting surface meshes are shown in Figure-2.2.3. The 
resulting surface mesh had 5,030 triangle faces and 3,062 quadrilateral faces after solver 
initialization using the FlightStream™ quad mesher (78.25% quad recovery). 



Figure-2.2.3: The NASA Trap-Wing surface mesh as generated by the FlightStream™ extruder tool (left) 
and after solver initialization through the quad mesher (right) for Configuration-8 


The solver was run in parallel on 2 processors and required -150 seconds for each angle-of-attack 
flow convergence. The Surface vorticity solution for 5° angle-of-attack is shown in Figure-2.2.4 
with the relaxed wake vorticity strands. Aerodynamic data for Configuration-8 are compared with 
FlightStream™ results in Figure-2.2.5(a-b). The numerical results are found to be in good 
agreement with the experimental values for all three deflections. 



Figure-2. 2.5(a) compares the lift coefficient with the angle of attack sweep. It is noted straight 
away that FlightStreamTM captures the correct linear slope of the curve. FlightStream™ has a 
small zero-offset with the experiment and over-predicts the lift slightly. These results are found to 
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be within acceptable limits for a potential environment solution. Similarly, the drag results shown 
in Figure-2. 2.5(b) compare favorably with experiment. The over-prediction in lift is carried 
forward into this figure as well and changes the shape of the lift-drag curve slightly. Even so, the 
results are found to be well within accepted error limits for a potential flow solver. Note that the 
angle of attack sweep for the numerical results were maintained within limits of the potential flow 
region of the geometry. Beyond 35° angle of attack, the flow is found to separate completely on 
the wing and the results cannot be evaluated using a potential flow solver. The separation and 
corresponding transition regions are clearly seen in Figure-2.2.5 beyond 30° angle of attack. 




Angle of Attack (Degrees) 


CD 


Figure-2.2.5(a-b): Numerical and Experimental aerodynamic results for the NASA Trap-Wing geometry 


2.3 DLR-F11 Geometry 

The DLR-F11 geometry was obtained from the AIAA 2 nd High Lift Prediction Workshop 
(HLPW-2), organized in 2013. The focus of this workshop is on the use of CFD to simulate the 
aerodynamics of typical commercial aviation aircraft in high-lift configurations. Similar to HLPW- 
1, the basic geometry data is provided. However, a surface mesh must be generated from the 
geometry which is available in CAD format. The DLR-F11 as tested experimentally is shown in 
Figure-2.3.1. The geometry includes both flaps and slats extended on a multi-section wing through 
the span. The configuration provided includes the slat at 26.5° and the flap at 32° deflections in 
their deployed state. 

The free-stream conditions are near incompressible with a very low Mach number (M=0. 175). The 
experimental angle-of- attack sweep used is from 0° to 25° at a Reynolds number of 15.1xl0 6 using 
the mean-aerodynamic-chord. The experimental data is corrected for free-air conditions. The 
reference area is 419.13 mm 2 and the mean aerodynamic chord is 347.09 mm. 
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Figure-2.3.1: The DLR-F11 geometry as imported into the FlightStream™ solver 


Once again, similar to the NASA Trap-Wing model, several simplification assumptions are made 
to the DLR-F11 simulation model. Namely, the symmetry plane condition can be enforced to 
mimic the tunnel configuration as seen in Figure-2.3.2. Additionally, the exposed mechanical 
components of the flap are not modeled in the FlightStream™ geometry and certain deviations 
from the solution are expected. 



Figure-2.3.2: The DLR Fll in the wind tunnel 
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Meshing in FlightStream resulted in the use of the extruder for aligned meshes on the slat, wing 
and flap. An unstructured surface mesh was generated on the fuselage of the DLR-F1 1 and the 
floor of the tunnel was simulated through the use of the symmetry plane boundary. The resulting 
quasi-valence surface mesh environment solved for in FlightStream™ are shown in Figure-2.3.3. 
The resulting surface mesh had 8,236 triangle faces and 5,355 quadrilateral faces after solver 
initialization using the FlightStream™ quad mesher (69.96% quad recovery). 




Figure-2.3.3: The quasi-valence surface mesh for the DLR-F11 as seen by the FlightStream™ solver 

The solver was run in parallel on 2 processors and required -230 seconds for each angle-of-attack 
flow convergence. Aerodynamic loads for the DLR-F1 1 are compared with FlightStream™ results 
in Figure-2. 3. 4(a-b). The numerical results are found to be in good agreement with the 
experimental values for all three deflections with two notable exceptions. The lift curve is found 
to become more accurate as the angle-of-attack increases, matching the zero-offset and slope of 
the lift curve almost exactly after 5° angle-of-attack. However, lift values at lower angles are found 
to be close to the experimental data but do not match the correct slope of the curve. It is speculated 
that the simplifications in the geometry lead to some variation in the results for the DLR-F11. 
Additionally, it is also possible that some non-potential physics may be at play which contribute 
to the deviations noted in Figure-2. 3.4(a). 
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Figure-2.3.5(a): Numerical and Experimental lift coefficient versus angle-of-attack for the DLR-F11 geometry 



CD 

Figure-2.3.5(b): Numerical and experiment lift versus drag results for the DLR-F11 geometry 
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3. BASELINE MODEL 


3.1 Geometry 

In order to incorporate a single-slotted flap, the basic D8.5 geometry was modified. The D8.5 was 
provided by NASA Langley to Research in Flight in the form of a VSP model (Figure-3.1.1). 
However, as the geometrical data required for modeling the details of the underside of the wing, 
the stow area for the flap and the flap geometry was supplied by NASA separately, certain 
modifications to the VSP D8.5 were made in FlightStream™. 



Figure-3.1.1: The initial D8.5 model without flaps 


Geometry changes in FlightStream™ included the use of the extruder to create extrusion sections 
between the major sections of the wing and flap. The single-slotted flap was incorporated in the 
wing between 10.3% of the span (measured from the centerline of the fuselage) to 79.4% of the 
span with a constant 28% of the local wing chord serving as the flap chord. The 79.4% span 
location also corresponded to the flap/aileron junction. The aileron was not modeled for this effort. 
The resulting flap data is summarized in Table-3.1 and Appendix A.2 of this report. 


Table-3.1: The simple-slotted flap geometry parameters for the D8.5 


Original NASA baseline flap design data 

Surface Type 

Clocal (Feet) 

Cf/Clocal 

b/2 

Description 

Flap 

11. OSS 

28.0% 

10.3% 

At fuselage edge where wing meets fuselage 

Flap 




At wing break 

Flap 

3.60455 

28.0% 

79.4% 

At flap/aileron junction 

Aileron 




At flap/aileron junction 

Aileron 




Attip of aileron 


The flap airfoil was evaluated for the nominal 28.0% local chord from the basic MIT NC-120 
airfoil used in the D8.5 wing. The chordwise refinement was based on the practices established for 
the validation studies previously described in this report. The refinement was further driven by the 
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need to have a suitable refinement in the opposing surfaces of the gap as shown in Figure-3.1.2 for 
the stowed condition. This nominal refinement was used in the extrusion surfaces to generate 
aligned meshes. The resulting surface meshes for the extruded flap and wing sections are shown 
in Figure-3.1.3. Note that FlightStream™ was successful in generating solutions even with very 
complex geometric features including the cavity in front of the flap. Repaneling this section of 
the lower surface of the wing to bridge the cavities was not required for the solver. 



Figure-3.1.2: Refinement nodes in the D8.5 flap mesh in stowed condition 



Figure-3.1.3(a): The D8.5 wing with flap in different positions during deployment 



Figure-3.1.3(b): Extrusion meshes generated for the D8.5 wing with flap 
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3.2 Flap Design 


The flap design as implemented in this study differs from an operational flap in practice in certain 
ways. The D8.5 main wing and flap are tapered. Hence rotation through compound angles are 
involved as the flap is simultaneously extended from a stowed position. Furthermore, rotation 
about the leading edge is required for the trade study. To accomplish this the desired positioning 
for each end profile was identified and a flap was extruded between the two end profiles. 

The motion of the flap has been further simplified by locating a fixed rotation axis on the leading 
edge within the flap once the two ends of the flap has been translated to its position for a given 
deflection angle. This rotation axis for the flap in the stowed condition is shown in Figure-3.2.1. 
Since the two ends of the flap translate independently, the flap encounters a negligible twist that 
would not occur in practice. 



FlightStream™ allows users to create user defined coordinate frames and to associate them with 
moving surfaces. This was utilized in this effort to mark the rotating axis of the flap in a local 
frame and to simply translate the frame as required. This decreases the simulation time for flap 
trade study by eliminating significant user inputs. 

3.3 Experimental Validation 

The baseline D8.5 with the stowed flaps was analyzed to ensure that no dramatic departures from 
the clean wing configuration would emerge in the analysis. The geometry was evaluated for 
aerodynamic loads for a range of angle of attack values as well as Mach numbers. The angle of 
attack sweep was from -3° to 4° and the corresponding Mach numbers were 0.2, 0.55, 0.6, 0.7, 
0.74, 0.77 and 0.8. Comparison was made with the Cart3D Euler CFD solver results for the same 
dataset and was provided by NASA Langley. A sample surface vorticity result for the D8.5 
baseline model with stowed flaps at 2° angle of attack is shown in Figure-3.3.1. The reference 
area for the geometry tested was 1162 ft 2 . The reference length was 170 ft. and mean chord was 
8.12 ft. 
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Figure-3.3.1: Vorticity distribution on the D8.5 with stowed flaps at 2° angle of attack 

A quantitative representation of the baseline solutions using the stowed flapped geometry is shown 
in Figure-3. 3.2(a-b). Generally good agreement with Cart 3D is seen to occur except at high Mach 
numbers for the lift as a function of angle of attack and this disagreement is mostly an offset rather 
than a slope issue. Because of this offset characteristic, the drag polar still agrees quite well with 
the Cart 3D results, indicating that the aerodynamic performance of the aircraft has been captured. 
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Figure-3.3.2(a): Lift versus angle of attack for the D8.5 baseline geometry 
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Figure-3.3.2(b): Lift versus drag for the D8.5 baseline geometry 
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4. FLAP MOTION DESIGN SPACE 


The results for this effort focus on the calculation of lift and drag coefficients for three different 
flap settings: 15°, 25°, and 35°. The trade study for the gap space is designed to take these three 
deflection angles into account. As the flap rotates through these angles, the translation design space 
must prevent the intersection of the flap with the lower surface of the wing. This criterion is critical 
in the selection of the design space. 

As previously described in Section-3 of this report, the translation motions of the extrusion flap 
surface was independent on the two ends of the flap. This allowed the creation of a unified, 
normalized design space as shown in Figure-4. 1 . The shape of the design space in the X-Z plane 
is such that the flap does not intersect with the wing at any position within it. This condition also 
determines the minimum offset needed between the upper bounds of the design space and the 
lower surface of the wing. 



Since the design space is unified for the two ends of the flap, it can be normalized using the local 
chord length to create the normalized design space shown in Figure-4.1. The translation space is 
restricted to the X-Z planes on either end of the flap. During the trade study, the flap ends are 
moved within this design space for each deflection angle to determine the optimum gap sizes. A 
collection of the modeled flaps within this design space for the stowed flaps case is shown in 
Figure-4.2. Contour plots for the trade study based on these discrete solution locations are shown 
in the next section of this report. 



Figure-4.2: Cascaded flaps at different positions within the design space for the stowed case 
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5. RESULTS 


The trade study involving looking at multiple flap locations within the established design space 
for each deflection angle of the flap. Table-5.1 shows the various locations that were each tested 
for the three flap deflection angles of 15°, 25°, and 35° and for each angle of attack sweep. Each 
flap is named on a two-digit number. The first digit specified its normalized location in the X 
direction and the second digit specifies the normalized location in the -Z direction. For example, 
Flap-24 is located 2 % of the local chord in the positive X direction compared to the baseline stowed 
position and 4% of the local chord in the negative Z direction within the design space. 


Table 5.1: Trad e stu dy flap position data within the normalized design space 


Flap offset [% chord) 

Actual Values (Root; Meters) 

Actual Values (Tip; Meters) 

X 

Z 

X 

Z 

X 

Z 

0 

-2 

0.000 

-0.068 

0.000 

-0.022 

0 

-4 

0.000 

-0.135 

0.000 

-0.044 

0 

-6 

0.000 

-0.203 

0.000 

-0.066 

2 

0 

0.068 

0.000 

0.022 

0.000 

2 

-2 

0.068 

-0.068 

0.022 

-0.022 

2 

-4 

0.068 

-0.135 

0.022 

-0.044 

4 

0 

0.135 

0.000 

0.044 

0.000 

4 

-2 

0.135 

-0.068 

0.044 

-0.022 

4 

-4 

0.135 

-0.135 

0.044 

-0.044 


The sweep of locations shown in Table-5.1 required 52 FlightStream™ solver runs for baseline 
compressibility validation sweeps and 276 FlightStream™ solver runs for the flaps trade study. 
Each run was conducted on parallel processors and required -180 seconds for each angle of attack. 
The overall trade study run time was therefore, -15 hours. 

Figure-5.1 shows the lift-versus induced drag distribution for Flap-02 at different deflection 
angles. The results are typical for a flap configuration on a wing with each lift-drag curve offset to 
the right as the deflection angle increases. 
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Figure-5.1: Lift versus induced-drag for Flap-02 in the design space for different deflections 
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Similar runs were completed for all of the flapped geometry cases and the results for nine flaps are 
shown summarily for an angle of attack sweep with a flap deflection angle of 15° in Figure-5. 2(a). 
Notice that, as expected the slopes are essentially identical for all of the flaps but the offset changes 
with flap position. The induced drag coefficient for the 15° angle setting is shown in Figure-5. 2(b) 
for the same set of flaps. Again the behavior is essentially identical for all of the flaps except for 
the offset depending on flap location. 
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Figure-5.2(a): Lift versus angle of attack for different flaps in the design space at 15° deflection 
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Figure-5.2(b): Induced drag versus angle of attack for different flaps in the design space at 15° deflection 


o 

Flap_02 

• 

Flap_04 

• 

Flap_06 

• 

Flap_20 

o 

Flap_22 

• 

Flap_24 

• 

Flap_40 

• 

Flap_42 

* 

Flap_44 




Q 0.03 








o 




* 

o 

• 


* * 

8 : 
• 

S 

o 

o 

o 


< 

o ° < 

o ° 

o 

i 








Q 

Flap_02 

• 

Flap_04 

o 

Flap_06 

■ 

Flap_20 

o 

Flap 22 

■ 

Flap_24 

* 

Flap_40 

■ 

Flap_42 

# 

Flap_44 






< 


[ 




c 

4 

? 

i 

i 

► 

! 




i 

i 

i 

i 

■ ! 

' 

P 

* 

i 



O 

* 

o 


i 

C 

> 



3 

0 

< 






o 






















20 


A summary of the lift and induced drag coefficients for the 25° deflection for the nine flap cases 
is shown in Figure-5. 3 (a-b). As expected the lift and induced drag are higher than for the lower 
angle settings but the character of the solution is unchanged. 
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Figure-5.3(a): Lift versus angle of attack for different flaps in the design space at 25° deflection 
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Figure-5.3(b): Induced drag versus angle of attack for different flaps in the design space at 25° deflection 

A still further increase in the lift and induced drag is associated with the higher flap setting of 35° 
as shown in Figure-5. 4(a-b). Again, we see linear aerodynamics at work for lift coefficient and a 
nonlinear induced drag with angle of attack as expected. 
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Figure-5.4(a): Lift versus angle of attack for different flaps in the design space at 35° deflection 
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Figure-5.4(b): Induced drag versus angle of attack for different flaps in the design space at 35° deflection 


The drag polars shown in Figure-5. 5(a-c) for 15°, 25° and 35° offer some perspective on how lift 
clearly increases at increased deflection angles and increases to the highest values for the highest 
combined flap deflection angles and angles of attack. This compound effect increases induced drag 
as expected and as represented by the fact that the curves move to the right as angle of attack is 
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increased successively from (a) to (c), all as expected. This data is a quantification of how this 
expected result is delivered by the D8.5 geometry. 
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Figure-5.5(a): Lift versus induced drag for different flaps in the design space at 15° deflection 
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Figure-5.5(b): Lift versus induced drag for different flaps in the design space at 25° deflection 
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Figure-5.5(c): Lift versus induced drag for different flaps in the design space at 35° deflection 

Figure-5. 6(a-c) is another way to package the data from a pseudo aerodynamic efficiency 
perspective. In many cases, maximum lift to drag does not necessarily lead to the globally optimal 
solution for a given design when other factors such as mass, takeoff distance, flutter or other 
performance factors are considered. Nevertheless, this is an appropriate way to present the results 
of this purely aerodynamic performance based trade study. Since all cases in the trade study 
involve significant flap deflection angles, the lift is positive for all angles of attack and the peak 
aerodynamic efficiency is near zero angle of attack. 


a 

AJ 30 


4 

• 

► i 

4 

► 




4 

* 

4 

> 4 

1 

4 

> 

4 

4 

► 

if 

i 




: * 

: ' 

* 

1 

1 ! 

4 

1 1 

r i 

I 

i 

i 

4 

1 

1 

1 

* 























Flap_02 

Fiap_04 

Flap_06 

Flap_20 

Flap_22 

Fiap_24 

Flap_40 

Fiap_42 

Flap_44 


Figure-5.6(a): Lift-to-drag ratio versus angle of attack for different flaps in the design space at 15° deflection 
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Figure-5.6(b): Lift-to-drag ratio versus angle of attack for different flaps in the design space at 25° deflection 
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Figure-5.6(c): Lift-to-drag ratio versus angle of attack for different flaps in the design space at 35° deflection 

The plots in Figure-7 are contour plots of the performance of lift (a), induced drag (b), and lift to 
induced drag (c). The geometry is planar and represents the x-z location of the leading edge 
location of the flap. The upper left corner is cut out of this two-dimensional plot contour plot 
because this would represent interference between the wing and the flap. Interestingly, the 
aerodynamic efficiency represented in Figure-7. 1(c) is a maximum with the flap at the minimum 
separation from the main wing. This is mainly due to the lower induced drag at these locations. 
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Keep in mind that this does not account for viscous affects. The lower flap deflection angles 
correspond to the higher efficiencies. Despite the obvious best values as represented in the contour 
plot, the variation in efficiency over a significant flap position range is not dramatic and it is 
entirely possible that the best design point will be driven by other factors in practice. 
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Figure-5.7(a): Lift distribution in the normalized design space for different flap deflections 
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Figure-5.7(b): Induced drag distribution in the normalized design space for different flap deflections 
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Figure-5.7(c): Lift-to-induced drag distribution in the normalized design space for different flap deflections 
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6. CONCLUSIONS 


A fundamentally new solver, FlightStream™, based entirely on surface vorticity has been used to 
investigate the NASA D8.5 geometry in a high lift (flapped) configuration. This investigation 
included a solution for a baseline configuration with no flaps, a solution with a flapped 
configuration with the flaps stowed, a validation study to ensure that the flapped configuration 
solutions would be reliable, and a series of solutions at flap angle settings of 15°, 25°, and 35° for 
a range of flap positions and angles of attack. The validation cases included comparisons to a 
flapped general aviation configuration, the first high lift workshop configuration with flaps only 
and the second high lift workshop configuration with both slats and flaps. The fidelity of the 
FlightStream™ solution with the available experimental data was quite striking for angles of attack 
up to a range at which flow separation is thought to be a likely component of the flow physics. 
FlightStream™ is a potential flow solver based on the assumption that the flow is entirely attached 
and therefore is not expected to maintain high fidelity with wind tunnel or flight data in which 
separated flow is present or anticipated. 

The D8.5 flap study showed that a near doubling of lift coefficient is possible with the flaps that 
were modeled. The trade study included a map of X-Z locations for the flap leading edge and this 
study showed that the best efficiency occurs for flaps that are in very close proximity to the wing. 
On the other hand, this does not correspond to the highest lift configuration. The highest lift 
configuration corresponds to the most fully extended position for the flaps, which also corresponds 
to the highest induced drag position, as expected. Depending on the design goals for the entire 
aircraft, different flaps, flap positions and possible flap settings might be of interest. 

This study has been designed to create an aerodynamic map of the D8.5 geometry, specifically 
with flapped wings for the purposes of both validating the method and generating an aerodynamic 
performance map suitable for use by an aircraft design team. This analysis does not address the 
structural requirements for the aerodynamic loading, flutter considerations, propulsion 
requirements, or any of a number of other design constraints and focuses instead entirely on 
providing a high fidelity characterization of the aerodynamic performance. 

It is highly recommended that additional work on this and similar geometries include the 
development of methods suitable for accurately estimating the value of CL m ax, which would 
involve the development of appropriate methods for analyzing separated flow over wings. The 
Research in Flight Company is very interested in pursuing methodologies in this genre. 
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APPENDIX 


A.l RESEARCH IN FLIGHT 

The Research In Flight (RIF) Company was established in Auburn, Alabama in 201 1 with the aim 
of developing new tools for Aerospace Engineering applications. Since its inception, RIF has 
developed numerous tools for design and optimization of aerospace vehicles which are at the 
considered state-of-the-art in multi-disciplinary aero-propulsive optimizations of preliminary 
designs. Each of these tools allows complete vehicle optimization using powerful bio-mimetic 
inspired optimization tools also developed by RIF. 

Research in Flight’s flagship software, FlightStream™, has been validated with both extensive in- 
house testing and in joint collaborative tests with NASA Langley. With tests confirming 
Flights tream™’s ability to generate accurate aerodynamic results much faster than equivalent 
Navier-Stokes solvers, it is now being used as the aerodynamics toolbox for a variety of uses 
ranging from generating load-maps of flap and control-surface deflections, multidisciplinary 
optimization of entire airframe outer-mold-lines, investigation of new vorticity phenomena, stores 
separation problems, wake-fuselage interactions, formation flying analysis of unmanned vehicles, 
missile control dynamics, wind turbine design, propellers, helicopter aerodynamics and much 
more. 

Research in Flight offers a range of aerodynamic analysis services primarily centered around the 
use of FlightStream™ to produce high-fidelity, rapid aerodynamic load analyses for flight vehicles 
in the subsonic regime. Accurate analysis of aero vehicle wakes, wake-related phenomena and 
close proximity configurations is available. 

For more information on RIF, FlightStream™, white-papers and validation efforts, visit the 
company website at www.researchinflight.com 
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